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ON EVALUATING THE QUALITY OF NONLINEAR 
AUTOMATIC SYSTEMS WITH RANDOM 
INTERFERENCE 


E. P. Popov 
(Leningrad) 


The method of statistical linearization [1, 2] and its further development from the aspect of harmonic 
linearization for self- oscillating systems [3, 4] are the mast effective means for an approximate investigation 
of the dynamics of automatic systems with typical nonlinearities (relay, saturation, insensitivity zone, cubic 
dependence, etc), Some other methods are described in [5]. However, the solution of the given problem, as 
was shown there, is still in an initial stage of development, 


A new application of a combination of the methods of statistical and harmonic linearization to the problem 
of investigating the influence of random interference on the stability and dynamic properties of nonlinear auto- 
matic control systems is given in this article, It is known, for example, that a decrease of the amplification 
factor of the useful signal under the influence of noise is observed in the presence of saturation nonlinearity (see 
(6), page 278). It is evident that if such an amplifier is included in a closed automatic system, the latter can 
substantially change its dynamic properties in relation to the useful signal under the influence of random inter- 
ference up to the point of possibie loss of stability, ultimately, In connection with this, an attempt to develop 
an approximate theory of the above-mentioned phenomena is undertaken here, and is illustrated with an example 
of the calculation of a concrete automatic system. 


1. Initial Equations 





Let a nonlinear automatic control system, having an arbitrary structure, be described by the differential 
equation 


Q(p)x + R(p) F (x) = Sy (p) f1 (t) + Sa (P) /2 ), (1.1) 


where F (x) is the given nonlinearity, f, (t) is random interference with zero mathematical expectation, f, (t) is 
a regular extemal perturbation, which slowly changes in comparison with the interference, and Q (p), R (p), 

S; (p) and S, (p) are operators with constant coefficients (p = d/ dt), whose form depends on the different distri- 
butions of the variable x and the external perturbations f, (t) and f, (t). 


Let us seek a solution in the form 


x = tu (t) + 2p (2), (1.2) 


where x,, (t) is the regular component or mathematical expectation (useful signal), Xp (t) is the central random 
component (interference at the input to the nonlinearity), 


Considering the interference f, (t) as a stationary random process, given by its spectral density s¢(w), we 
will also consider the solution for the component xp (t) as a stationary random process, This supposition is based 
on the slow change of the mathematical expectation x,, (t), in so far as it can be considered constant for possible 
"periods" of change of xp (t) corresponding to the given frequency spectrum s¢(w); that is, the case is examined 
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when the interference frequencies are high in comparison to the rate at which the useful signal can be changed, 


By using, as in the case of harmonic linearization, the property of a filter in the linear part of the system, 
we will approximately characterize the random processes in the nonlinear system (1.1) by the two first statistical 
moments; the mathematical expectation (mean value) and the dispersion (or mean square deviation), In 
correspondance with this, we will produce the following statistical linearization of the nonlinearity: 


F (t)= Fut Qp2p, (1.3) 


where F,,, is the regular component (mathematical expectation) of the nonlinear function,gp is the amplification 
factor of the random component in the nonlinear link, where for single- valued nonlinearities 


+00 


Fy = M[F (1% + 2p) = \ F (2 + 2p) w(x) dx 


for the normal distribution law 


w (2) = — ho, a9.“ 





where oy is the mean square value of the random component xp, We will select the simpler second method [1] 
for determination of qr: 


+o 


IR =; \ F (ry + 2p) rw (x) dz. 


—co 


Analogous formulas can also be written for loop nonlinearities, 


As a result, we have ready expressions for any typical nonlinearity 


Fy (Zu, Sx), (1.4) 
IR (Zu, Sx), (1.5) 
which are depicted graphically, as is shown, for example, in Fig. 1, for saturation nonlinearity [1]. 
Putting (1,2) and (1.3) in the given equation of the nonlinear system (1,1), we will separate it into two 
equations: 


Q(p) tu + R(p) Pu = Si (PAO, (1.6) 
Q(p)t% + R(P)g TR = S2(P) fo (2), (1.7) 
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respectively for the regular component (mathematical expectation), 


which is the useful signal, and for the random component (central), 
which corresponds to interference. 


The solutions of Eqs, (1.6) and (1.7) are associated with one 
another in a nonlinear manner, according to (1,4) and (1.5); this 
is essential for the ultimate manifestation of the nonlinear pro- 
perties of the given system, 











. 
Fig. 2. On the basis of the equation for the random component (1.7), 
and knowing the spectral density of the external interference s¢(w), 
it is possible to determine the dispersion of the interference at the input of the nonlinear link a? according to the 
known formula . 





-Loo 
‘ 1 S (jo) » 
Dee ae 
ox 2n \ Q (jo) + R (jo) g R (Tys Fy) sf (w) dw 
—oo 
or 
a2 = Alyn (Xu, Sx), (1.8) 


where a is a constant factor (including the ordinarily given quantity o ¢) placed outside the integral sign, and 


+00 
1 G (w) 
In= ae \ ray = 


—oo 


in which G (w) includes the square of the modulus of S (jw) and the part of the numerator of $¢(w), which depends 
on w, while |H (jw)|* is the whole denominator of the integrand; the denominator of s¢(w), transformed in an 
appropriate manner, is also included in the term | H (jw) 2. In the literature (see for example, [2], page 643) 
there are prepared expressions for 1, up ton = 7, where n is the degree of the polynomial H (jw). 


From Eq. (1.8), it is possible to find the mean square value of the interference 0, at the input of the non- 
linear link, as a function of the magnitude of the useful signal 
Sx (Tm) (1,10) 
graphically, in a manner similar to the solution for A (x*) in [7]. In order to do this, we write 


o2=C, aln(tu, 2x) =% 


instead of (1.8), and we construct a parabola 1 (Fig. 2) on the plane €, oy from the first equation, and a series 
of curve segments 0 , (€) for various x,, = const from the second equation, Their intersections with curve 1, when 
translated to the plane Xm, 0, (curve 2 in Fig, 2) give the solution sought for (1,10), 


2. Investigation of the Properties of a Nonlinear System in Relation to the Useful 





Signal in the Presence of Random Interference 





Having determined the dependence o , (x,,,), in such a manner we will place it in expression (1,4); asa 
result, 0, is eliminated, and we get the function 


Fy = D (zy), (2,1) 
which can serve as a characteristic for determining the passage of the useful signal through the nonlinear link in 


the presence of random interference, Having placed it in (1,6), we obtain the nonlinear equation for the useful 
signal of the system as a whole in the form 


Q (p) tu + FR (p) D (ru) = Si (PA), (2.2) 
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which differs from the true Eq. (1,1) in that here the interference is discarded, but a new function © (xp), 
which was found for the useful signal and accounts for the effect of interference, is substituted for the given non- 
linear ity F (x). 


It appears that the function © (x,,) ordinarily has the shape of a smooth curve even for a step~ function 
nonlinearity F (x) (relay, loop, etc)— the effect is that of smoothing the nonlinearity under the influence of 
random interference, Therefore it is possible to subject it, in many cases, to the ordinary linearization 


D (2) = kurt, (2.3) 


where ky is determined either graphically as the slope of the curve (xj), or analytically: 


D 
hols -) (2.4) 
Then the nonlinear equation for the useful signal of the system (2,2) is replaced by the linear equation: 


[2 (p) + kak (p)) am = Si (p) fi (0), (2.5) 


which significantly simplifies all further investigation of the dynamic properties of the system in relation to the 
useful signal (in the presence of random interference), since any known methods of linear automatic control theory 
can be used to this end, 


However, in spite of the use of linear methods, all the essential nonlinear properties of the system are 
manifested, due to the particular method, presented above, of determining the amplification factor of the useful 
signal in the nonlinear link kay in the presence of random interference, Indeed, the quantity ky, like the function 
© (Xp), is not simply a given number distinct from the amplification factors of ordinary linear links, According 
to (1.4) and (1.8), the quantities ky; and @ (xm) depend, firstly, on the spectral density of the external random 
interference s¢(w), on the basis of its mean square level o ¢ and frequency spectrum; secondly, on the shape of 
the nonlinearity F(x); and thirdly, on the structure and on ail other parameters of the system itself. 


All the statistical and dynamic properties of the system, which were determined from Eq, (2.5), including 
the system stability, will depend on these factors. Their dependence on the level of the extemal interference is 
an especially important distinction from ordinary linear systems (see example below), 


The dependence of the quantity k;; on the structure and on all the parameters of the system (manifested 
through the magnitude of 0, from which kj; is determined) is very important to consider in synthesizing the 
system, that is, in the choice of the structure and parameters of the system, according to Eq, (2,5), for desirable 
dynamic properties, If, for example, the magnitude of any parameter, k, of the system is selected, then it is 
necessary preliminarily to determine the dependence ky; (k), which is not difficult to do in the method for finding 
ky described above (see example below). 


The method of determining the quantity kj; can be significantly simplified in the particular case when the 
linear part of the system, determined by the transfer function 


R (p) 
W = 
does not pass the frequency spectrum corresponding to the given spectral density of the external interference 

sf (w), that is, it is considered in the particular given case that R (jw)/Q (jw) » 0 is in the limits of the spectrum 
of the random component xp (t), Then, dividing Eq, (1.7) by Q (p) and considering the above-mentioned pro- 
perty, we obtain a simple formula for the dispersion of the interference at the input of the nonlinear link: 


-oO 


— S (jo) |? 
=>, \ Q (ja) 


—o 





Sf (@) dw, (2,7) 


which, unlike (1,8), does not depend on the magnitude of the useful signal xyp. 


1284 








k, According to (2.1) and (1,4), the function (x 
| will concur with F,, (Xm, 0x) in this case, that is, 

it will be determined directly from a graph of the 
type shown in Fig. 1a, at an appropriate value of 
Oo, (2,7). Consequently, the amplification factor of 
the useful signal in the nonlinear link ky; also will be 


m) 








0 & determined simply as the tangent of the slope angle 
x of the appropriate curve of Fy, (xq) in Fig. 1a, for 
Fig. 3. a given 0, or 
a b OF 
: Instability ny 






k, =(— (2.8) 
Instability - (ir, sive 











Stability stiles directly from expression (1.4), 
In other systems it is often observed that al- 
though the linear part (2,6) does pass the frequency 
0 Z 9 3 spectrum s¢(w), it is not passed by some other separate 
Fig, 4 block (most often, the object being controlled), Then 


a simpler nonlinear equation of the system without 
this block is formed for the interference, instead of 
(1.7). In such a case, one is obliged to use the general method of determining 0, according to (1.8), but for a 
simpler equation (lower order), 


In both of the particular cases described, which are often encountered, the simplification of the solution 
is concerned only with the determination of how @ (x,,,) or ky; depends on the interference, after which the com- 


plete equation of the system (2,2) or (2.5) is used in all cases for an investigation of the properties of the system 
in relation to the useful signal, 


If the linearization of the function @ (x,,,) in the form (2.3) is undesirable for some reason and it is 
necessary to evaluate the effect of the virtual curvature of the function © (x;), one is obliged to investigate 
the nonlinear equaticn (2,2) directly with any of the known nonlinear methods of control theory, In particular, 
it is possible to find the stability of the system from Eq, (2,2) by the method of harmonic linearization or har- 
monic balance [8, 9], and also to find the self- and forced-oscillations in relation to the useful signal (10, 11), 
and to evaluate the transients [12, 13]. In order to do this, it is essential to produce a harmonic linearization of 
the function @ (x,,,) for Xm =A sin # in the form 


® (zu) = 9 lA: Ox) Du, 


where 


an 


a ( F,(Asing, ox) sing dd. (2.9) 
_0 


Here © (x;,,) is replaced by the expression Fin (%ms Ox) since it is necessary that the linear part (2,6) does 
not pass the higher harmonics of the useful signal, in order to use harmonic linearization; this means that the 


whole interference spectrum must not be passed, * Therefore, the quantity o, in expression (2,9) is determined 
from the simple formula (2.7), 


As a result, we get the following harmonically linearized equation for investigating the properties of the 
system in relation to the useful signal, instead of (2,2): 
Q(p) +94, (A, 2x) R(p) tu = Sy (P) A) (2,10) 
at the particular values f, = 0 and f, = B sin Qt. 


*In those cases when the interference frequencies are higher than the possible self-oscillation frequencies in 
relation to the useful signal. 
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Ready expressions for the coefficient of harmonic linearization q (A, 0x), obtained after statistical lineari- 
zation which was done earlier, were given in [3] for different types of nonlinearities. They now can be used for 
the calculation of a system in relation to the useful signal in the presence of interference according to Eq, (2.10) 
with exactly the same methods as in [8-13], but with the use of a new expression 4, (A. 0 x) instead of q (A); 

Oy is determined from formula (2.7). 


3. Example 











A nonlinear automatic system is described by the equation 
Pp? (Tsp + 1)2+ (kak oP" + Tykkyp + kk) F (x) = kp® (Tip + 1) fe (t) 


with a saturation nonlinearity F (x), for which the expressions Fr, (xm, 0x) and ap (Xe Oy) are shown graphi- 
cally in Fig. 1, The random external interference f, (t) is given by the spectral density 


2B oF 


<p lee 22)? 2p? 
(w; — a?w*)? + p*o 





The parameters of the system and the spectral density are such that the transfer function (2.6) 


kak a + Tokkyp + kk, 
p® (Tap + 1) 





does not pass the frequency spectrum s¢(w), It is necessary to determine the stability of the system in relation 
to the parameters k and T,, depending on the interference level of, 


We find geometrically from Fig, 1a the quantity kj; in the function from 0, (Fig. 3). According to (2.5), 
we write the characteristic equation of the system in relation to the useful signal 


pP® (Tap + 1) + (kakog p+ T,kkyp + kk) ky (o,) = 0, 
the condition of the system stability according to the Hurwitz criterion will be 


i 














ky (o,) > Ts (kek oa T ykkiT 1) ; (3,1) 
We have from formula (2.7) 
ae ik (Tyjo + 4) |? 28a? 
ot. \ Tw + 1 (w? — 2%)? + 2%? ’ 
whence 
So, = ko, V 28s, (3.2) 
where 
4 Ly (T?0* + 1) da 
fom Be » | a®7's (jo)® + (a* + 17's) (jo)* + (@?7 2 + p) jo + oF? 


a® + »7',— Tio? 
~ Qpw? (a? + pT's + T2w?) © 





Using the stability condition (3,1), the expression for 0 , (3.2), and the dependence k;, (0 ,) (Fig. 3), we 
can plot graphs of the system stability in relation to any desired parameters, depending on the interference level. 
In particular, the stability limits have the form shown in Fig, 4a and b, in relation to the parameters T, and k. 
Here the essential dependence of the limits of system stability (in relation to the useful signal) on the random 
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interference level can be seen; this is not observed in linear systems, The possibility of loss of stability at a 
definite interference level, or, in any case, the essential dependence of all statistical and dynamic properties of 
the system on the external interference is an important peculiarity, which must be taken into account when 
nonlinear automatic control systems are designed, 


SUMMARY 


The paper deals with a new application of a combination of statistic and harmonic linearizing methods 
to the investigation of the random disturbance effect on stability and dynamic properti2s of nonlinear auto- 
matic control systems. In this connection there was made an attempt to develop the rough theory of these 
phenomena with a control system design given to illustrate it. 
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L. Ss PONTRYAGIN MAXIMUM PRINCIPLE IN THE 
THEORY OF OPTIMUM SYSTEMS, I 


L. I. Rozonoér 


(Moscow ) 


Questions are discussed, which are associated with the proof and use of the 
L. S, Pontryagin maximum principle in the theory of optimum systems. The work 
also contains some new results, The problem of optimization for the case of a free 
right end of a trajectory is examined in the first part of the work, In the second 
part, the maximum principle is formulated for boundary conditions of a more 
general type, The connection between the method of dynamic programing and the 
maximum principle is established in the third part; a method of solving optimi~- 
zation problems in linear discrete systems is given, and a number of considerations 
conceming the use of the maximum principle in the solution of a definite class of 
problems associated with the theory of dynamic accuracy of control systems are 
also discussed, 


The problem of creating systems, which are optimum in some prescribed sense, is one of the most im- 
portant problems of automatic control quality theory. Questions about obtaining optimum laws of control were 
first seriously posed, apparently, in the theory of the motion of a rocket (see, for example, [17-19]). However, 
the results obtained here are standard from the viewpoint of general control theory. 


The majority of the most important and interesting practical problems lead to nonclassical variation 
problems, which require new mathematical methods for their solution, The circumstance comes from the fact 
that the variables in the equation of a real system are (or should be) limited in some way or other, The first 
most general and important results in this sphere concern theories of systems which insure minimum time of 
control processes, In 1952-1955, the bases of the theory of processes in linear systems, optimum in relation to 
rapid action, were developed, mainly in the works of A. A. Feld*baum and A, Ya. Lerner. In particular, the 
solution of the problem of synthesizing systems, optimum in relation to time, for the case of a linear aperiodic 
object of arbitrary order with one controlling perturbation, was given in [1-2]. A detailed review of results 
obtained up to 1956 can be found in [3-5]. 


In 1956, a principle, leading to the solution of the general problem of finding a control process, optimum 
for rapid action, was hypothesized by L. S. Pontryagin on the basis of the results of work performed by him, 
V. G, Boltyanskii, and R, V. Gamkrelidze [6-7]. This principle, which received the name “Maximum Principle," 
was verified at first for individual types of systems, and in particular, was proved in [8] for the case of linear 
systems, V. G. Boltyanskii [9] fully proved that the Pontryagin maximum principle was a necessary condition 
for optimality in relation to rapid action, R. V. Gamkrelidze [8-10] proved theorems of existence and uniqueness 
and examined the problem of synthesizing controls for linear systems, optimum for rapid action, In [11], the 
maximum principle was extended to the general case of minimizing an arbitrary functional of the integral 
function of variable systems. In [12], there is a detailed presentation of basic results obtained by L. S. Pontryagin 
and his associates, 


In examining questions concerned with the theory of optimum systems, it is necessary to note the numerous 
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works of R, Bellman, which are systematically presented in [13] (see also [14]), The method of “Dynamic Pro- 


gramming,” developed by R, Bellman, gives a new tool for the solution of nonclassical variational problems, 
which are closely associated with the Pontryagin maximum principle, 


The works of N. N, Krasovskii [15, 16] are also concerned with general questions on the theory of optimum 
systems; discrete systems, described by finite difference equations, are examined in [16]. 


The short review presented above is concerned only with papers connected with the maximum principle 
in some way or other, Hence, many interesting papers were not mentioned, 


This paper contains a presentation of the most essential problems in the theory of automatic control which 
aré associated with the proof and use of the Pontryagin maximum principle, The presentation of the majority of 
the known results (in particular, the proof of the maximum principle) differs somewhat from that in the original 
works, References to the original works are not given in the text as a rule, so as not to complicate the pre- 
sentation, Some new results are also contained in this work; a formula for the increment of a functional is 
based on these results; and also, the sufficient conditions of optimality, formulation of boundary conditions for 
a number of problems, the connection between the maximum principle and the method of dynamic programming, 
the formulation and proof of a principle for optimum processes in linear discrete systems, analogous to the maxi- 
mum principle, result from this formula, The work is published in the form of three separate articles, The 
optimization problem for the case of a free right end of a trajectory is examined in the first part, In the second 
part, the maximum principle is formulated for the general case of arbitrary boundary conditions, and the problem 
of synthesizing optimum systems is discussed. In the third part, the connection between the maximum principle 
and the method of dynamic programming is established; a method of solving some problems of optimization in 
discrete systems is given; also, a number of considerations concerning the use of the maximum principle in the 
solution of a definite class of problems connected with the theory of dynamic precision of control systems are 
presented (in particular, B, V. Bulgakov's problems about the accumulation of deviations for nonlinear systems). 


The most complicated proofs have been placed in appropriate appendices, 


1. Statement of the Problem 





Let us examine an object of automatic control, which has r controlling elements and is described by a 
system of differential equations of the n-th order: 


ay = f(t, «-->Tni Uy. -+sUrs b) (i=4,...,n), (1) 


where x,..., X, are the parameters of the object, and Uy,..., Uy are the positions of the controlling elements, At 
each moment, the positions of the controlling elements uy (t),... U, (t) must satisfy the inequalities 


95 (Ur (t),---,Ur(t)) SO = 4,-- +. m), (2) 


which reflect the restrictions which are imposed on the control system, In the future, we will use a geometrical 
terminology, and we will consider an r-fold space R (uy,..«, Ur). Some point (vector) u of the space R corresponds 
to any given combination of the positions of the controlling elements (Ugy..e, Ur), The vector U (t) = (tg (t)seeestly(t)s 
given as a time function, will be called the “control” of the object, The control u (t) will be considered piece- 
wise continuous, that is, each of the functions u;,, (t) can have a finite number of discontinuities of the first type 

at the end of an interval of time, 


The inequalities (2) isolate some closed*set of points, whose boundaries are given by the equalities 
9; (Ugr-+-» Ur) = 0, in the space R, This set will be designed by U, and the inequalitities (2) can be briefly written 
in the form 


u(t)eU, (2°) 
that is, the vector u(t) must belong to the closed set U of the space R at any given moment. In this case, our 


*The set is called closed if it includes its own boundary. The set being examined is closed in so far as con- 
ditions (2) permit a sign of the equalities such that the point u can be found on the boundary of the set, 
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choice of the possible controls is limited by condition (2°) and by the requirement of piecewise continuity of 
the functions uy (t),..., Ur (t). The control u(t), which obeys these two conditions, will be called admissible. 
When speaking about controls from now on, we will have admissible controls in mind, 


We will examine an n-fold phase space X (x,,..., X,) of the parameters of the object, along with the space 
R (Ug see0, Uy), in which the control can be changed, The state of the object is described then by the point (vector) 
X (Xqeee09 Xp) Of the space X, If the control u (t) and the initial state of the system x’ = oe Xn) are given, then 
the behavior of the system (trajectory x (t) in the phase space X) is uniquely determined from Eqs, (1), 


The basic problem of the theory of optimum systems is the choice of a control u (t) such that the behavior 
of the system will be “best” in some prescribed sense, We will introduce a number of examples of precise state- 
ments of problems, in which the sense of the concept “optimum control" is made concrete each time, 


It is possible to require that the control u (t) be so selected that the system will go from one fixed state x* 
to another fixed state x’ in minimum time. The control u (t), which satisfies this requirement, is called opti- 
mum for rapid action. The time of transition of the system from state x" to state x', which determines the 
criteria of optimality in the given problem, is a functional* given on the controls u(t). The requirement of 
optimality of the control u (t) means minimization of the value of the functional, 


Let us examine another problem, Let it be required that a control u(t) be chosen such that in a given 
time T the system will go from an initial state x° into a state in which one of the variables (for example, x) 
will become as large as possible, and the remaining coordinates will be fixed values, This problem arises, for 
example, in the calculation of the control law of a rocket, or in putting an artificial satellite of the Earth into 
a given orbit [20], if it is required that the horizontal velocity of the rocket be maximum at a fixed height and 
at a given time, for zero vertical velocity and for a given flight range. In the present example, the functional, 
which determines the criteria of optimality, is the value of the generalized coordinate x, (T) (horizontal velocity 
of the rocket) at the time t = T. In a number of cases, there is no necessity to fix the final values of all the re~ 
maining coordinates x,,..., X,. If, for example, only the height and not the flight range plays a role, there is no 
necessity to impose any restrictions on the latter, Then an analogous problem about the maximum of the func- 
tional x, (T) arises, but it is not necessary for all the remaining coordinates to assume prescribed values, Let us 
examine, finally, one more problem, It is necessary to choose a control u (t) in such a manner that the integral 


T 
Sa \F (x1,..., 20 iy, -.., Uy) at 


0 


becomes minimum (or maximum) during the time of movement T. Here again, different variations are possible, 
It can be required that the initial and final states of the system be fixed beforehand (for example, in the problem 
about the flight of an airplane from one fixed point to another, with the least fuel consumption; the function 

F (x, u) represents fuel consumed per unit time, and the integral represents the full consumption of fuel during 
the flight), In other cases, the final state of the system does not have any value; it is only important that the 
integral become the maximum (or minimum) value of all those possible (for example, if it is required that the 
aggregate being controlled give the maximum quantity of production in the absence of any other conditions; in 
this case, the function F (x, u) represents instantaneous producitivity, and the integral represents the full pro- 
duction output), The time T, during which optimization of the system is required, can also be either fixed be- 
forehand, or not fixed, If the airplane moves between two points on a schedule, then its flight time is fixed, and 
minimization of fuel consumption must be required for this time, fixed beforehand, If time does not matter, but 
only fuel economy is important, then the time must not be fixed beforehand. 


The number of examples of such optimization problems can be increased further, However, all these 
superficially different problems prove to be closely associated mathematically, and require completely analogous 
methods of investigation, Moreover, one of them can quickly to transformed into others, 


*A functional is a correspondence established between functions (or combinations of functions) and numbers, In 
the given example, a number— the transition time from x’ to x'- corresponds to each control u (t) [that is, the 
combination of functions u, (t),..., u, (t)) that insures that the system gets to x’, 
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We will show, for example, that the problem of optimizing a system with respect to an integral leads to 
the problem of optimizing with respect to coordinates, Actually, let it be necessary to optimize the magnitude 


T 


\ F (2, ..+, Un; Uy,..., Up; t)dt (the time T, and also the final values of the coordinates can be either 
0 


fixed, or free) in an object described by the differential equations (1), with restrictions on the control (2), We 
will introduce a new variable 


t 
In41 ()=(F i, ..., 20s Uy,.++, Up; t)dt, In4, (0) =2r4, = O. 
0 


Then one more differential relation can be added to system (1) 
Int = F(x, «++, Tn; Uy,--+, Up; b), 


and the problem of optimizing the integral leads to the problem of optimizing the n + 1 st coordinate x, ,4 (T) 
at the final moment of time, 


The problem about the optimization of some function from the final value of the coordinates © [xq (T)s.«-, 
Xp (T)] also leads easily to the problem about optimizing the final value of one coordinate, if the function @ is 
differentiable, 


To show this, it is sufficient to introduce a new coordinate, 


Tn+1 (t) = O[z, (¢), » an. Zp (t)], In+1 (0) =20,, 7 ® (2°, °r ee x) 
and to add the relation to system (1) 
- 8® (x, eoesg 


° x.) 
Inti = >) a, — fol, «+5 Snr Bayeecy Mes Spe 
8=1 





Let us further note that the problem of optimization with “integral” restrictions of the type 


T 
| F (ey. .-, oni Uy,-.-, Ur, thdt <A 
0 


easily leads to ordinary statements of the problem, Introducing the variable xp 4, (t) and adding the relation 
Int, = F(%4,.--, Eni Ur,--+, Urs 4), 

to (1), we arrive at the ordinary problem, with the additional restriction on the final value of the n + Ist 

coordinate: Xn44(T) =A. 


The problem about the minimum transition time of the system from one fixed state to another is also a 
particular case of the problem about the minimization of one coordinate. 


Indeed, by combining equation 


with (1), we may convince ourselves that minimization of time means minimization of the coordinate x,, , 4. 
In this case, the final values of the remaining coordinates are given, and the time T, evidently, should not be 
fixed, 


Hence, a broad class of optimization problems leads to the problem concerning optimization of one 
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coordinate of the system, It will be more convenient for us, however, to examine a somewhat more general 
problem (although one which easily leads to the problem about optimization of one coordinate), To wit, we 
will pose the problem of optimizing a linear function of the final values of all the coordinates of the system, 


n 

that is, the quantities S = > Ctx (T), where cy, are certain constants, The quantities S can be treated as 

ket 
the scalar product of the vector x (T) = [xq (T) sees; X,, (T)], which gives the end point of the trajectory x (t), and 
the given vector c =(C4,..., Ck). In other words, we require optimization of the projection of the vector x (T) 
on the given direction (cy,..., c)). By the same token, the requirement of maximum (or minimum) size of S$ 
means that we attempt to move the system as much “farther® in the direction of the vector c as possible (or as 
much *farther® in the direction of the vector — c as possible), 


As we have seen, various restrictions can be imposed on the final state of the system (that is, on the 
position of the final point ‘Of the trajectory x (T) in the phase space X), In general, these restrictions can be 
formulated as a requirement of a transfer of the system to some fixed set G of the phase space X, described by 
some combination of equalities and inequalities, In particular, if the final values of some coordinates of the 
system are precisely fixed, and the remainder free, then the set G is a linear manifold in the space X (*linear*, 
“plane,” etc). 


Hence, we will examine the following problem. It is necessary to select a control u (t) from the group of 
admissible controls which transfer the system (1) from the point x (Ty) = x’ to the fixed closed set G of the phase 


space, in such a way thatthe sum S = S cyai(T) at the given time t = T will assume a minimum (or 
i=1 


maximum) value, 


In order to be definite, we have formulated a problem in which the control time T-T, is fixed, However, 
we will see consequently that in those problems where the time is not given beforehand (in particular, in the 
problem about minimizing the transition time), completely analogous results obtain, 


n 
The control u (t) which makes the functional S = Dj cirs (7), minimum (maximum) will be called 
1 


min-optimum (max-optimum) in S, for brevity. 


Mathematically, the problem examined by us is a generalization of Mayer's classical problem in the 
calculus of variations, However, the presence of restrictions on the positions of the controlling elements, lead- 
ing to the requirement that the controls u (t) must belong to the closed set U, generally speaking does not permit 
one to use the methods of the classical calculus of variations, Therefore, the solution of the problem formulated 
above requires the development of new mathematical methods, 


It should further be noted that the theory of optimum systems is not to any degree exhausted by the solution 
of the problem formulated above. Firstly, various complications in the statement of the problem are possible. 
In particular, accounting for restrictions, which can be imposed on the phase coordinates x;,..., X,,, has great 
practical value, Secondly, in a number of cases problems arise which do not lead to optimization of the final 
value of the coordinates of the system and which require the solution of variational problems for functionals of 
another type (for example, the problem about the minimum of the maximum deviation of some coordinate from 
a given level), These and analogous questions are not discussed in this article, 


2. The Problem with the Free Right End of the Trajectory 





In this section we will examine a particular case of our basic problem, To wit, we will consider here 
that the right end of the trajectory x (T) is “free”, that is, no restrictions are imposed on the final values of the 
coordinates (the set G occupies all of the space X), The results, which pertain to the case being examined, 
prove to be completely analogous to the results obtained in the solution of the general problem, differing from 
the latter only in details which concern the boundary conditions (transversality conditions) and at the same time 
permit more simple proofs, Moreover, problems of this type are important in themselves, for example, in a 
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number of questions about the theory of the dynamic precision of control systems or in solving problems of control 
in mathematical economics, Therefore it is well to examine this simpler problem separately, 


Below, we will formulate and discuss results leading to the solution of the problem being examined. Full 
proofs of the theorems are in the appendix, 


We will bring into consideration a combination of n time functions py (t),... Py (t), which form a variable 
vector p (t) = (py (t)seees Py (t)), Which has a direction, at rtime t = T, opposite to the direction of the vector 
C =(Chyeeey Cn), that is 


p(T)=—c, =1,....%). (3) 


We assume here, for simplicity, that the modului of the vectors p(T) and c are equal, The variables 
pi (t) are subject to a set of differential equations of the following form 


n 


: é oes Day Ug, oe yl; 
p(t) = —)'p, fg (iy +++) Dys May «++ 9M it) re (4) 


" Oz, 





by the same token, we give a law which the vector p (t) can be changed in accordance with, Let us note that if 
any control u (t) is given, the vector p (t) is uniquely determined from Eqs, (4), where conditions (3) play the 
role of boundary conditions, Indeed, if the control u (t) is known, then, by solving the system of equations 


x = fi(z, u, t) inf, . ows n), 


for which the initial conditions [xj (T,) = x}. i=1,..., 0)] are given, we can find the functions x; (t). Putting 
these functions in (4), we would have a set of linear differential equations with variable coefficients and with 


boundary conditions (3), for the magnitudes of pj (t), The solution of these equations would permit us to deter- 
mine p (t) as a function of time, 


The vector p (t) will be called an *impulse* which maps points at time t. 


n 


At each time t, let us examine the scalar product >) Pi (t) zi (t) of the vector p (t) and the velocity 
1 


vector x (t) which maps the points. Bringing Eqs. (1) to our attention, we will express this quantity in the form 
of a function of 2n + r+ 1 variables Xqye00y Xni Pysevos Pri Ugeooes Ups t: 


H(z, p, u, i= >) Pofa (Say - 2-5 Bns By, +s 09 Be; &). (5) 
-s=1 


By using the function H (x, p, u, t), which plays a basic role in all further discussion, Eqs, (1) and (4) can 
be written in the form 


° oH . oH . 
7 Op, ’ ie se (i = 4;..<6,8) (6) 


i 


with the boundary conditions 


a(T,)=2, p(T)=—a (i=4,...,0). (7) 


Let us turn our attention to the fact that Eqs. (6) concur in form with the canonical equations of Hamilton 
in analytical mechanics. The function H is analogous to the Hamiltonian, and the vector p(t) to the impulse 
vector, which explains the name given by us to the vector p(t). The connection between the equations of the 
theory of optimum systems and the equations of mechanics is discussed also in appendix 4, 


Now let u (t) be some admissible control, and x" (t), p" (t) be the position and d impulse which maps the 
points at time t for the control u(t), Let us put the values of the variables x} (t), Py (t) in the function 
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H (x, p, u, t). Then we will obtain some function 
K (t, uy,...,Uy)= A (a“(t), p*(t), u, 2), 


which is a function of r variables u =(u,..., u,) for every fixed time t, that is, points which belong to the closed 
set U of the space R, We will say that the control u (t) satisfies the maximum (minimum) candition if the 
function K (t, u) reaches an absolute maximum (minimum) in the set U at any time t (Ty =t = T), for values 

of the variables equal to the values of the control at the same time, that is, at uj, = uj, (t)(k =1,..., r). It appears 
that fulfilling the maximum (minimum) condition for u (t) is a necessary condition for min-optimality (max- 
optimality) of the control u(t); this leads to the theorem expressing the Pontryagin maximum principle for our 
problem, * 


n 


Theorem 1, If the control u(t) is min-optimum (max-optimum) in § = citi (7), then it satisfies 
1 


the maximum (minimum) condition, 


Hence, the optimum control u (t) must be chosen at each moment so that the function H (x" (t), p" (t), 
u, t), which is a scalar product of the vector p" (t) and the velocity vector xu (t), is minimized (maximized), 
Therefore, the geometrical sense of the theorem is that the optimum control tries to “drive away" the image 
point optimally in some direction determined by the vector p" (t), at each moment, At the time t = T, when 
p" (T) =~ c, the direction of *maximum dispersal” concurs with the direction of the vector c if we are seeking 


n 
the maximum of the functional 5S = dois (7), and in the opposite direction to c if the minimum of the 
1 


functional is sought, This property of the optimum control becomes especially clear in the following example, 
where the equations of the object have the form: 


a = fi(uy,.-+,Ur) (i= 4,...,0), 


that is, the right halves of the equations do not depend on the coordinates, Forming the differential equations 
for the impulse p" (t), we find that the vector p" (t) is constant in our case, and is equal to 


p’(t)=—c 8 (T,<t<7). 
at any time, 


Indeed, Of,/8x, = 0, and it follows from Eqs, (4) [or (6)] that p" (t) = 0 and, in view of conditions (3) [or 
(7)}, we obtain the equality which is written, Therefore, the direction of “maximum dispersal” is always con- 
stant in the given case, and at any time, the control is chosen so that the maximum possible velocity is im- 
parted to the direction of the vector c for max-optimality, and the maximum possible velocity in the direction 
of the vector — c for min-optimality. Naturally, such a simple control law does not prevail in the general case, 
since if the right halves of the equations of the object depend on the coordinates, then interval “forces® will act 
on the image point, and these must be taken into account, The “internal” forces are accounted for by introducing 
a variable vector p" (t). 


In the future, we will omit the index u in the designations p" (t), x" (t). 
The maximum principle permits one to form a set of differential equations whose solution determines the 


*We note that fulfillment of the maximum condition is essential if the functional S=)c,x,(T) is to be 

1 
minimum (and vice-versa), This is connected with the fact that the final value of the vector p (t) was defined 
as a vector opposite to the vector Cc =(Cy,s0.6, Cp). One can postulate that p(T) = +c, but in this case the for- 
mulation of the theorems would be the reverse of analogous formulations assumed in [6-12], which would lead 
to confusion, 
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optimum trajectory and optimum direction, Indeed, the quantities ly seee, Uy are determined from the maximum 
(minimum) condition, at any time, as functions of the variables Xpreces X35 Ptovore Pras £3 


Ux = by (Z1,-++52n; Pry,---,Pns t) (e=1,...,7F). 


Putting the expressions obtained in the right halves of Eqs, (6), we will get a set of 2n differential equations 
with 2n unknown Xq (t)seoos Xn (t); Py (t)reees Pr (t), Whose solution, under the boundary conditions (7), will give 
the optimum trajectory x (t) and impulse p (t). Knowing x (t) and p(t), we obtain the optimum control u (t). 


Let us take the simplest example, Let it be required to find the control u(t), which makes the integral . 


- 


1 
=\@ + u*)dt, maximum, if the equation of the object has the form: 
0 


t= —axr+u, z(0) = 2°, 


and no restrictions are imposed on the control u (t) (that is, the set U occupies all of the space R— in the given 
case, the whole real axis -@ < u< +), 


t 
Introducing the variables 2, (t) = -2/(t), z, (= \ (x? + u*)dt, we get the system 


0 


. . 4 1 
= —antu, m=tat+ tut 


The functional which must be minimized is S = xg (T), that is, cy = 0, cg = 1, According to equality (5), 
we write the function 


H(z, p, u) = —ap,x, + = Pat? + pu + = py? 
and the equations for the impulse (second group of Eqs. (6)) 
Pr = ap; — P2%, Ps = 0 
with the kind of boundary conditions of equalities (7) 
z,(0)= 2°, 2,(0)=0; p(T)=0; p,(T)=—1. 
We notice at once that p, (t) is constant and p, =— 1 always, so that the function H has the form: 
H = —ap,x, — 42 + pu — +0". 

According to the maximum principle, the function H must be maximum in u for any values of x and p, In 

the given case, the maximum of H is reached at the value u = p, (t). Putting the value of u obtained in the 


equation of the object, we obtain the following equations for the determination of x, (t) and p, (t) (the function 
X_ (t) does not interest us): 


t= — ar, + pr, Pr =2, + ap, 
with boundary conditions xq (0) = x’, Pi (T)= 0, The equations are easily integrated; 


Z(t) = Cye* + Cre", =p, (t) = Dyc™ 4+- Dye, 


where ) = /a? +1 is the root of the characteristic equation, and the constants Cy, Cg, Dy, Dz depend on x*, Since 
u (t) = py (t) the optimum control is also determined, 
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The problem which was examined can evidently be easily solved without the use of the maximum principle, 
since there are no restrictions on the controls. The given example was introduced only to illustrate the procedure 
of the solution coming from the maximum principle, 


In a very broad class of practically important problems it appears that the maximum of the function 
K (t, u) = H [x (t), p(t), u, t) is reached at the boundary of the set U, and the optimum control takes its boundary 
values. In particular, this circumstance always occurs if the functions wy (t),..., u, (t) enter the system of 
equations of the problem linearly, For example, let the equations of the problem have the form: 


z, = f1(%, 2) + U4, em = fo(%1, T2) + Us, 


and the functions u, (t) and u, (t) be limited in modulus so that; | uy] = Mg, |ug| = Mg (the set U is a rectangle), 
Then, forming the function 


H(z, p, u) = pif (21, 23) + Pofa (Zr, Ze) + Pils + Pate, 


we find that the maximum of H is reached at uy = My sign py, Ug = Mg sign pz. Indeed, H is maximum in uy, and 
Ug when uy and py, Ug and p, respectively have identical signs (so that the products p,uy and pylly are positive); 

in this case, uy and ug reach a maximum in modulus,* Therefore, the functions uy (t) and ug (t) can only *switch" 
from one boundary position to another, no matter how py, (t) and py (t ) change with time; p, (t) and p, (t) are 
determined only at the moment of switching. 


Theorem 1, formulated above, which expresses the Pontryagin maximum principle, gives only the necessary 
conditions of optimality, It is easy to show examples in which the control satisfies the maximum (minimum) 
condition, but never the less is not optimum, From the viewpoint of the remaining conditions, it appears that 
the control, which is subject to the maximum (minimum) conditions is optimum, generally speaking, only in a 
definite local sense, if some additional conditions are not fulfilled, which are connected with the *sharpness® of 
the maximum (minimum) of the function H [21], However, for a system whose controls are linear relative to 
the coordinates 


Li = >) Gis (t) te + 9: (uy,---, Ur), i=—1,...,m, (8) 
1 


the following theorem obtains, 


Theorem 2, The necessary and sufficient condition of min-optimality (max~optimality) of the control u (t) 
in system (8) is the fulfillment of the maximum (minimum) condition for u (t). 


It must be noted, that if the maximum (minimum) condition fs not sufficient, generally speaking, even so, 
it permits one to determine uniquely the optimum control in the majority of practically important cases, Indeed, 
if the optimum control exists (which is often evident from physical considerations), and the control satisfying the 
maximum (minimum) condition is unique, then this unique control is optimum. 


APPENDIX 
1. Increment of the functional during a change of the control, Let x; and p; satisfy the equations 
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* We designate the function 


s={ 1 for y>9, 
~ (—1 for y<0, 


by the symbol z = sign y, which is not defined when y = 0, If at any time p, (t) = 0 and p, (t) = 0, then H does 
not, in fact, depend on wu, and u,, and the maximum principle does not permit one to find the value of the 
optimum control (see appendix 3.3), 

















n 
H (x, p, u, N= Pil (z, u, t). (10) 
1 


The functions f; (1 = 1,..., 1) will be considered continuous in all of the arguments (x, u, t), and to have 
continuous partial derivatives in the arguments (x, u), up to and including the second derivative, 


The equality 


T 


n 
I iz, p, wu = \ [ Sae—Hee. p,u, »| dt = 0, (11) 
ae 


follows from (9) and (10), and {s justified if the first group of Eqs. (9) is satisfied under arbitrary controls and 
boundary conditions. We will choose some control u(t) @ U and examine its increment u(t). Let x (t), and 
p (t) be a solution of (9) for the control u(t) and some boundary conditions, and x (t) + 6x (t), p(t) + dp (t) be 
a solution of (9) for the control u(t) + u(t) with the same boundary conditions. Let us examine the difference 
4 = I1(x + 5x, p + dp, u + u)— 1 (x, p, u), in which A= 0 always. It follows from (11) that 


T n T n 
a=\ ) (pba, + Bp, z,) dt + ( ) Bp,ba,dt — 
Tel f, 1 
3 [H (x + tz, p + Bp, u + bu, t) — H(z, p, u, t)] dt. (12) 


We note that the functions x (t), p(t), 45x (t), Sp (t) are continuous, and their derivatives exist and are 
piecewise continuous, ; 
We will integrate by parts: 
T T T T 
‘ T ; ; T , 
\ p,dzdt = p,bz; [ + p,bz,dt; \ bp,dz,dt = bp,tz, i _ \ bz, p,dt. 


Hence, keeping (9) in mind, and introducing the vector y = (Ygse00s Yan) (Xi = Vie Pi = Yn 4ie 1 = Loess 2) 
for brevity, we obtain 


ee rs T .m aoe u, ). 
21 1 | 
Ta n T on 
. i OH (y, u, t 
Sl spat = i sp,tz, |. » +e re oy by, dt. (14) 
| 1 T, 1 





Using Taylor's formula, we put the last integral in (12) into the form 


T 
Jom u + tu, t)— HAH (y, u, t)] dt 1=\ (Gy, w+ tw, t)—H (y, u, t)] dt + 
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T, 1 


yy st t) 








¢ 2 cay d,by, u + Bu, ¢) 
by, +z) 3 ys (y + dy, u + bu, (15) 
8.q= 


Oy, PY, a 











where 0< &(t)< 1, Since 


1 0H | 6H (y + by, u+ bu, t) 1 0H (y, u, t) 
2 yt ae es > oy, We 





by putting (13)-(15) in (12), we get 
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T 
A=>\(p,+ >t, )ee, |. — \u (y, u + bu, t)— A (y, u, t)] dt — 
1 T 








T on 
i 0H (y,u+tu,t) dH (y, u, t) | 
ra \ xl dy, <2 1’. 2 
1 in 0H (y + Ody, u + Bu, t) 
1 , ’ 
-> \ x ayy, by bygdt + 
08, Q=1 


7 
es 4 (OH iy + by, u + bu, t) 0H (y, u + bu, t) 
5 das 


Oy, dy, | by, dt. 





-b43 
ee) 


Using Taylor’s formula once more for the integrand of the last integral, and keeping in mind that A = 0, 
we will have 


T 
1 T 
Yi (r+ TH) 8 4 = Sur, u + bu, t) —H (y, u, t)] dt + N (16) 
1 . 


where 9 =, +, and 





T on 
i 6H (y, u + bu, t) 0H (y, u, t) 
n= 5 a> iy, a | syyae. 


T on 
aot 1 [ PH (y + Ody, u+ du, t) OH (y + by, w+ bu, #) 
m= \ > oy “i dy,t, | by, by,dt, (17) 





6,q=1 


in which 0 < 8 < 1, 0 < 8, < 1. Fixing the boundary conditions in the form 


' 2, (T.) = 2%, p,(T)=—<¢, (i= 1,...,8), (18) 


in view of the fact that 5x; (Ty) = 5p; (T) = 0, we get, finally, 


a T 
BS = Die,t2,(T) = — \ (Hp, u + bu, 1) —H(e, pu, t)] dt—, (19) 
i T, 


2, Evaluation of the remaining term in the formula for the increment of the functional, First, we will 
get an evaluation of 5x, (t) and 5p, (t). 


It follows from the first group of Eqs, (9) that 





tz,= f; (2 + Ba, u + du, t)h—f,(z,u,t) (é=1,...,R). 


Since the initial conditions are fixed and 5x; (T») = 0, then 
t 


b2,(t) = \ {4 (2 + ta, u+ bu, t)—f, (2, u, t)] dt. 
tT. 














By using the Lipschitz conditions, which always hold for our assumptions, we obtain the inequalities 


t n t r 
18, 1< | KY 06,8) [de + | x Ym conta (i= 1,...,n), 
1 


° ~ 1 


where K is the maximum of the Lipschitz constants for the functions fj. 
Summing and strengthening the inequalities, we have 


bs t 
x0) <Kn\ >} | Bury (2) [ds + (nx (ode, 


‘ 
* 


n 
where X(t) = }}|%2,(t)|, It is easy to show that the inequality 
1 : 
X(t) [at (ox (2) dz 
T, 
follows from the inequality 
X(HN< act—T.) 


(a and b are positive constants), 


T 
Therefore, putting a=Kn \ |3u,(s)jaz,o=—Kn and substituting t for T, we get 


~Ms 


Tr 
X () <Knekmt—10 { >}! bu, (2)|ds (To t<T). 
T, 1 


Thence, we have analogous inequalities for 5x, (t): 


T + 


82, (0) < Mi \ Si bu, (2) | ds (i=i,... 


T, 1 


> 
= 
2° 

a\ 
4 
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where My = KneX%™(T~T9) is a constant which does not depend on §u (t). 
Examining the second group of Eqs, (9), we obtain for p, (t) 





BP; _ = 


dH (x + tz,u+ bu, p+ p,t) OH (z, p, u, t) 
i Ox; ae Oz, ° 


in which 6p;(T)=0, Therefore, by now using the evaluation for 5x4 (t), we get, in a completely analogous 
fashion, 


T r 
bp, (t) < Ma | ‘| Bu, (2) | de (i= 1,....", Te Qt 7). 
T, 1 


Hence, the evaluation 


Tr 


by, (I< m;\ dj | Buy (2) Ide (a= 4,...,2n, Te Qt ST), (20) 
T.1 











occurs, where M is a constant which does not depend on u(t). Applying (12) to the remaining members n,; 
and mz in(9), using the Lipschitz conditions for the functions 0H/ dy, and having the boundedness of the 
functions 8°H/@ys®yq in mind, we obtain 


. 
ini <4(\ >) | Bu, (2) | az ) (21) 
| 


for n = ny + M2, without difficulty, 


If the function 5u (t) is different from zero only in the interval [t,, tz] of length t,— t, = r, then, using 
the Cauchy-Bunyakovskii (Schwartz) inequality, formula (21) can be re-written in the form 


t, tear 


of 2 
lnl< Ax\( \ jou, (1) di or |nl< Br \ >} buf, (4) at (B =rA). (22) 
t, 1 i; 1 


3. Proof of theorem 1, To be definite, let us examine the case when the minimum of the functional S 
is sought. For any change 6u (t) of the control u (t), the following must hold; 





bS > 0. (0) 


Let us assume that the maximum condition does not hold at t = t*, that is, a vector u* =(Uf,..., u,) 6 U 
exists such that 


H [z(t*), p(t"), wu’, O]>A[x(e’), p(t’), u(t), &). (2) 


In view of the piecewise continuity of u(t), and the continuity of the functions x (t), p (t), and H (x, p, 
u, t) some interval [t*, t"], which includes t*, is always found, on which the functions F,; (t) = H [x (t), p(t), 
u*, tj) and F,(t) =H (x(t), p(t), u(t), t) are continuous, and consequently, are uniformly continuous (if t = t* 
is a discontinuity point of the control u (t), then, by agreeing that u(t*) = u(t* — 0), one can take t* = t*), 
Therefore, the existance of a number a> 0 follows from inequality (°), such that the inequality 


H [x(t), p(t), u®, t]—A[x(t), p(t), u(t), t]>a (22) 


holds for all ¢€[t;, t2], where [ty, t,] are sufficiently small intervals which belong to the interval of continuity 
[t*, t"]. Let us now examine the varied control u (t), which is equal to the control u (t) everywhere outside the 
interval (ty, t,], and equal to u* at [ty, tg). The control u(t)@U is admissible, Then, using (19), (::), and 
(22b), we get 


te 
85 == — \ [H (y(t), u*, t)— A (y(t), u(t), HJ) dt—n< 


ty t, 


<—|fad—|n1] <-S[s-eymto]a ; 


r 
Since }'buj(t) is bounded, it is always possible to choose a sufficiently small r = t,— ty, for which 
I 


the integrand is positive and 5S < 0, which satisfies condition (0), The case of the maximum of the functional 
§ is examined analogously, The theorem is proved, 


4. Proof of theorem 2, The necessity of the maximum (minimum) condition has been established by 
theorem 1, We will prove the sufficiency of this condition in the case of linear systems, The Hamiltonian of 
the system has the form: 





n n 


H= > Ay (t) Py 2, +- >) Pi % (u), 
i, Re] 1 














and the impulse satisfies the relations 


Py = — D} Opi (8) Pp (i=41,...,n). 
1 


Let us examine the remaining terms (17) in formula (19) and prove that 7 = 0. Indeed, a H/ 0x; 0 x; = 
=0, 2 H /dpi op; = 0, and the quantities ay; = 3 H /0p; 9x; clearly do not depend on x and p. Therefore 
Nn; = 0. In the expression for ny, the differences OH(y,u+ du, t)/ dx, OH (y, u, t)/ 8x; disappear identically, 


0H he 
since oe, =— dee (t) p, clearly does not depend on u, On the other hand, the change of the impulse 


5p (t) satisfies the relations 


n 


py = — DY} ay (t) Bp, (i= 1,--- 9m), 
1 


which do not depend on 6x and 6u, But since the boundary conditions for the impulse are fixed and 5p (T) = 0, 
then 5p (t) is identically equal to zero for any t. Therefore, ny = 0, Hence, the formula for the increment of 
the functional takes the form: 


ss = —| [H(y, u+ bu, t)—Aly, u, #)] de. 
T 


If the maximum (minimum) condition holds, then the integrand is not positive (not negative), and 
6S 20(6S 0), Q.E.D. 
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THE PROBLEM OF OPTIMIZING SYSTEMS WHICH 
CONTAIN ESSENTIALLY NONLINEAR ELEMENTS 


E. P. Merkulova 


(Moscow) 


A method is given for finding the optimal weight function of the linear portion 
of a stationary system which contains an arbitrary number of essentially nonlinear 
elements, The problem is solved for various methods of connecting the linear portions 
and the nonlinear elements. The input is a stationary random function of time con- 
sisting of signal plus noise, The linearization of the nonlinear elements is carried out 
by a statistical method, 


The increased accuracy required of automatic control processes in the design of optimal systems makes it 
ever more compelling to turn our attention to the problem of taking nonlinear characteristics into account. 


The literature contains a number of works [1-3] which present various methods for designing optimal non- 
linear systems. However, either the final results are given in a form which is difficult to use in practical cal- 
culations (1), or the nonlinearity considered is a continuous function [3] or must be given as a power series [2], 


In the present work, we present a method of designing optimal dynamic systems which contain essentially 
nonlinear elements which are either relay-type elements with dead zones or elements with limited zones of 
























































linearity. 
Mt) . Wt} Wt) The design of the nonlinear elements is carried out by 
4 J Kazakov's method [4], Le,, a nonlinear element is replaced by 
a linear one in such wise that the mathematical expectation 
WT, Lae, and the dispersion of the element's output signal remain as they 
were previously. 
Fig, 1. 1) Nonlinear clement No. 1; 5) 1. Posing of the Problem and the Problem'‘s 





nonlinear element No, 2; 3) linear 
portion No, 1 and 4 is linear portion No. 2. Solution 


We consider a system consisting of any number of linear 
portions and nonlinear elements which are so connected in a forward circuit and in a feedback path that the 
system may be represented in the form shown in Fig. 1. 


We shall denote by @(p) and ¥(p) the transfer functions of linear portions No, 1 and No, 2 respectively. 
We assume the following with respect to input signal X (t): 
1) X (t) = m (t) +n (t); 


2) m (t) = mg + m' (t) is a stationary random function of time with mathematical expectation mp», correlation 
function R,, ( r) and spectral density S,, (w); 


3) n (t) is a stationary fandom function of time with zero mathematical expectation, correlation function 
Ry (7) and spectral density S, (w); 


4) m (t) and n (t) are noncorrelated; 

















5) the input signal*s spectral density, Sg ( w) = Sp (w) + Sp (w), is a rational- fractional function of w. 





We shall denote by w (t) the weight function of the system as a whole, and we shall assume that 





w(t) =0 (0O>t>T), 


where [0, T] is the time interval during which we locate information about the input stimuli, i.e., the system's 
memory. 


After the application of statistical linearization, nonlinear element No, 1 is replaced by coefficients 1, 
and 1;, and nonlinear element No, 2 by coefficients ky and ky (by constant and random signals, respectively). 
These coefficients are functions of the mathematical expectation and the dispersion of the signals applied to the 
inputs of the corresponding nonlinear elements, 


As a result of these replacements, the system as a whole will have two different transfer functions: one, 
%; (p), for the random signals and one, , (p), for the constant signals, for which we can write: 








Io (p) 11 (p) 
= = 1 Cee 1 
PoP) = Trike FR)’ = TFTLOM TO - 
Since the types of the nonlinear elements are given, the problem reduces to the finding of the optimal 
transfer function for the linear correcting element, considering the transfer function of the system's invariant 
portion as being given. We shall assume that ¥ (p) is given, Then,by eliminating @(p) from the two equations 
in (1), we obtain 


Op) 
Do (P) = TF () Oi)’ (2) 





where 


q=7 1 v= ho—kyg. 
By assuming that q is close to unity and y is close to zero, we can present (2) in the form of a series 


Dp (Pp) = = {D, (p) — ¥¥ (p) D5 (p) + VE" (p) OF (p) + - 
vee + (— 1)" vB (p) OF (p) + - - -}- (3) 


If we go from transfer functions 4» (p) and @, (p) to the corresponding weight functions, we (t) and wy, (t), 
we get 


t 
wg (t) = — {ws (t) — v\ gs (¢ — 9) w, (0) 0 +... 


(HA gua (t— 8) wy (0) A9-4...1, (4) 
where 
‘. t,—t, ty—ty—..-—tn_4 ts 
Wy (t;) = \ W, (ta) dt, \ W, (ts) dts... \ W,(tn) Wy(t;—te —o.-. —tn) dtn = \ Wy, (ta) Ln—; (ty — tg) dt, 
0 0 0 9 


ty tt, t—ts—---—tn_y 


gn (ts) =| g(t) dts \ g (ts)dt,... \ 8 (tn) g (ty —te— ..- —ty) dtn. 
0 


It is obvious that 


L, = 1, &o (ts) = 1, 
L, (%) = w, (*), 8: (41) = g (é), 


1 


t 
L(t) = \~, (t2) w, (t, — t,) dt, 82 (ty) = \s (ta) g (t1 — ty) dts, 


Ty (t) = \ U', (t,) dt. \ W, (Ts) Wy (t, — tT, — T,) dt, 
0 0 


t t,—t, 


83 (t;) = \e (t,) dt, \ & (ts) g(t; —t, —ts) dt, ,etc. 


0 0 
where g (t) is the weight function of linear portion No, 2, 
Thus, the system's output, Y (t), will have the following form: 


T T 
Y(t) = Mme \ wo (2) de + \ [m’ (¢ — x) +n(t —2)] w, (2) de. 
0 0 
Then, the error, € (t), in signal reproduction will consist of two parts: 


a nonrandom part 


4 
> = H (p) mg — mg \ wo (2) de; 


and a random part 


T 
€1 (t) = H (p) m’ (t) —\ [m' (¢—2) + 0 ¢—2)] 4 (0) dr, 


0 
where H (p) is the transformation operator. 


If we introduce the static coefficient, Cy», of dynamic error and use approximately the same reasoning as 
was followed in [5], we obtain the following expression for Cy by making use of formula (5); 


T 
Cy) as Hy—\wo (t) dr, 


where Hy = H (p) | pmo 


If we substitute the expression for we (t) from (4) in (6), we obtain a condition relating to the weight 
function wy (t): 


C= H,—— ie (t) tx—fdel ee m+ see 
T 


eee (—1)" 7 v" (atl go (t — 8) wy, (8) dd +-.. +} : 
















We shall characterize the random component of the signal reproduction error by its mean square value 





my ‘ ) T > 
eo = lim \ {H (p) m’ (1) —{ [m’ (¢—t) + n(t—*)] wy, (x) de} dt. (8) 








Thus, the problem of optimizing the system reduces to the finding of a weight function wy, (t) which 
minimizes expression (8) and which satisfies additional condition (7). 





This typical problem of variational calculus leads to the following integral equation 


T T—t 


| n(t—) + Ralt— my (d= 2 ft—a | HOA —t—Hde4 
0 — 0 
+3v { Ly (t) fo(T —t—t)dt+... 
” mt 
ot (—1)"" Vn \ Int (*) faa (T —t —1) de +...) + 
oo \ Rm (t — t) x (t) de, (9) 


where y, isa Lagrange multiplier, 


fn(u) = \ en (9) d8, x(t)= — \ H (ja) ei do. 


0 


Equation (9) can be solved by the method of successive approximations, For this, wy (t) must be given in 
the form of a power series in the small parameter v: 


W, (t) = Vo (t) + vv, (t) +- vv, (t) +... (10) 
By substituting (10) in Eqs, (9), we obtain, for each approximation vq (t), an integral equation of the type 


T 


\ [Rin (t — 2) + Rn (t —2)] %: () dt = F (2), 


methods for the solution of which are presented in [6], 


By thus determining weight function w, (t), and then @, (p), we find, from formulae (1), the sought- for 
transfer function,  (p), of the linear portion, 


2. Practical Design Method 





Since q and vy are functions of the mathematical expectation and the dispersion of the signals applied to 
the inputs of the nonlinear elements, the problem can only be solved by the method of successive approximations, 


The following method of solution is suggested, 


1, Starting from the conditions of the concrete physical problem, one writes Eq. (9) and finds its solution, 
limiting oneself to the necessary number of approximations, The expression obtained for w, (t) will contain the 
unknown constants A; which are found by substituting wy (t) in Eq. (9) and in condition (7), Since the unknown 
quantities q and y also enter in (7), the A, will be functions of these quantities. 


2. Starting then with some initial values of q and vy (knowing the type of nonlinear elements, the roughly 
determined statistical characteristics of the signals applied to the nonlinear element inputs, etc), one computes 
the constants Ay and, consequently, wy (t) and 4 (p). 


1306 





3. By finding concrete expressions for the functions 


1, (m ’ Dy) 
q = hime Dey 24 v¥,= ky (my, Dy) — k, (my, Dy) q (my, Dy) 





(from [7], for example) and by taking (3) into account, one then finds the values of my, Dy, mz and Dz from 
the following equations: 


my = mo {1 — ~{®, (0) — v¥ (0) 2 (0) + v™¥* 0) @ (0) — lf, 
Dy = 2 \ S4(w)|A—K¥ (jo) ®; (Fo) de, au 
0 
mz =~ {¥ (0) D, (0) — v¥* (0) Dz (0) + v°¥? (0) Ds (0) —.. .}, 


@ 


Dz = 2\ S_(w)|¥ (ja) ®, (jo) [Ado 


0 


These equations can be solved graphically, 


4, From the values found for my, Dy, mz and Dz one computes new values of q and y and, consequently, 
new values of the constants Aj, Otc. 


These computations continue until the latest values of the A, differ only slightly from the previous set of 
values, 


5. One then finds the final value of wy (t), then ®, (p) and, from formulae (1), @ (p), the sought-for 
transfer function of the system's linear portion. 


3. The Case when the SignalgContains a Regular Function of Time 





In the case when the signal contains a regular component, the method described above cannot be given a 
rigorous basis since the coefficients of the statistical linearization will be functions of time and, consequently, 
formulae (1) will not, in general, be valid, However, in many practical cases, particularly when the regular 
component is a slowly varying function of time, one can cite approximate methods of finding the optimal weight 
function, these methods, however, making no pretense of rigor. 


Let X (t) =a (t) + m(t) +n (t) be the input stimulus, where a (t) + m (t) is the signal component, We 
assume that a (t) has r first bounded derivatives, The assumptions with respect to m (t) and nm (t) are as before, 


In order to use the method previously presented, it is necessary that, in formulae (1), ky (t), ky (t), 29 (t) 
and 1, (t) be replaced by some average values of them, kg, ky, 2g and ly, Then, the only quantities subject to 
change are condition (7), Eq. (9) and the expressions for the mathematical expectations in formulae (11) which, 
for the given case, take the form: 


T T t 
C, =U! A,+ (— yt — ff tw, () de — vi t! de| gx (t — 8) wy (0) d0 + amie 
E.. 0 0 
P +A ede Leaner a+. | (i=0, $00 9?) (12) 
0 


T—t 


T r 
{ Rn (t= +) + Ra (t— 9) 0, (2) de = + You ft!—av) LAT, t+ dat 
é l=0 











Z T—t 
gu oe +34 | LAT, t+)dr—... 


| 0 


.+(—1)"*nv"" 


1 





t 


Ln—1(t) fny(T, t +t) de +.. ; + 








Rin (t — t) x (2) de, (13) 


beg . 


Fig. 2. 7 


T—t 
where /,(7, 1) = ( (t+ 6)' gn (6) d6, ~, are Lagrange multipliers, 





rg 0 r go o 
l=0 l=0 


where , (p) is taken from formula (3). 


The expressions for Dy and Dz in (11) are not changed, since with the given choice of limits of inte- 
gration from 0 to T, the correlation functions of functions V (t) and Z (t) depend on the difference of their argu- 
ments and, consequently, their dispersions remain constant, 


The method of solving Eq. (13) is the same as for Eq. (9) since they do not differ in principle from each 
other, 


Since, in (14), the mathematical expectations depend on time, it is suggested that the interval [0, T] be 
decomposed into a series of discrete values of t, tj, and that graphical solutions of Eqs, (14) be found for each 
of these, The resulting functions, q (t) and y (t) are then averaged, and these average values are substituted in 
the equations for determining the cofficients A, in order to then render them more exact. 


All remaining computations are carried out analogously to those for the case of a constant signal. 


Example 1, We consider the systern which consists of a linear correcting element and one relay element 
connected in the feedback path, the characteristics of the relay element being shown in Fig. 2. 


Let H(p)=1,m(t)=mo, Co=0, R,(t) = ee". 
Then, 


1 2 
5, (0) = 5, (0) =— ao y=, v = ko — ky. 


If we limit ourselves to the zero*th and first powers of y, we obtain Eq, (9) in the form 


T T—t 
\ R,,(t — t) w; (t) dt = vy, [1 — 2v \ w, (t) dt]. 
0 0 


We substitute the first two terms of expansion (10) in Eq. (9). Then vo (t) and vy (t) must satisfy the follow- 
ing equations: 


T 
\ R,, (t — *) 0 (t) dt = Yo. (15) 
0 
T T-t 
\ R,, (t — 1) v1 (*) dt = — 2%0 \ Vy (t) dt. (16) 
0 0 


1308 








If we recall the expression for solving the equation of Zadeh and Ragazzini [8], we will write the ex- 
pression for v» (t) in the form 


Vo (t) = Ag + A,B (t) 4+ Ag’ (t —T). (17) 


By substituting (17) in Eq. (15) and requiring that it be satisfied identically, we obtain the equations for the 
coefficients and the expression for y9: 


A, =—, As» a 1s to= Ao: (18) 


Then, vo (t) will have the form: 


1 4 i 
v9 (t) = Aol 14+ 33 () + > B0—2)]. (19) 
By taking (18) and (19) into account, we can write the right member of Eq, (16) in the following way: 


T—t 
Act 1 he? 
— 2%, \ 0 (1) dz =— (T+) Ay + — Apt = Bo+ Bt. 
0 


Consequently, we obtain v, (t) in an analogous fashion: 
v, (t) = Ag + Agt + Asb (t) + Agd (t — T). (20) 


By substituting (20) in Eq. (16) and requiring that it be satisfied identically, we obtain the following 
algebraic equations for the determination of the coefficients A,: 





As AqT A, As Ay 
ate tert % Oo ane. 
2c? 4c? 28 4c? ‘ sx 
= A, = B, = Ab —q As = Bo =— (T+ =) AB 
If we express all the coefficients from (21) and (18) in terms of Ag, and take into account that 
w, (t) = v9 (t) + vo, (2), 
we obtain the following expression for wy, (t): 
1 i 2vApo / 2 * i 
w(t) = Ae {1 —2v4g(74+2)| +2404 [= —AB(7 4 =) 4+ ose—}. (22) 


To determine the coefficient Ag, we use condition (7) for the weight function wy, (t), this condition taking 
the following form in our example; 


T i t 
i= ( w, (t)dt — v \a \™ (+ — 0) w, (8) dO. (23) 
0 0 


0 


By substituting (22) in (23), we obtain a fourth-order algebraic equation for Ag 


QsAp + QsAh + Q2A5 + Qido + Qo = 0, (24) 
where , 
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Since we agreed to retain only terms in the zero’th and first powers of vy, we set Q, = Qs, = 0 in (24). Then, 
the expression for Ay can be written as 








—Q:4+V Q— 4020 
a 2Q2 





For further computations, it is necessary to give numerical values to the Qj, Let T = 10 seconds, a = 1 sec™! 
y = 0,05, Then, 


’ 
Ag = 0,0907 sec -}, 
w, (t) = 0.0816+4-0,0008¢ +- 0.08078 (t) + 0.0907 8 (t —T). 
For purposes of comparison, we give the value of Ag for the case when py = 0: 
Apo |,a-9 = 9.0833 S&C -1, 


In order to verify that the chosen value of v is correct, we compute the mathematical expectation my and 
the dispersion Dy at the input of the nonlinear element for the weight function wy, (t) obtained. 


For our example, 





my, = me [M; (0) —v(m,, D,) ©? (0)], (25) 
co 
D, = mae \ st | @; (jo) |? do. (26) 


0 


For simplicity of the calculations, we choose c” in (26) in such fashion than Dy *™ 1, Then, by setting 
1 (0) » 1,029 and 7 (0) = 1,058 in (25), and by assuming that m, = 1, we get 


m,, = 1.029 — v(m,) 1.058. (27) 


For the example chosen by us, the dependence of v on my can be expressed in the following way: 


2 — 1 a 4 2 
v= kp — ky, ho = 5 Olm,y)» ky = 7 V1 — 408 (m,) + Fe q tim, 





u 
© (u) = —_le 2 dt. 
V 2x 
0 
If we solve Eq, (27) graphically, we can convince ourselves that the chosen value, y = 0,05, is completely 
applicable and that further sharpening of this value is senseless since it changes only very slightly the value of 
coefficient A, and, consequently, weight function w, (t) and transfer function 4; (p). 


Thus, the optimal transfer function for the linear portion, @(p), will equal 


®; (p) 
® (P) = 7=0, 630, ’ 





since the coefficient k; = 0.63 for the value y = 0,05 chosen by us. 
Example 2, We solve the previous example when the imput signal contains a regular component, a (t) = 
re mg + t. 


. 


By again limiting ourselves to terms of the zero'th and first powers of v in the expansion of wy (t), we ob- 
tain the equations for vg (t) and vy (t): 





T 
| Balt) r(x) de = yo + wut 
0 





(28) 
T T—1 Tt 
\ R,, (t — t) v; (t) dt = —2 {re \ Vo (t) dt + 71 \ (t + t) a (7) ax}, 
1 0 0 0 
the solutions of which have the forms: 
Vo (t) = Ao + Ajt + Agd (t) 4+ Agd (t —T), 
V; (t) = Ag + Ast + Agt* + Apt® + Agd (t) + Apd (t —T). (29) 


By substituting (29) in Eqs, (28) and requiring that they be satisfied identically, and also taking conditions 
(12) into account which, in the given example, take the following forms: 


¥ t 
J mi (2) dv de | w, (+ —0) (0) 20 = 4, 
0 0 


wae { ws (t — 6) Ww, (8) d6 = C;, 
0 


tw, (t)dt—v 


ee | owns 


county 


we obtain the algebraic equations for the coefficients Ag 


Ao A, 
A= a ae 
Ae. the ! 
A= S4+F(7T +3). 
1 2 
A= —2(T +=) AR — 2(1*— 5) AeA — T*A?, 


As =2.A5 + a(r ~ =) AoA, + (3s + r) A’, 





A? 
As = 0, A= — 
2 2/2 z 1 (27% 6 
ee =) @+e(a-G—T) aa — 2(F 4745 =) AB, 
2 


‘aie + AoA;B; + A?By + ABs + AB, — 1 =0, 
ASEo + AgAyE; + ATE + AoE: + A\E,—C, = 0, 


where 
372 67 . 4 T i 
bo=—+(F+ 5+) B= —vr (> +4) 
47° 47 2 4 4T 2 
Ba=—(s- + -a—3) B-—w (51+ i — aw): 


7T* 27* a 32 
Ey=—vT Go - = a ~a@ as) 


=. 8=7($+2). ,=1(5+>). &=7(3+—+2)- 


By solving these equations for T = 10 seconds, a = 1 sec, C; = 0.1 second and y = 0,05, we obtain the 
following values for the A;: 


Ao = 0.300 sec=, A, = — 0.040 sec -2, 
A,=0.340, . Ay=—0.140, 
Ay = — 0.670 sec-!, Ay = 0.130 sec +t, 





































A,=0, A;= — 0.0005 sec-*, 
40 = 
ont A, = — 0.803, Ay = 0.035. 

t,°5 

6 

Thus, the weight function sought will equal 
wy, (t) = v9 (t) + vr, (t) = 0.266 — 0.034t a. 
0 7 My + 0.3008 (t) — 0.1388 (t — T). 
Fig. 3. By finding 


®, (0) = AT -+- As + Ap + As = 1,12, ©? (0) = 1,26, 
@; (0) = — (A, +. 4+ Ay 7. +AT)=—06, [@ (C0) =0,36, 
we can write Eq. (14), for mg = 1 and m, =~ 0,02, in the form 
m,, (t) =n (t), n(t) = (1,132 — v-1,253) — (0,022 — v-0,025) ¢. (30) 


The dispersion (exactly as in Example 1) we make equal to unity by properly choosing c*, 


By choosing three values of t, t, = 0, 5 and 10 seconds, we find the corresponding solutions of Eq. (30) 
graphically (Fig. 3). For each point of intersection, we obtain py (t,): 


v (0) = 0.03, v (5) = 0,06, v (10) = 0.10. 


By averaging, we find the new value, y = 0,06, which must be substituted in conditions (12) to determine 
the new values of the coefficients Aj. As a result of this computation, we find that 


Ao = 0.300 sec -', A; = — 0,039 sec~, Ay ~ 0.339, As = —0.129, 
A, = — 0.675 sec~!, As = 0.127 sec ~?, A, = 0. A; = — 0.0005 sec -*, 
Ag = 0.806, As => 0.033. 





As is obvious, the values of the coefficients A; changed very little, so that it is not necessary to sharpen 
them any further, 


Thus, the final value of the weight function wy, (t) is found to be 
w, (t) = 0.260 — 0.032t + 0.2918 (t) — 0.1278 (t — T). 


From this we find 4, (p), after which we determine the optimal transfer function of the system's linear 
portion: 


®; (p) 
® (2) =F —=0.65 Dip) ° 





where 0,65 is the value of coefficient k, for v = 0.06. 
For purposes of comparison, we give the value of w, (t) for vy = 0 (the problem of Zadeh and Ragazzini): 


w; (t) = 0.300 — 0.039¢ + 0.3398 (t) — 0.1298 (t — T). 


In conclusion, the author wishes to express her gratitude to V. S, Pugachey and Yu, P, Leonov for the aid 
they have rendered in carrying out the present work, 


















LITERATURE CITED 
[1) L, A, Zadeh, Optimum nonlinear filters, J. Appl. Phys. vol. 24, No, 4 (1953). 


(2) P. L Kuznetsov, R. L, Stratanovich and V. L. Techonov, “The passage of random functions through 
nonlinear systems,” [In Russian] Automation and Remote Control (USSR) 16, 3 (1954), 


[3] G. C, Newton, Compensation of feedback-control systems subject to saturation. J, Franklin Inst. 
vol, 254, No, 4, 5 (1952). 


[4] L E, Kazakov, “Approximate probabilistic analysis of the accuracy of operation of essentially non- 
linear systems,” [In Russian]. Automation and Remote Control (USSR) 17, 5 (1956). 


[5] V. V. Solodovnikov and P, S, Matveev, “The synthesis of correcting devices for servo systems with 
noise present with given requirements on dynamic accuracy,” [In Russian) Automation and Remote Control (USSR) 
16, 3 (1955). 


[6] J. H. Laning and R, H, Battin, Random Processes in Automatic Contro) [Russian translation) IL (1958), 


(7] V. S. Pugachey, The Theory of Random Functions and its Applications to Automatic Control Problems 
[In Russian] Gostekhizdat (1957), 


(8) L, A. Zadeh and J, R, Ragazzini, Extension of Wierner*s theory of prediction, J. Appl. Phys. vol. 21, 
No, 7 (1950). 


Received May 18, 1959 








OPTIMAL CONTROL OF AN OBJECT WITH TWO 
CONTROLLING STIMULI®* 


E, A. Rozenman 


(Moscow) 


The problem considered is that of the optimal transient response in an object 
at whose input the product of two independent controlling functions acts, namely, 
an electric motor controlled by the independent variation of the armature current 
and the excitation current, It is shown that the extremals are combinations of 5- 
functions, 


In recent times a great deal of attention has been given to the question of optimal control of regulated ob- 
jects and, in particular, to electric motors used as executive organs in various automatic devices, There are 
works [1-4] devoted to the determination of the most advantageous law of armature current variation for which 
the duration of the transient response is minimal, With this, the excitation current is assumed to be invariant, 


A significant portion of the results of the works cited can be obtained on the basis of the general theory 
of optimal control, as developed in (5, 6, 7]. In the present paper, as in [8], we consider the problem, of 
practical interest, of the optimal transient responses in electric motors whose control is implemented by the 
simultaneous and independent variation of the excitation and the armature currents, This problem is not im- 
mediately contained in the general theory mentioned above [5, 6, 7] since it leads to extrema from the class of 
generalized functions. 


Thus, we consider an electric motor wherein, as controlling stimuli, the armature current and the ex- 
citation current are used, With the static load and saturation ignored, and with the proper normalization, the 
equations of motion can be written as ; 





dz dy 
a). Uae. (1) 


where x is displacement, y is velocity and u and y are proportional, respectively, to the armature current and 
the excitation current, We consider the transition from the initial state (x = xX», y = yg) to the given state 


(X ™ Xho Y * Yu) 
We shall haye the following two problems [8] in mind, 


1, To so determine u (t) and y (t) that, for a fixed duration t, of the process, the energy dissipation 9, 
will be minimal, i,e., 


t. 
= \(u + v*) dt + min. (2) 


*Presented April 15, 1959 at the all-Moscow seminar on nonlinear problems of automatic control theory, held 
at the IAT AN SSSR, 








2. To so determine u(t) and v(t) that the duration of the transient response will be a minimum for a 
bounded dissipation of energy, i.¢., so that tg + min for 


YS dk. (3) 
Both problems, although having different practical meaning, lead to the same extrema, 
We will consider initially the first of these problems, We write the boundary conditions in integral form 


t, t 
Yr—Yo= | unde, Lk — Ly — Yoto = ie t) uvdt. (4) 


0 0 
These may be considered as isoperimetric relationships which limit the choice of the function u (t) and v (t). 


We now set up the auxiliary functional: 


daw J(u + v) + hyuv + dg (tp — 2) ur) dt, (5) 


0 
where ) 4 and \¢y are constants, On the basis of the Euler equations, we find from (5) that 
2u + yw +g(fg,—t)v=0, 20 +rAuw+,(t, —t)u =0. (6) 


System (6) of equations is homogeneous and has the unique solution u = v = 0 which can not be extermal, 
We replace the zeroes in the right members of (6) by 4 € +0. System (5) can then be satisfied for 


u=v+0, (Ta) 
if the signs of the right members are the same, and for 
u=—vU- 0, (Tb) 


if the signs of the right members are different. 





The validity of this is attested to by the fact that the functional 9, depends in an identical fashion on the 
controlling functions u(t) and v(t) (Cf, (1) and (2)). 


However, the functions u (t) and v (t) may not be defined by (6), It is possible to assume that the extrema, 
if they exist, lie on the boundaries of the domains of the given functions u and v and, consequently, go to in~ 
finity, since no bounds were placed on the instantaneous values of the controlling functions, Condition (2) re- 
quires, however, that the integral 9, be finite and minimal, If the assumption just made is valid, then u (t) and 
v (t) can be functions of the 5-function type, i.e,, for all values of t > 0, u (t) = v(t) = 0 with the exception 
of one or two points where these functions may take the values +o or ~ a, We now show that such a law of 
variation of u(t) and v (t) is in fact optimal, 


We consider the simple case of the transition of the system from one equilibrium state to another, Le., 
we set x (0) = y(0)=0 and y, = 0. In this case, the transient response is symmetric with respect to time 
t =t,/2. We shall seek extermals in the class of binomial functions, and shall set 


u(t) = v(t) =e(1— 2)” for 0<t<, 


u(t)=—v(t) for B<t<ty. 








By integrating, we find that 





y(t) =49(t) = Nery el ea ak 


(9) 


ct te 2t \2m-+2 to 
a(t) = mat t('—reetrnt—-z) —TeeFH!: 


If we take into account that, 


= 
—>§% 
—~4@ 














ay s(t) 
wit) 














Fig, 1, 


for t = te/ 2, x = x, /2 and $= 9,/ 2, we find that 





in * AP... & ad. ty \2 
%= sare te =a 5a(F): (10) 
Whence 
4x, 2m+2 
~= eit ° (11) 


An analogous relationship can be obtained by considering the time 
4z, 2 


09, 
interval ty/2< t < tg. It is easily shown that —=— i, (@m+i?~< 0 


for any m. Consequently, as moo, the magnitude of 9 tends to a least 
value: 


li 42, a 
im 3, = = =% min, (12) 


m-—+oo 


As m +o, the controlling functions u (t) and v (t) reduce to 6- functions, 
Indeed, for t = 0, the binomial in expression (8) equals one, but the initial 
value of 


c= —Vax(m +1)—co as m-—>oo. 


For t > 0 and m -+ @, the binomial tends to zero and, by eliminating 
the indeterminate expression w-0, we can convince ourselves that u (t) and 
v (t) tend to zero, 


For the functions y (t) and x (t) we obtain, correspondingly, y (t) = x,/t) = const and x (t) = xyt/ ts (Fig. 1). 


We can obtain analogously for the interval tg/2 < t < ts that, during the entire interval, u(t) = v(t) =0 
except for the points t - ts, where u (t) + 0, and v (t) +o, 


We note that, in the example considered, the choice of the class of the functions u (t) and v (t) is without 


value, We will be led to the same 


result if we seek the extrema, for example, in the class of exponential functions 


u(t) = ae~™ (13) 


or step~functions (Fig. 2). In the latter case, the least value 9 pin is attained for r +0 and a +a, 


The second problem formulated above, namely, the problem of finding the least duration t, of the transient 


response with a bounded value 9, 


J,= 





= 9, is related to the minimization of a functional of the following form: 


t, 
2 VPs (u® + v%) + dguv + dg (t, — t) uv] dt. (14) 


0 
































Analogously to the previous case, we find that dt, /@m < 0 for any m, and 
limt, = —* = 
a int, = —— = 
Alpe t 
—— tf, ——_> The extrema in this case also consist of 6~ functions, 
we) The expression 
pm 4 tominSomin= 42% (16) 
i : 
0 } t relates the least duration of the process with the least dissipation of energy. 
| The optimal control law obtained can not, obviously, be realized, since an 
infinite controlling stimulus can not be implemented, Therefore, for the deter- 
Fig. 2. mination of a practically applicable optimal process, there should be imposed on 
at least one of the controlling functians, in addition to condition (3), a modulus 
limitation 
|u| <um. (17) 


For u (t) = Upp = const, one obtains a linear optimal law of variation for v(t)(Cf, (3), (4)): 


v(t) = »(1—<). (18) 


It can be shown that, in the problem considered with the limitation expressed in (17), the optimal transient 
response is defined by the expressions: 


v(t) = %(1— +) (0<t<t,), 


u(t) = Um = const 


v(t) =u(t)=0 (tg <t<to—t,), 
tp—t 


v(t)= —r(1—; = ) 
o—"a —t.<t<t 
u(t) = Um = const adliaedatlitad 





The magnitude of tg is so defined that integral (2) will be minimal, The magnitude of vs is determined 
by the boundary conditions, 


In the problem considered, both controlling functions, u (t) and v (t), were independent, and the minimi- 
zation of functional (2) was determined by the choice of these functions, Another formulation of the problem is 
possible, wherein one of the functions, u (t) for example, is chosen beforehand on the basis of some consideration 
of another, Then, functional (2) depends only on the one argument vy (t), and the extremal is defined by (6): 


v(t) =—ruiy—4 (Hu, (19) 


and the multipliers \, and \, are defined on the basis of the boundary conditions, 
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A METHOD OF SYNTHESIZING LINEAR SAMPLED-DATA 
AUTOMATIC CONTROL SYSTEMS TO ACCORD 
WITH DYNAMIC CRITERIA 


L. N. Volgin 


(Moscow) 


A method of synthesizing sampled-data systems is proposed, The transfer 
function of a system is written as the product of the so-called “minimal polynomial," 
to which the given portion of the system is reduced, by an artificial portion which is 
introduced in order to obtain the desired quality. The method presented permits the 
use of a digital computer for the design of linear sampled-data automatic control 
systems, 


The use of sampled-data correcting filters permits significant improvements in the quality of automatic 
control systems, The analysis of sampled-data control systems is widely covered in the literature (Cf, for 
example, [1-5)). Works [6-8] are devoted to the problems of synthesizing optimal sampled-data® systems in 
accordance with dynamic criteria, Work [6] emphasizes an important factor which must be taken into account 
in synthesis, namely, the attempt to reduce the zeroes of the desired system's transfer function which are located 
in the unit circle on the plane of the complex variable z disturbs the *roughness® of the system, The existence 
of such nonreduced zeroes leads to the formation of a *minimal polynomial." The present work gives a method 
of synthesizing sampled~data servo systems which have stable continuous portions and which possess m*th-order 
astatism and minimum control times, In [7] there was given a clear presentation of the idea of factoring the 
polynominals which enter into the transfer function's expression, Here there is presented a purely analytical 
method of designing a system in accordance with the criterion of minimum mean square error, In [8] the design 
of sampled-data control systems is considered as a problem in the theory of dynamic programming, The authors, 
using a matrix representation of the system and the signals, compute the parameters of a linear sampled-data 
correcting device which provides a minimization of the square quality criteria, 


In the present work we present a method which allows one to synthesize optimal linear systems by using 
the techniques of nonlinear programming. This method permits one to obtain systems which satisfy the “rough- 
ness" conditions and the given requirements on process quality, 


1. Posing of the Problem 





The transfer function  (z) of a single-loop linear sampled-data system is 


N (2) W (2) 
L (2) + M (2) W (2)’ (1) 





® (2) = 
where z is the lag operator, related to the differentiation operator p and the pulse repetition time T by the 


*A digital system with a large word-length and operating in accordance with a linear program can also be con- 
sidered as a linear sampled-data system, 
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relationship z = exp (~ pT), N (z), L (z) and M (z) are polynomials in z with real coefficients which depend on 
the system parameters, and W (z) is the transfer function of the sampled-data filter, and is the ratio of two 
polynomials; 


C (2) 
D\2)° (2) 


W (z) = 


The sampled-data"s transfer function  (z) is a rational- fractional function of z, The locus of the zeroes 
and poles of © on the z plane is a complete characteristic of the system properties, 


The introduction of the sampled-data filter in the system changes the number and locus of the initial 
system's singular points, allows its poles to be taken out to infinity, liquidates its zeroes by reduction and intro- 
duces artificial singular points, The interdiction, imposed by the conditions for maintaining system “roughness,” 
on reducing the zeroes placed in the circle |z| = 1 limits the possibility of completely changing the properties 
of the original system, The fundamental distinction between singular points lying on the different sides of the 
unit circle necessitates the factoring of all polynomials which enter into the expression for the system's transfer 
function, The factoring operation consists of decomposing each polynomial, A, into the product of two factors: 


A(z) = A* (2) A (2), (3) 


whose roots are distributed on different sides of the unit circle, The factor A* is formed of the roots for which 
|z| 21, and A’ of the roots for which |z| < 1. In the case of an m-fold root for which |z]| =1, one root can 
be assigned to A* and the others to A™.* 


After all the polynomials entering into @ have been factored, it is necessary to choose the free polynomials 
in such manner that the system's properties are changed in the right direction. This operation is carried out by 
the algorithm described in section 2 which follows, 


2. Algorithm for the Synthesis of Linear Sampled-Data Automatic Control Systems 





The zeroes of the original system's transfer function which cannot be eliminated without disturbing the 
system’s “roughness” form the “minimal polynomial," 


Thus, the system's transfer function can be brought to the form 
D(z) = D, (2) F (2), (4) 
where 4, (z) is the minimal polynomial and F (z) is the artificial portion, a rational- fractional function 


F(z) = ~ (5) 





The algorithm in accordance with which one chooses the transfer function of the sampled-data filter, 
W (z), can be written in the form of the following set of rules. 


1, We factor the polynomials N (z), L (z) and M (z}: 
N=N*N-, L=L*l, M=M*tM. 


2, We decompose the transfer function W into two parts; W = W,W* and, with the choice of W*, we carry 
out the possible reductions which will not disturb the system's “roughness,” and we also introduce an artificial 
zero, A, 


3, By the choice of Wy = Cy/ Dy we reduce the remaining zeroes, and also introduce the artificial pole B, 
This leads to a polynomial equation which is solvable by some polynomials Cy and Dy of minimal degree 
(cf,, the Appendix), 





*One can easily convince oneself that reducing (eliminating) a multiple root for which |z| = 1 would disturb 
the system's * roughness,” 









































4, So as not to infringe the condition of physical realizability of the sampled~data filter, W (0) < a, we 
fix certain coefficients of the low powers of z in Cy and Dg and, in order for the system of equations for Cy, and } 
D, to remain consistent, we increase the degrees of polynomials Cy, and Dg. ! 


5. Starting from the quality criteria, we choose the artificial portion, The quality criteria can either be 
formulated analytically or given in the form of a system of requirements on the transient response, The method 
of presenting the quality criteria affects only the method of seeking an optimum which will be used for pro- 
gramming, and does not affect the method of presenting the initial data for setting up the program, 


Since, in the general case, the minimal polynomial depends on the artificial portion, one must, to choose 
the optimal system, solve a joint system of equations for finding the minimal polynomial and estimate, 


In the case when the minimal polynomial is independent of the artificial portion, the calculations are 
essentially simplified, 


3. Illustration of the Method as Applied to the Simplest Structures 








As illustrations of the algorithm presented, we give the computations for two standard sampled=data auto- 
matic control systems, 


wv }+{eh System 1, The block schematic is shown in Fig. 1, where G (z) 
% is the system's continuous portion and W (z) is the sampled-data 


correcting filter, 








The system's transfer function has the form: 


W (2) G (z) 


®(2)= 1 + W (z) Giz) * 





We set G (z) = P (z)/Q (z), where P (z) and Q (z) are polynomials 
in Zz. 





We rewrite the system‘s transfer function in the form 


Fig. 2, @ — — 2) Pe) 
~ Q+ W (s) P(s)* 





The further calculations are carried out in the following order: 
1) we decompose P and Q into factors: 
P=P'P, Q@=@'C; 
2) we choose W* = AQ*/ p*, and shorten by Q*/ pt; there remains 


AP-C, 


fy 0-D,+A-P-C, * 





3) by equating the denominator to B, we obtain the polynomial equation ® 
QD,+ APC, =B8. 
As the result, we have 
age 
@ = P Com . 
It follows from this that the minimum polynomial equals 


D, = PC, 





*On the solution of the polynomial, cf, the Appendix, 








Shortening of Cy is impossible, since it would lead to the equation 


QD, + (AP” — B)G, = 0, 
the solution of which is 


Cm @, Dom APB. 
Thus, shortening of Cy, would disturb the system's “roughness, 


System 2, The block schematic is shown in Fig. 2, where W, (z) and W, (z) are the system's continuous 
portions and W (z) is the sampled-data correcting filter. 


The representation of the transient response has the form: 








G (2) 
7? 1+ H(z)Wi(z) ’ 
where 
G (2) = W,X (2) = +, H(z) =W,W, (2) = 
1 Qs,’ (z) 1 2 (2) SiS ° 


The polynomials S;, Sy and Q are the denominators of the z-images of Wy, Wz, and X. After substitution, 
we obtain 


PSs, 
QS:S:+QRW  * 


The further computations are carried out in the following order: 


Y= 





1) we carry out the decompositions: 
P=P*P, Q@=Q2O, S=SS, 5,=S)S;, R=RR; 
with the condition imposed on the disturbance, 
Q =1; 
2) we choose W* = = /(QR*A); after shortening there remains 


PS ADg 
RS,S, AD, + R-Cy 





3) we obtain Cg and Dy as solutions of the following polynomial equation: * 
Q5,5S; AD, + RC, = P*B. 
The representation of the transient response takes the form: 
- A 
The necessary quality can be obtained by the choice of polynomials A (z) and B (z). 


4. Choice of the Artificial Part Providing Optimal Quality 





Computing the number and locus of the artificial zeroes and poles which are necessary to attain a given 
*On the solution of the polynomial equation, cf, the Appendix, 
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quality is, in the general case, a complicated task which can only be solved by searching on a digital computer, 
However, in the simplest cases, the computations can be carried out analytically throughout, 


We consider the case when the minima! polynomial is a first-degree polynomial which is independent of 


the artifical part: * 


Y,=1+ az, jaj>1. 


We introduce one artificial zero; ** 


A =1—az, |o| <1. 


We obtain the resulting representation of the transient response 


Y =Y,A = (1+ az)(1— oz) = 1+ (a —o)z—aoz?, 


The criterion 
co 
J= 2 y? 
in the given case gives 


J =1+ («—o)* 4+ a%o*, 
The optimal placement of the artificial pole is defined by the condition 





with which 


Ji+ee 
—itrya° 
The gain in comparison with the case when o = 0 is 


which is the greater, the larger is a, 


A more powerful means of improving the transient response is the introduction of artificial poles, 


For B= 1~ oz, we get 


ee ee + (a—o)z+ 0(¢—a)2* + o* (a —<) 2+ ... 


B i—oz 





*An example of this case would be System 1 with a stable continuous portion with one zero for |z| < 1, 
**The introduction of the artificial zero is carried out by choosing the transfer function of the sampled-data 


filter in accordance with the method presented in Section 3 above, 





co 
The estimate by the criterion J = >) y? equals 
{=o 





(a —o)® 
J=1+5—4 
Optimatization gives 
{ 
o=-— s,s 
a 


With this, the gain is unity. 


5. Example of the Design of a Sampled-Data Automatic Control System 





We consider the system represented in Fig. 2, the transfer functions of whose continuous portions equal: 


p+05 jc? 4 
Wi()= FF 0Sp+0* War) =—p— 5-7: 








Let the disturbance applied to the system have the form of a jump, Le., X (p)=1/p. Wy, (p) corresponds 
to the stable oscillations of the system, The poles of W, (p) are: 


Ppi=—0.4+4 70.3, pp =—0.4—]-0.3. 


The function W, (p) contains fixing and unstable elements, 


By setting up the functions C (s) and H (s) and making the transition to the variable z (by well-known rules, 
presented in [5] for example), we get 


—2+142 
G= a) — 1.282 + 0.452%) * 
0.6 2 (4 — 0.64 z) (1 + 1,3 2) 
H = G—3.7 1 — 1.282 + 0.45 24) * 








We carry out the factorings: 


Pt =—2+4 1.42, P- = 1. 
Q+=i1—z, Q- =i. 

R+ = 1— 0.642, R- = 0.62 (1+ 1.33), 
S| =1—1.2824+0.452%, ST =1. 

SS =1. S, =1—2.78. 


To improve the transient response, we introduce one artificial pole: 


B=i—oz, A=1, jo| <1. 


We choose 


1 
(i — 2) (1 —0.643) * 





Wt = 


The polynomial equation has the form: 


LD, + MCy=N, 
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L =1—524+ 3.129 + 1.3 2% + 1,4 s#, 
M = 0.62 + 0 78 z?, 
N =—2 + (1.4 + 2c) s—1,4 2%, 


We write C, and Dg, in the form 


Co = Co + C32 + cge® + cgs®, 
Dz = d, che dys. 


By expanding the polynomial equation in an algebraic system and introducing the homogeneous notation 


@ =}, y= 6, Te=Co, Me= Cy, Me Cq, Be = Cy, 
we obtain the system of limitations 


2, — 225 + 0.6 25 = — 8.6, 
5a, + 1.42, + 0.78 2, + 0.6 2, = 6.2, 
3.1 2, + 0.78 2, + 0.6 2, = 2.6, 
1.3.2, + 0.78 2, +- 0.6 2g = 2.16, 
1.42, + 0.782, = 0, 
—1< 2g <i. 


The representation of the transient response has the form: 


—2 + (5.4 + 2) s— 2.7 2, 2* 
1—2gz 





Y= 


bos] 
Wetake J =) yj. as the quality criterion, We compute this estimate; 
i=o 





32.2 + 10,8 2 — 21.6 zg + 8.3 2 — 33.2 ayy + 10.8 223 — 5.427 xy 


J= 
1— 2 


The technique of nonlinear programming [9] can be used to find the minimum of this function given the 
limitations imposed, 


APPENDIX 
A polynomial equation of the form 
LD, + MC, =N 


is expanded in a system of algebraic equations 


lgdg + Moly = No 
Iydq + lg dy + mylg + Mots = ny 
Igdg + by dy + lode + macy + myc, + mocg = ng, 


which has a solution if the following conditions hold 
S>p—i, x>A—i, voA+b—=p+x, 


where X, , 5, K, vy are the degrees of the polynomials, respectively, L, M, Dg, Cy and N, 








The condition of physical realizability of the sampled-data filter requires that dy # 0. If mg = 0, it is then 
necessary to choose N such that ny ¥ 0. 


The polynomial equation is equivalent to the congruence LD, ® N (mod N), and Euclid's algorithm can 
also be used for its solution. 
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AC SERVO SYSTEMS IN WHICH THE SIGNAL DEPENDS ON 
THE ERROR AND ITS DERIVATIVE 


N. P. Vlasov 
(Gor’ kil) 


AC systems are considered, Under the assumption that the signal depends on 
the error and on the speed with which the error changes, we find the transmission 
function for a linear four-terminal network and for servo systems which are stabilized 
by differentiating four-terminal networks and velocity feedback, 


In work [1) it was assumed that the function at the output of the modulating device for a servo system can 
be presented in the form 


aU 
v(t) = x (t) == cos (wot + 9). (1) 
Here, 


x (t) = 0, (t) — 4, (t) (2) 
is the error and 0, (t) and 6, (t) are the functions at the input and output respectively of the servo system, The 


function v (t) is the oscillations of the carrier frequency ws modulated by the error, 


Actually, at the output of a modulating device which is implemented, for example, by means of two 
selsyns, with the assumption that the following conditions hold 


6, (t) << @,, 4, (t) <p, (3) 


there appears, as can be shown [2], the function 


ux (t) = = [= (t) cos (@gt + 9) + . sin (@>t + ®) | ; (4) 


In the present paper, on the basis of the methods developed in [1], we consider servo systems in which the 
transmission function is specified for the four-terminal network connected in the control loop, as well as systems 
which are stabilized by a velocity feedback, with the condition that function (4) appears at the modulating de- 
vice’s output. 


1, We consider initially a linear four-terminal network V with transmission function Y (p), connected 
between modulator M and ideal phase detector D [1] (Fig. 1). The Laplace transform of the function uy (t) will 
be 


4 0 





m 


Ux(P) = > Ge 


i [z (p — jw,) e'* — x (p + ja) e~*]. (5) 






































# x The transform of the function u (t) at the phase detector's 
Ml uU;| ¥ | uD input is 
u (p) = ux (p) Y (p). (6) 
rig. 1. (p (p) ¥ (p) 


Let the phase detector be ideal, and let it implement the 
P multiplication of the functions u (t) and cos (wet + ¥), Consequently, 


PTHY ul y Eo at the output of the phase detector we shall have 
1 4% 7 [fe 














up(t) = u (t) cos (gt + 4). (1) 
Fig. 2, The transform of function (7) is 
Up(p) = = [4 (p — jw) e** + u(p + jmp) e—*}. (8) 


By substituting the value of u (p) from (6) in this expression we get 





4 9U fa) , genereeenniene 
TD = 452 [PHM PHF Foy ete» — 


1% 


_ {pth Y (p — j,) eto) z(p)+ 


]® 





IP — ¥ (p — ja) eX? x (p — j +2) — 





— PAI Y (pF ja.) eet) Z (pF 7 -2a)} 
If we discard from this expression those terms which correspond to oscillations modulated at twice the 


carrier frequency, which will not be passed by the low-frequency filter [1], we obtain the following transmission 
function for the chain of links between the input link and the phase detector: 





77 Un?) 1 OU», I a a 
TO) = “BE = SP [SEY EF fojeto — 


— (P—I%) F (p — 70) e-He-0) | (9) 


1% 


If the four-terminal network is inertialess and Y (p) = ky then, for Yp (p) we shall have 


a au , 
Pw) = "+ G2 | cos (p—4) + Z sin (p—4)] ao) 
Thus, by selecting the phase angle ¥ of the reference voltage applied to the phase detector, one can not 
only affect the gain, but can also obtain a component which {s proportional to the velocity of the error. 


2, We now consider an ac servo system with a single-phase commutator motor as the servomotor, con- 
trolled on the armature side, 


In this case, the phase detector is the servomotor and, for the transform of the function at the phase de- 
tector’s input, i.e., for the voltage on the servomotor's armature terminals, we take expression (6) (Fig. 2). 


With certain simplifying assumptions, formulated in [1], the servomotor equations can be written as: 


the equation of the armature circuit 


u(t) = ue(t) = ir + LE + Cig; (11) 
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the equation of armature rotation 





















































~ ,)U, } d*6, er, * 3 . d6, 
LA t Y pe J | Me J ar = Cigi — Me = Cei— ks}. (12) 
. 3 ©. " 
5 21 44 J The nomenclature used in Eqs. (11) and (12) was intro- 
ce Dene eae d duced in [1). 
Fig. 3. For the current in the servomotor’s excitation winding 
we take the expression 
lq 
ic = Im cos (Wot . 13 
20-089 Y m= c m COS (Wot + >) (13) 
pare ee ee / We now take the transforms of Eqs, (11) and (12), taking 
Eq. (13) into account. We find 6, (p) from the system ob- 
Fig. 4. tained, taking Eqs, (5) and (6) into account, having first dis- 


carded the terms corresponding to oscillations modulating 
twice the carrier frequency, The result is 

















1 Clm 0m { (p+ fod F PETA) ,-y0-»)_ P= ied TOTS) ro- } op 
4G ja Oz i _ (+—9) 
8, (p) —4_ioe 92 UO +L tion * PAT weg zp) 
. { (Cl,,)* (r + Lp) . 
TP ++ 2 e+ Dp + (Loo J? 


If we go to the argument p/w, = 8, we will have, for the transmission function of the open-loop system 


Bay _ 1 VAT + N+ NVC —HPjOt+MTas—N+ NV ETD 
z(s) @ Ty8 (hen (T m8 + 1) ((T 98 + 1) 4 79) + 78 +4) ' 











(14) 


Here, 5 =¥— gy The values of the constants Ta, Tp, km and Ty are given in [1). 


We now consider the systems which is stabilized by an ac differentiating four-terminal network [3, 4) with 
transmission function 





Y (s) = ko; (8), (15) 


where 





s* + ks +1 
Yi (8) = yet * (16) 


For the transmission function of the open-loop system we get 


is) ai A (8) {cos 8 — [T,, (s* + 1) + 8] sin’ 8) + {[7,, (s* + 1) + 8] cos 8 + sin 8) (kg — h,) s* 








17 

za) T 43 (hey (T gy + 1) ((T gs + 1) + 73] + Tas + 1) B (a) vis 
Here, 

A (s) = s? (s* + 4) + (hy + hg) (s* + 2) 8 + Aik (8? + 4), (18) 

B(s) = 2 [s? (s* + 4) + 2ky(s* + 2)s + 43(s? + 4). (19) 
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We now consider a system stabilized by velocity feedback (Fig, 3). In the control loop let there be four- 
terminal network 1, with transmission function ¥ (Pp) and in the feedback path, let there be four-terminal net- 
work 2, with transmission function Y_(p). The voltage u, (t) is applied to the input of four-terminal network 2 
from a bridge circuit in the motor’s armature circuit. For this voltage we take the following expression: 








u, (t) = hci S. (20) 
The transform, u, (p), will be 
us (p) = 5h Clm [(p — jo) 4, (p — ja) ef + 
+ (p + jo) % (p + jm) e~ 54). (21) 
The transform, Up (P)> of the voltage at the output of four-terminal network 2 will be 
Uo (p) = u, (p) Ys (Pp). (22) 


Since the transform of the voltage u, (t) at the differential output is defined by expression (5), the trans- 
form of the voltage at the input of four-terminal network 1 will be 





Uy, (p) = Ux (p) — Uo (P)- (23) 
Here, Up (P) is given by Eq. (22), 


In this case, the transform of the voltage on the motor*s armature winding will be 





uw, (P) = Win (p) ¥ (p. (24) 


By again using Eqs, (11) and (12), taking the transforms of these equations, and also by taking Eqs. (5) and 
(21)-(24) into account, we obtain the following transmission function of the open-loop servo system, written for 
the argument s: 


Be) _ 1 ie— MIT, (8+) +11 @=7) e* —7 (+f IT, @—D +N YETI 














Tr — (25) 
a(s) 2 T 48 mn (T m8 +4) [(T 98 + 1)® + 72] + 7,8 + 1+ kA (s)} 
Here, 
H(®) = +(%,06—NY¥@—s)iTa(s+ i) +11+ 
+¥o(s+ ))¥ (s+ /)(Ta(s—/) + 1)}- (26) 


3. We now consider a servo system with a single-phase commutator motor, controlled on the side of the 
excitation winding (Fig. 4). The phase detector is once again the servomotor, but here expression (6) is the 
transform of the voltage supplied to the excitation winding. With the simplifying assumptions made in [1, 5), 
the servomotor equations can be written in the following way: 


the equation for the excitation winding circuit 


; di 
u, (t) =ir+ Ls; 









































the equation for the armature circuit 














u(t) =iere + Lege +ci %; (28) 














the equation for the armature rotation 


d0. 


> Cita — M e=Ciig—k—=- = (29) 


dt? 





be The nomenclature used in Eqs, (27)-(29) is given in 
(1). 
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We make the assumption that the last term in the 
right member of Eq, (28) is small; this condition holds 
for sufficiently small servomotor rotational speeds, Then, 
with a sinusoidal voltage u (t) applied to the armature 
terminals, the armature current ig will also be sinusoidal, 








We set 


ig = Im C08 (Wot + 9). (30) 


We now take the transforms of Eqs, (27) and (29), keeping Eq. (30) in mind; we find 6, (p) from the system 
obtained with Eqs, (5) and (6) taken into account, having first discarded the terms corresponding to oscillations 
modulating the doubled carrier frequency, The result is 








_ ae = {Lip— joe) Ye — fa) nw _ UP + fod FT 0} oo 5 
m Gz _\_ [F+L(p— joo)| [F+Z(p + fool <a 
4m (J p* + kp) 








6, (p) = 


By going to the argument s = p/w», we will have, for the transmission function of the open-loop system, | 





(8 1 (s—DIT +) +1 Y 6A ®F —j (s+ /) IT. —)N+UY EF) | 
a: 2 KT 8 (Ty + 1) (Ts + 1)*+ 72] ; 


a= 





Here, the constants T,, Tm, ky. Ty, and 6 have the same values as in Eq,.(14), 


We next consider a system stabilized by a differentiating four-terminal network with transmission function 
(15). For the transmission function of the open-loop system we obtain 








63 (8) k A (s) cos 8 — [T,, (s*+-1) + s] sin 8 +({[7,, (s*+-1) + s] cos 5 + sin 6}(kg—k;,)s* 
i kT ',8(T m8 +4) {(T 98 + 1)? + T2] Bs) 











Here, the functions A (s) and B (s) are defined by expressions (18) and (19), Asan example of a system 
stabilized by a velocity feedback, we consider a system with an asynchronous tachometer generator (Fig, 5). It 
can be shown that, if a bridge scheme is used for velocity feedback (Fig. 6), the systern will be described by a 
set of nonlinear differential equations with periodic coefficients, for which the method developed in [1) is in- 
applicable, 


For the voltage taken off the tachometer generator, u, (t), we take the expression 


U 
u, (t) = ke = COS (Wot + «) . . 








For the transform, u, (p), we obtain 


— k ee 3 
Us (P) = 5 Um I(p — joo) 5 (p — jg) e'# 4 
+ (p + jo) 92 (p + J) e~i2). (31) 
By again using Eqs, (27)-(30), taking their transforms, and by also taking Eqs. (5), (31) and (22)-(24) into 


account, we obtain the following expression in the argument s for the transmission function of the open-loop 
system: 


G(s) _ j *—/) ¥ (@— +7 se +7] e¥—*) _j (947) ¥(s+j)[1+T7, (s— je 
# (s) Ts {(T18 + 1) ((1 + 7,8)? + T2] + Tok (s)} 





Here, 


E(s) 





5 FE—DV6— Ut + Tals + few + 


ll 





+¥@FAY CF IM + To(s— fee). 


The values of the constants T;, T; and Ty are given in 
20:00 TARE: (1. | 
/ 














P th 4, We now consider a servo system with a two-phase 
9, 5% © | & | asynchronous motor as the phase detector (Fig. 7). With 
bn cacmeswess J pc certain simplifying assumptions, formulated in [1], we can 
Fig. 7. write the servomotor’s equations of motion in the following 
way: 


the equation of rotor rotation 





dq, < 
<a Ja = M (t)— MM, (tt) ss 
£7 1 ; d8 
Cae ee Tae ee Se = + [Petts — Gelte — FP (42 + 42) | — RSP (32) 
Fig. 8. 
© the equation of the control winding 


dy 


u(t) = a: (33) 
the equation of the second stator winding 
dhe 
Fr = M2 (t) = Usmsin (wot + 4)- (34) 


The nomenclature used in Eqs, (31)-(34) was given in [1). 


By assuming that the amplitude of the controlling voltage is small compared with the amplitude of voltage 
u, and, consequently, that v« vw, we discard the term v from Eq, (31), We now transform Eqs, (31)-(33), 
find 63 (f) from the system obtained with Eqs, (5) and (6) being taken into account, having first discarded the 
terms corresponding to oscillations modulating the doubled carrier frequency, The result is 








ss Yam U mm li (2 — 12) Y (p — jaje¥—*) — (5 + j2) ¥ (p + joo) ie) 


6, (p) = 4rwo Ox (, { yu, 
Pk + 5 r 2) 


By using the argument s = p/w», we obtain the transmission function of the open-loop system 





6 (s) _ j (s— j2) Y(s—7) e*® —j (8s +72) Yetne* 
>) Tas (Tis + 1) 





The values of the constants Ty and T, are given in [1). 


We now consider a servo system stabilized by a differentiating four-terminal network with transmission 
function (15), In this case, the transmission function of the open-loop system is 


BC) _ y, (2+ js) (ACW) + j (ks — hs) st] e* + (2 — js) [A — I (hs — ba) A e* 
a (s) T 28 (Ts + 1) Bis) (s) 








Here, functions A (s) and B(s) are defined by expressions (18) and (19), 


We now consider the servo system stabilized ized by a velocity feedback (Fig, 8), In the control loop is four- 
terminal network 1 with transmission function Y (p) and, in the feedback path, four-terminal network 2 with 
transmission function Y3 (p). Voltage ug (t) is applied to the input of four- terminal network 2 from the bridge 
scheme in the control winding circuit. It can be shown that voltage ug (t) is proportional to the rotational speed 


of the motor's rotor, or to voltage ug (t), but shifted in phase by #/2 with respect to voltage u, (t). We take 
the following expression for voltage u, (t): 





k 
u, (t) = = U om C08 (Wot +- >) 7. (35) 
The transform of ug (t) is 
u(p)= > 5 Vm l(p — ]@o) 9, (p — jw») e%% +- 
+ (p + jo) 9, (p + ju) e~*). (36) 


To find the transform, u, (p), at the input of the control winding, we use formulae (22)-(24), 


Taking the transforms of Eqs, (32)-(34) and taking Eqs. (5), (36) and (22)-(24) into account, we obtain the 
following expression in the argument s for the transmission function of the open-loop system: 


Os (8) _ [i (s — j-2)¥ (s— j) e* — is +4-2) VFA e*) (+1) 
z (8) T'9s [(T 18 + 1) (s* + 1) + gk b (8) 








Here, 





£ (s) = + ls—j-2)(s + N¥C—Nhe—N+ 
+ (s+ 7-2)(s—/)Y¥ (s+ )¥a(s+7)) 





and 
= ieee 
2rkw@y + Uem 
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Thus, the methods presented in [1] can be used to obtain the transfer functions for ac systems whose signals 
depend on the error and on the velocity of the error’s changes, 


In this paper we obtained the transfer functions for linear four-terminal networks and for servo systems sta- 
bilized by differentiating four-terminal networks or by velocity feedback. For this, we considered as phase de- 
tectors both single-phase commutator motors (with control on the side of the armature or on the side of the ex- 
citation winding) and two-phase asynchronous motors. 
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GRAPHICOANALYTIC METHOD OF COMPUTING 
STEADY-STATE SPEED 
OF PNEUMATIC DRIVES 


V. N. Dmitriev and V. I, Chernyshev 


(Moscow) 


A method is considered for computing the steady-state motion of a pneumatic 
piston executive mechanism, A comparison is given of the computed and the ex- 
perimental characteristics of such a mechanism, 


Pneumatic pistons are widely used as the executive mechanisms in pneumoautomatic systems, However, 
very few works, either domestic or foreign, in the literature are devoted to the theory and design of this type of 
mechanism, Because of the compressibility of air and the variability of the air which is charged into, and dis- 
charged from, the cylinder voids in the executive mechanism, the motion of this mechanism is described by a 
system of complicated nonlinear differential equations, The majority of the authors therefore consider this 
problem in a linear approximation [1-3], There are also works [4, 5) which consider the motion of a pneumatic 
piston executive mechanism without feedback by stages, Such a solution to the problem is complicated, 


The method of designing pneumatic piston executive mechanisms which {s presented in the present work 
does not pretend to applicability to an arbitrarily chosen mechanism of this type. However, it can be success- 
fully applied for pneumatic piston executive mechanisms whose pistons have a significant distance of travel, 
which do not thereby develop large accelerations, and whose load is constant. The class of such systems in 
industry is quite large, This is a different type of drive than drives for clamping devices, drives for feed 
mechanisms and terminal organs of regulators, for moving gates and dampers, etc, The pneumatic piston executive 
mechanisms alluded to frequently have either constant, or slightly varying, loads along their entire path of motion, 


1. Unloaded Piston Executive Mechanism 





The schematic of a piston executive mechanism is given in Fig. 1, Supply pressure P, is kept constant, 
pressure P, = 1 atmosphere, It is easily seen that, in the case when there are no loads on the piston and the piston 
has zero mass, the piston executive mechanism is an ordinary pneumatic flow chamber, By bearing in mind that 
the design of valves and other forms of distributing devices always provides equality of the geometric areas of 
the valve apertures, and making the additional simplifying assumption that in such mechanisms the discharge 
coefficients for the apertures for charging and discharging the cylinder voids equal each other, i,e,, that py =p, 
we obtain the equality of the effective areas: yf, = sf = f. 


From the general theory of pneumatic flow chambers [6] we know that if the effective areas of the valving 
apertures are varied in such fashion that their ratio always remains constant and if, moreover, pressures Py and Pg, 
are also unchanged, then the pressure in the intervalye chamber P, will also be constant, 


Bearing in mind the remarks made above, we write the equations describing the motion of an unloaded 
piston executive mechanism, For writing the equations we shall use the approximate expressions for the air dis- 
charge through the valve apertures which differ only slightly in their accuracy from the well-known formulae of 


















































Fig. 1. Schematic of a pneumatic 
piston executive mechanism: Py is 
the supply pressure, Ps is atmospheric 
pressure, f, and f are the geometric 
areas of the apertures, F., is the 
piston area, Py and Py are the pres-~ 
sures on the cylinder surfaces, P; is 
the load on the piston rod, x is the 
piston coordinate, a is the piston 
thickness and 7 is the cylinder 


Sen-Venana Vantsel’ for an adiabatic process, The equations cited 


are obtained by equating the air discharge through the valve apertures 
to the air discharge along the cylinder cross section, * 


Then, for the case of connected modes of flow through the valve 
apertures ** "sub-sub,* we shall have 


for the intake chamber 


P 
1V 76, P, (Py — Pi) = Et Poi a) 


for the discharge chamber 








1V & (Pi—1) = FE pe (2) 


where g is the acceleration due to gravity, R is the gas constant, T is 
the absolute temperature of the air ***,F {s the piston area and x is 
the piston coordinate, 


By eliminating pressure Py from Eqs, (1) and (2) and then sub- 













































aan stituting its value in one of them, we shall have; 
- Vy — Tt a —1 fob}. (3) 
Here, 
~ * — V P,—1 TVD T me 
. In the case of combined modes of flow “sub-super,” we have 

1V Fe PilPo— Pi) = the (4) 
1P1V shn= a *pi (8) 

By carrying out an analogous transformation, we get 
z=. (8) 


It follows from (3) and (6) that the velocity of the piston executive mechanism for subcritical flow through 
the first and second valve apertures of the distributing device depends on the area of the valve aperture cross 
section, on the supply pressure and on the coefficient k and, in the case of subcritical flow through the first valve 
aperture and supercritical flow through the second, depends only on the effective area f of the valve apertures 
and on the coefficient k, and does not depend on the supply pressure Pp». 


We now determine the areas in which Eqs, (3) and (6) are applicable, 


*We shall assume that, prior to the moment under consideration, some motion of the piston had occurred, and 
that the transient response connected with establishing the pressure P,, caused by the given pressure Po, had ter- 
minated, i,e,, we consider the steady-state motion, 

**Here “sub-sub* and further *sub-super” denote, respectively, subcritical flow through the first and second 
valve apertures, and subcritical flow through the first, and supercritical flow through the second, 

***We assume that the transition from the chambers before the valve apertures to the chambers after the valve 
apertures is isothermal, i,e., that the air temperatures in these chambers are all equal, 




















If we solve Eqs, (5) and (4) jointly, we find the following relationship, which is constant and equal to 
Py/P, = 0.8 for subcritical flow through the first valve aperture and supercritical flow through the second, 


We now investigate the function g, which enters into expression (3), 


By differentiating this function with respect to P,) and equating the result to zero, we determine the pressure 
P, for which the function assumes its maximum value, This pressure will equal 2,5 atmospheres, and the value of 
function @, for this pressure will be 0.5. Thus, with P, = 2,5 atmospheres, expression (3) turns out to be identical 
to expression (6) for the case “sub-super,” As the source pressure is increased above the value of 2.5 atmospheres, 
the magnitude of function ¢ will be decreased although, naturally, the combined modes of flow through the 
valve apertures will be “sub-super,” 


Therefore, for all values of P, greater than 2,5 atmospheres, expression (6) should be used or, what amounts 
to the same thing, instead of gy, the value 0,5 should be substituted in (3), The combined modes of flow through 
the valve apertures “super-sub" and “super-super® can not exist. This will be shown below for the more general 
case of a piston executive mechanism with a constant load, 


2. Executive Mechanism with a Constant Load on the Piston 





We now obtain the equations for a piston executive mechanism found under a constant load, 


The weights of the air (Fig. 1) in the left and right chambers of the servocylinder will be, respectively, 
P P 
Qi = ay yr and Qa = (—2—a) Fy er 
By differentiating both equations with respect to time, we obtain the equations of the intake and discharge 


cylinder voids, By adding, in addition, a third equation which expresses the equality of the forces on the piston, 
we will have a system of equations which completely describes the law of motion of the piston rod: 


dQ, __ Fp ¢; ; (7) 
t= ph (Pi + =P), 
F . ° 
£Qs _ 5B [—aP, + (l— 2 — a) Py), (8) 
P 
Pity Pam ye Ae. (9) 


Here, the derivatives dQ,/dt = f¥ (Py, Pg) and dQ, /dt = fH (Pz, Ps) are the mass per second discharges of 
air through the valve apertures of the distributing device, 


By analogy with the unloaded piston executive mechanism, and by also taking into account that the dif- 
ference P,-P, is a constant quantity, we shall assume that the pressures P, and P, are themselves constant,* In 
other words, of all the possible solutions to the system of (7), (8) and (9), we shall choose such for which Py and P, 
are constant, 


Taking into account the assumptions made, we rewrite Eqs, (7), (8) and (9), for example, for the case 
“sub-sub,” in the following form: 





2 F ; 
1V Te Ps(Po— Py) = Pat, (10) 
Rar 
1V pp(Ps—t) = phePy, (11) 





*Here, a3 previously, it is assumed that prior to the moment of consideration of the system's behavior there was 
a certain displacement of the piston under the action of a given supply pressure, and that the transient response 
connected with the establishment of the steady-state pressures Py and P, is terminated, 































ta hkY Po —4y, (13) 


-  » VP,—AP—i 
z=k P,— AP }. (14) 





























It might be possible here, by eliminating pressures Py and P, 
from the system of (10), (11) and (12), to find the equations of motion 
of the piston executive mechanism, However, such a route is a 
complicated one, It is more reasonable to use a graphicoanalytic 
method, This method reduces to the following. From the system 
of (10), (11) and (12), one determines AP as a function of Py and 
the supply pressure Py, After this, one constructs the graphs of this 
"A function and, for given values of P, and AP, one finds Py, sub- 

Is stitution of which in (13) or (14) can give an expression for the 
Fig. 2. velocity of piston motion, 


Gj 


We now find the dependence of AP on Py and Py. For this, we shall replace the pressures, for greater 
generality, by their ratios, namely, 


P14 P, P, a4 


ise "1 = Pp,’ “a, * wre FP, 


Then, in accordance with (10), (11) and (12), and adding (16), we shall have 


1 





Fy = APrs +4, (15) 
r = PsP, (16) 

J 

— “ = FFs. (17) 


If we now eliminate r, and rg and then solve the resulting quadratic equation for AP, we get 


_ y(2(1—r;)—r] Vr: (57; —4) 
4P= 2r (i — ry) = 2(i1—rn) ° (18) 





By considering expression (18) we can easily convince ourselves that the minus sign before the second term 
is meaningless, Indeed, if expression (18) is analyzed with the minus sign before the second term, one can see 
that, with increasing r, and, consequently, with increasing Py, the drop across the piston AP is decreased, which 
contradicts the physics of the phenomenon, By equating the second term to zero we obtain the minimum value 
of r for which the drop across the piston, AP, still has a real value, This magnitude will be r, = 0.8, By its 
definition, the greatest value of r; cannot be greater than unity. Thus, for combined regimes of flow “sub~sub,” 
the magnitude of r varies from 0,8 to 1, 


For joint regimes of flow in the intake and discharge valves of “sub-super,” the system of equations of 
motion of the piston executive mechanism will be: 





: P. 
rau 22-1, (19) 
z= sf, (20) 
P,— P, = AP. (21) 


It is clear from this system of equatims that the velocity of the piston in this case depends neither on the 
load on the piston rod nor on the supply pressure, 
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By solving (19) and (20) jointly, we find that ry = 0.8, 


Thus, pressure Py and, consequently, also pressure P,, does not depend on the velocity of the piston’s motion, 
The joint modes of flow “super-sub” and *super~super® through the valve apertures of the cylinder*s intake and 
outlet voids do not exist, It is possible to convince oneself of this by investigating, for the cases cited, the 
equations of motion of the piston executive mechanism, Then, for a joint regime of flow “super-sub” one ob- 
tains a complex drop on the piston and, for °super-super,” ry = 1. 


Since neither the one nor the other is possible, these modes of flow, with the assumptions made,* are ab- 
sent in piston executive mechanisms, 


Using Eq. (18) with respect to joint regimes of flow “sub-sub,® and also taking into account the fact that, 
for the case “sub-super,” r, = 0,8, we construct the graphs of Fig. 2, 


All curves on the figure are for joint modes of flow “sub-sub," and all modes “sub-super® lie on the 
vertical line passing through the value of ry of 0.8. 


By using these graphs we can, without previous knowledge of the joint modes of flow through the valve 
apertures in the intake and discharge chambers, but knowing only the supply pressure P, and the fall across the 
piston, AP = P,/ F.,, determine the value of. With this value of rm we can, either by formulae (13) or (14), 
compute the piston velocity. 


To simplify even further the computations connected with the determination of the yelocity of motion, x, 
of the piston, we can construct a graph (Fig. 3) whose axis of ordinates represents the function ¢, equal to the 
radical in expressions (13) and (19), and whose axis of abscissas represents supply pressure Pg, The parameter is 
the pressure drop AP, 


The graph was constructed by means of Fig, 2, from which the values of r, were determined, but the 
function g, was computed from the formula @, = \(1/m— 1, All curves on the graph relate to the joint regi- 
mes of flow “sub-sub,” and the horizontal line g, = 0.5 relates to the joint regimes “sub~super.” It is also clear 
from the graph that there will be no piston motion if Py =< (P;/ Fp + 1) atmosphere, 


Knowing the value of ¢, we can easily compute the velocity of the piston of the pneumatic piston exe- 
cutive mechanism from the formula 


3, Experimental Verification of the Method for Determining the Velocity of Piston 





Rod Motion 





To provide a basis for the correctness of the assumptions made, we experimentally determined the char 
acteristics of a loaded piston executive mechanism, The piston executive mechanism investigated had a cylinder 
diameter of 31,6 mm and a total travel of 114 mm, The load was connected to its piston rod, With the motion 
due to gravity taken into account, ** the load was so chosen that the pressure drop on the piston, AP, equalled 


*Here we have kept in mind the equality of the effective areas of the valve apertures and the constancy of the 
load on the piston rod, 

* *For the piston executive mechanism investigated, the force of gravity was about 300 grams, Such a com- 
paratively small gravitational force is explained by the use of rubber packing rings fitted on the piston with a 
significant initial gap. 
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Fig. 4, Oscillogram showing the dependence of the pressures in the cylinder voids 
and of the piston travel on time, d is the diameter of the valve diaphragm, 1 is 
the pressure in the discharge chamber, 2 fs the pressure in the intake chamber, 

Py = 2.1 atmospheres, d = 0.618 mm and x = 4,25 cm/sec, 














t=267sec 7 


£cm/sec 
16 


1 kilogram/ cm’, The piston rod load was implemented by 
a bob suspended by a fine steel wire drawn over light 
rollers, As valve elements we used impedances of the 

type “apertures in fine walls." In the pressure range in- 
vestigated, the discharge coefficient of such valve ele~ 
ments is constant, and equal to 0,8 [8]. 


Two series of experiments were conducted, The first 
series of experiments was carried out in such fashion that 
the initial pressures in the cylinder voids equalled atmos- 
pheric pressure, and the supply pressure P, was applied to 
the chambers by jumps, However, in some case when 
there were large piston accelerations, caused by com- 
atm paratively large valve diaphragm passage-way cross sections 
and large supply pressures, the pressure in the cylinder 
voids, with the lengths chosen for them, could not assume 
Fig, 5. Speed characteristics of the piston its steady~state value during the time of piston travel, 
executive mechanism. 


~~ = &® @ &§ 





sub-sub sub-super 


Since the methodology given in the paper is to be 
used only for computing steady~state piston motion, a 
second series of trials was made, The trials were carried out with an artificially established computed pressure 
in the cylinder voids. 


The pressures, P, and P,, in the cylinder voids, and the piston travel, were taken down on oscillograph 
paper, 
In the course of the experiments, a significant number of oscillograms were collected for various supply 


pressures and different dimensions of the valve diaphragms, The results established that the pressures in the cy- 
linder voids can be taken as constant, not depending on the dimensions of the valve apertures, 


Moreover, the fact was also established that the supply pressure variations in the region “sub-super™ do not, 
in practice, vary the magnitude of the piston velocity in the piston executive mechanism investigated, 


As an illustration, Fig. 4 shown an oscillogram, The parameters with which the experiments were con- 
ducted are also shown there, 


Figure 5 gives graphs constructed from the experimental results which express the dependence of piston rod 
velocity x of the piston executive mechanism on the supply pressure Py for various values of valve diaphragm 
diameters, with a pressure drop on the piston equal to 1 kg/ cm’, For purposes of comparison, the figure also 
shows, via the dotted curves, the computed curves, 


As is obvious from Fig. 5, the maximum error does not exceed 7.6% * 


*The maximum discrepancy between the computed and the experimental data (we are speaking now of the case 
when there were large valve diaphragm passageway cross sections and large supply pressures) when the experi- 
ments were carried out without the establishment of the initial pressures in the cylinder voids was about 30% 
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ON ONE METHOD OF DESIGNING TWO-INPUT 
ELECTRONIC FUNCTION GENERATORS 


V. B. Smolov 


(Leningrad) 


This paper deals with circuits for electronic function generators for simulating 
functions of two arguments, these circuits being constructed on the basis of standard 
operational amplifiers with controlled impedances in the input and feedback circuits, 


Concrete block schematics are given for generators to realize two- argument 
functions which are characteristic of standard analog (simulation) technology. 


1, Modern electronic simulating (analog) devices contain, as their basic standard elements, the so-called 
operational amplifiers, on the basis of which one can implement various mathematical operations on the input 
de voltages, However, in simulating complicated functions of two variables [1-3] 


Jou = © (Uin1, Vino) (1) 


the electrical circuits of the corresponding electronic function generators, using standard operational amplifiers 


and constant impedances, are significantly complicated, leading to reductions in a number of the usage and 
technical indicators applied to analog devices, 


In the present work we consider one of the possible methods of constructing function generators which are 
useful in simulating a number of important particular cases of function (1) with an accuracy, speed and re- 
liability sufficient for practical use, The basic special feature of the generators to be considered is the presence 
in their circuits of standard elements from electronic digital technology by means of which the necessary switch- 
ing of the operational amplifiers’ input and feedback circuits is implemented, 


2, The standard circuits for operational amplifiers [4, 5) provide two types of transformation of functional 
characteristics: 


a) a standard operational amplifier with anode~grid negative feedback 


—— 
S) 23 

2= —T 2 rr.” (2) 
1 25 


b) an operational amplifier with differential inputs 
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per rity) reowry rfz.y~const) In formulae (2) and (3) we haye used the notation: 


= X4k and xj are the input voltages of the amplifiers, ry 
J, 
| > 


| and tj are the impedances in their input circuits and typ 
is the impedance of the feedback circuit, 
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If the amplifier impedances are implemented in 
the form of variable (“controlled”) ohmic impedances 


% 2D whose magnitude depends on the two input voltages 
a 


Fig. 1. ik =Tie(Xix, Yr), oj = Pej (Tej» Yrs), (4) 

















then the output voltages z of the amplifiers are complex 
functions of the input voltages; particular cases of these 
functions include the simulation of transformations of the 
type in (1) which are widespread in practice, 


Thus, the use of controlled impedances with the 
characteristics of the form given in (4) in the block 
schematics of the operational amplifiers permits the 
solution of the problem of designing function generators 
with two inputs, 








3. One can consider three fundamental circuits, 
constructed with controlled ohmic impedances, which 
are used in computing and measuring technology. 








The first of these, shown in Fig. 1, (5) is used to 
obtain a functional impedance which varies stepwise in 
dependence on the input voltage x with constant values 
of voltage y: 








retains its stepwise character, but is gradually shifted 
along the x axis, 





The second fundamental circuit [6) for a con- 
trolled impedance contains linear transformer I which 
converts input voltage y into the digital quantity N which 
ti ——_——-8 is statisticized in counter IL The output buses of the counter’s 
flip-flops, one for each order (bit), excite the correspond- 
ing relays III whose contacts provide the successive 
switching of constant impedances IV (Fig, 2). 


a J r=r(z, y = const). 
-H With a variable value of voltage y, the impedance 
es 
i 
' 
' 
se 














Fig. 3. 


In the general case, the total magnitude of the digitally controlled impedance r can depend functionally 
on voltage y if, instead of a linear transformer of the input voltage to a digital quantity, one uses a nonlinear 
transformer, 


The third circuit, shown in Fig, 3, is made up of constant ohmic impedances IV, connected in parallel in 
a common circuit by means of the contacts of the relays III which are switched by the output voltages of the 
order flip-flops of pulse counter II, 


By combining the circuits considered with their impedances controlled by diodes (Fig, 1) and digitally 
(Figs. 2 and 3) with constant capacitive impedances in various block schematics of operational amplifier con- 
nections, one can provide simulation of various functions of type (1). 


Below, we consider the block schematics of concrete modes of operation of nonlinear two-input transformers, 
constructed of the most widely used operational amplifiers with anode~grid negative feedback, and we provide 
certain information related to the questions of using block schematics with differential inputs, 
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4, Standard operational amplifiers with anode-grid negative feedback, given the proper choice of type 
and circuit for connecting the diode~ and digitally-controlled impedances, allows one to simulate the following 
particular cases of function (1) which are widely used in the design of analog computers, 


a) A multiplier 
z= xy. (5) 
For the multiplication mode it is necessary to provide the following conditions for the transition from 
general relationship (2): 


Too =Y, In =X, Pe =const, (Xo; )j>1 = 0, 


(4 = 0. (6) 
2j /j>1 


b) A divider 
(7) 


z 
z2=>—-——. 
y 


The conditions for the transition from relationship (2): 


Too = const, %g3 = TZ, Toy =—Y, (%2j)j>1 = 0, (8) 


t=) ait 
"93 /j>1 


Thus, the block schematic of the multiplier must contain one linear digitally-controlled impedance in the 
amplifier’s feedback circuit, while the block schematic of the divider must contain the same impedance but in 


the amplifier’s input circuit. 
c) A multiplier-divider 


t= —-w. (9) 


For the transition from relationship (2): 


=) = @. (10) 
j>1 


Toy =U, To =Y, Toi =X, (Xj )j>1 = 0, eS 


The block schematic contains linear digitally-controlled impedances in the amplifier’s input and feedback 
circuits, 
d) A multiplier for the function 


z= yf (2). (11) 


For the transition from relationship (2): 


'o=Y, n= f(x)’ Ty =, (25 )j>1 = 0, hae hac 


(12) 





The block schematic contains one diode-controlled impedance in the input circuit and one linear digitally- 
controlled impedance in the feedback circuit. 
e) A multiplier for the function 


z= f(z) ¢(y). (13) 








For the transition from relationship (2): 


zx = ir { ~ 
7ta) ’ Too = Oy), (%23)j>1 = 9, (=), .=0 (14) 





yt n= 


The block schematic contains a diode~controlled impedance in the input circuit and a functional digitally- 
controlled impedance in the amplifier’s feedback circuit, The latter can be constructed in accordance with 


the block schematic of Fig, 3, or by combining diode-controlled and linear digitally-controlled impedances in 
the amplifier’s feedback circuit. 


a) If, in the general case, the amplifiers input circuit is formed of a group of diode-controlled im- 
pedances connected in parallel, f, (xk), and the feedback circuit by a group of series-connected functional 
digitally-controlled impedances, ; (y;)s then the block schematic of the function generator under consideration 

can be used to simulate a complicated multiplier- divider 
— relationship on the sums of the function fj, (xy) and oy (Y;). 


This latter circumstance allows one to use the generators 
for simulating functions of many variables by using the method 
of superposition, 





5. If the function of two variables to be simulated 


; Zz = f(z, y), 
A Bs 


does not belong to one of the previously considered particular 
Fig. 4. cases, the following technical paths can be recommended for 
implementing it by means of digital-analog function generators 
based on operational amplifiers with controlled impedances, 














a) One carries out a piece-wise linear approximation of the given function by one of its arguments and 
a picce~wise step~function approximation by the other argument, In correspondence with the method of approxi- 
mation, the function generator must have, at its input, n diode-controlled impedances 


zr . 
i= Tey) (j == 1, 2, poy © (15) 


whose connections in the circuit are made in dependence on the ordinal number, j, of the function, i.e,, on the 
value of the second argument, since 


= V(y). 


The function generator’s block schematic is shown in Fig, 4. 


b) A piece~wise linear approximation of the given function is made with respect to each of its arguments, 
for which the current value of the function is determined by the formula 


—%, 
~ —_ 6 
2m f(x, yj) +1) (2, yin) —/ (2, ys) = 9; (16) 
where f (x, y;) and f (x, Yj 41) are the values of the function z for the ends of the intervals of decomposition by 
argument y and rj = f (x, Yj)» Yj+a and y; are the values of the interval of decomplsition by argument y, and x 
and y are the current values of the arguments. 


Since the function generator which simulates the function z = f (x, y) in correspondence with formula (16) 
must provide the operation of multiplication, it is advantageous to connect a linear digitally-controlled im~ 
pedance for argument y in the amplifier’s feedback circuit, 


In order to simplify the generator circuit, the decomposition should be uniform for y and nonuniform for x, 
so that the multiplication of the function increment Af by the argument increment Ay within the limits of the 
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constant portion of the decomposition, yj 44~ yj, is simply 
implemented technically by the block schematic described 
7 ov in section 4, a of this paper. 








The over-all block schematic of the function 
generator which implements the simulation of functions 
by the method of piecewise approximations, and which 
realizes a variation of relationship (16) 











t= f(x, y;)( —*)+7(@, yin) 92, (17) 


_ SS —-— => 


—— 

















Ay =y— Yi, @ = Yj41— yj = Const, 


is given in Fig. 5. The circuit contains three operational 
amplifiers, two of which (OA~1 and OA~2) serve to 
simulate relationships of the form 


ij (1— AY) and fi, = 








while the third carries out the operation of summation, 


The input voltage y is applied to linear transformer 
LT-1 which transforms the voltage to a digital quantity 
which is applied to pulse counter CO-1. The number of 
orders (places) in counter CO-1 corresponds to the con- 
stant interval of decomposition by argument y, so that 
Fig. 6, counting of pulses occurs only in the interval a, after 

which the counter is cleared, and the overflow pulse, 

fixing the ordinal number j of the y decomposition, is 
applied to counter CO-2 of the segments a and then to the switching of the relay for the choice of the ordinal 
number of the function 














13 =] (2, yj). 


Multiplication of fj by a~ Ay, and of fj ,, by Ay, in the limits of the segment of the y decomposition is 
carried out by applying the inverted and the direct codes of the increment Ay to linear digitally-controlled im- 
pedances, respectively, DCI-1 and DCI-2, 


c) The relationship to be simulated, z = f(x, y), is presented in the form 


z= f(x, y1) + (M2, ya) — f (@, y)I Psy). (18) 
If we choose 
j= (2), /(2, wi) =9, (19) 
then 
= 1 (2, Ya) 95 (y). (20) 


The function generator which realizes (20) must contain, in the amplifier’s input circuit, a group of diode- 
controlled impedances % (y), switched as a function of the magnitude of argument x and, in the feedback circuit, 
a functional digitally-controlled impedance, tgp. = f, 


One of the possible block schematics for this function generator is shown in Fig, 6. 


The functional digitally-controlled impedance is formed here by using a diode function generator, not con- 
taining an amplifier, in the x line, 


1346 


























d) In many cases, the given function z = f (x, y) can, with sufficient accuracy, be represented in the form 
of a sum 


z= xf/(y) + (x,y), (21) 


where 


D (x, y) < Zax: (22) 


In this case, the first term is simulated by the generator described in section 4, a and the second term, for 
example, by means of the block schematic of Fig. 6 which, by virtue of inequality (22), contains an insignificant 
number of elements, 


6, For constructing a complicated integrating-differentiating two-input function generator, it is necessary 
to have in the amplifier circuits, in addition to controlled ohmic impedances, constant or controlled capacitive 
impedances, 


For example, to simulate a function of the form 


, a 
= + Kw + 


a K,\ = at (23) 


z2= Ky x 4. K, 


the following transition from relationship (2) is carried out: 


1 
In =T7, Ta =Y, 2 = — = const, "9 = Aw + =, 
2 


(24) 





(72; )j>1 = 9, ( ) = 0. 


Proj /j>1 


The block schematic of the function generator must contain linear digitally-controlled impedance fq in 
the input circuit, shunted by capacitance Cy, and, in the feedback circuit, a series-connected linear digitally- 
controlled impedance Aw plus capacitance Cy. 


For the construction of an integrating-differentiating function generator of higher order, it is necessary to 
use the block schematic of an amplifier with differential inputs; this same circuit is optimal for simulating 
functions of two variables which reduce to fractional-nonlinear relationships of the form 


sus fi (2) fe (y) 
~~ They) + fs (th @) + fay) * 





(25) 
As follows from (3), relationship (25) occurs for (xgj)j > 4 = 9. 


SUMMARY 


The method considered for constructing electronic function generators with two inputs uses, as a rule, one 
operational amplifier and standard hardware of contemporary digital machines: linear coding transformers, 
counters, logical gates, which can be implemented with semiconducting and ferrite elements, 


The accuracy and stability of these generators are determined basically by the accuracy and stability of 
two components, namely, the operational amplifier and the coding transformer, 


The method possesses sufficient flexibility and universality from the point of view of the class of functions 
to be simulated and the accuracy of simulation, 


The function generator circuits considered are significantly simplified if one argument {is given in the form 
of a voltage, and the other in the form of a numerical equivalent, 


The method permits the linking of the components of digital and analog electronic computers without 








intermediate coding and decoding transformers, which can lead to a simplification of the construction of hybrid 
machines, Therefore, the function generators might find application in the design of controlling and measuring 
systems, and in instruments of various purposes, 


The method can only be used for simulating nonlinear integral- differential functions if one of the argu- 
ments is given by a numerical equivalent. 
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ON THE DETERMINATION OF CERTAIN ERRORS 
OF ANALOG COMPUTERS 


N. L. Sosenskii 


(Moscow) 


The paper considers the effect of the nonlinearity of the frequency characteristics 
of electronic analog computers’ (EAC) operational elements in the solution of linear 
differential equations with constant coefficients, To determine the errors in solution 
we use a graphicoanalytic method based on the use of the logarithmic frequency 
characteristics, An example is given of the determination of the error in solving a 
second-order equation for the case when the operational amplifier has a frequency 
characteristic of complicated form and possesses a number of parasitic elements, 


One of the important sources of errors in solving differential equations on an EAC is the nonlinearity of the 
frequency characteristics of the EAC elements, This error shows up most strikingly in solving differential equations 
with a high degree of conservatism (Cf,, for example, [1)). 


Works [1-3] are devoted to the error engendered by the frequency characteristics of the EAC elements in 
solving linear differential equations with constant coefficients, Determination of the error by the methods 
described in [1-3], for operational elements based on multistage operational amplifiers with arbitrary fre- 
quency characteristics and a host of parasitic elements (capacitance between the summing points and ground, 
output impedances, etc) leads to such horrendous computations that, in the cases considered, the use of these 
methods for computing the error turns out to be impossible in practice, 


In this paper, for the determination of the frequency errors of an EAC, we employ the method, widely used 
in automatic control theory, of logarithmic amplitude (LAC) and logarithmic phase (LPC) characteristics, This 
allows one, by using the corresponding graphicoanalytic computations, to estimate the errors in cases where the 
operational amplifier’s frequency characteristic is complicated, and also to take its parasitic elements into 
account. 


Determination of EAC Errors Caused by its Operational Elements’ Frequency Char- 





acteristics 





We consider a linear differential equation with constant coefficients 


n 
y + 3 ay = f (0). (1) 
i=1 
The block schematic for the model corresponding to Eq. (1) is shown in Fig, 1, The use in this circuit of 
operational elements with nonideal frequency characteristics leads to the result that the model solves, not the 
given Eq, (1), but an equation of higher order, and with other coefficients, We now estimate the difference of 
the natural motions corresponding to this equation from the natural motions corresponding to the original Eq, (1), 
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Fig. 1, Block schematic of a model for solving an n'th order linear differential 
equation with constant coefficients, 


Initially, let the circuit of Fig. 1 consist of ideal elements with transfer functions 
{ a 
Ayu = Hy= — —, Hgi = (—1)*Ha\*"', (2) 


where Hyyj is the transfer function from the k'th input of the i*th integrator to its output, Hg; is the transfer 
function of the ith scale amplifier, p is the complex frequency and r is the integrators’ time constant. 


A plus sign in a transfer function Hg; indicates that the corresponding scale amplifier must be replaced by 
a resistance, We shall consider the circuit of Fig. 1 as a system with a multichannel feedback at whose input 
(at point 0) is applied the stimulus f(t). The system's output is the point at which the solution of the equation 
is obtained (point n), We shall denote the transfer functions of the portions of the circuit of Fig, 1 from point 
0 to points 1, 2,..., nm by Yy, Y2..... Y,, and we shall assume that all the poles of these functions lie in the left 
half of the complex p plane, 


Now, let the circuit of Fig. 1 consist of nonideal elements (Fig. 2). Then, the elements’ transfer functions 
are defined by expressions (III)-(XIII)(cf. Appendix I), The logarithmic frequency characteristics of the transfer 
functions of the nonideal eleme nts can be easily constructed from the transfer function, A (p), of the operational 
amplifier (which is considered as being given) by means of nomograms for closed-loop systems [4], We denote 
the transfer functions of the portion of the circuit with nonideal elements from point 0 to points 1, 2,..., n by 
Yu eds Yn and we assume that all the poles of these functions lie in the left half-plane, and also that the 
LAC and the LPC of the functions Yj, for the frequencies w < w*, differ little from the LAC and LPC of the 
functions Yj. Here, w* is the boundary frequency which satisfies the conditions, w* > 1/T;, w*> w;, and 
1/Ty and — nywy * ju; {T= 74; are the poles of function Y;, These assumptions are equivalent to the re- 
quirement that the error in solving the equation on the EAC be small, 


The difference between Y,, (p) and Y,, (p) may amount to a shift of the already existing poles, or to the 
appearance of new poles or zeroes, We denote by p; the poles of function Y, (p), amd by pj; and py, respectively, 
the shifted former poles and the new poles of function Y,, (p)(i=1, 2,..,n; k=n+1,n+2,..,n +m) In the 
case of Eq, (1), to each of the poles p; there corresponds, in the natural motions, a term of the form DePit + %i, 
For the natural motions of the model, there corresponds to each of the poles p; and p, a term of the form 
DyePit + A or DyePkt +k, With this, the quantities Dj and g; are uniquely defined by the initial conditions 
and by the poles pj, and the quantities Dj, $j, Dy, and $,, by the initial conditions, by the poles p, and p, and 
by the new functions Y,, In the present work we do not consider the errors arising from the discrepancies between 
D, and Dj, or between #; and g;. Thus, to determine the error in solution, it suffices to know the values of the 
new poles and the magnitude of the displacement of the old poles, This information can be obtained by a com- 
parison of the logarithmic characteristics corresponding to the function Y; and Y; which are constructed by 
successively applying the nomograms for closed systems to the circuit of Fig. 1. For this, in correspondence with 
the assumptions made above, 


| (o)|=|H (@)|(1+8()], arg H (w) = arg H (w) + @ (0), (3) 
LY: (o)! =|¥i(o)| [1+ e:(o)},  arg¥;(@) = arg Yi (@) + 1: (), (4) 
where 5, a, €4, 7; are small deviations from the ideal characteristics, 








to determine them from the corresponding curves, 
Relationships are given in Appendix II which per- 
mit the determination of the quantities 6(w) and 
a (w) from the frequency characteristics of the 
operational amplifier for the frequency range in 
which | G (jw)/H (jw)| > 1. 


The functions G (jw) and H (jw) are de- 
fined by expressions (IV), (VII), (X) and (XIII) of 
ar Appendix I, Also found there are the expressions 
for determining €4(w) and y, (w) from the 
vf frequency characteristics of the functions 
= + Y¥j-3 (jw) and Y4-3 (jw) for the frequency range 


Fig. 2, Equivalent circuit of the EAC's operational ele- este: « | ¥j-1 Hyg; » 1, 

ments, We now consider how, from a comparison 
of the frequency characteristics Y, and Y,,, we 
can determine the new poles, and also the trans- 
lations of the old poles, The appearance of new 
poles is expressed by the fact that, at frequencies 
greater than boundary frequency w*, the LAC of 
function » will have a slope larger than 

20 n db/decade and may also have a maxi- 
mum (the LAC of function Y,, at fre- 
quencies greater than w®* has a slope of 20n 

db/ decade and has no maxima), It is easily 
shown that all the new seal, —1/T,, and complex, 
— Nk Wy, + jw, {T= Te poles of function & 
are characterized by the relationships 


ae It may turn out that certain of the quantities 
: * 5, aor €, and y; are so small that it is difficult 
fe) 

l 

k 



































Fig. 3, Block schematic of the model for solving 
a second-order equation. 


i = i 
a-~>Q, o =e == 
i > K kms Nk 2 es ’ 





(5) 
where Q is the angular frequency at which the slope of the LAC changes, W,, is a frequency at which a new 
maximum of the LAC is observed and Yum is the value of this maximum, 


It.is also easily shown that a shift of a real pole ~ 1/ T, leads, in the LPC, to a variations with maximum 
value, for w » 1/ Tj, equal to 


aT, 


d om = TT,° (6) 


This same shift, at frequencies w > 1/T;, leads to a variation in the amplitude charactéristic: 





The shift of a pair of complex conjugate poles,~; wi # ju; {I- i leads - a variation in the magnitude 
of Yim and in the position, Wims Of the corresponding maximum in the LAC, For ni « 1, these variations are 
defined by the formulae 





a | Yim! dn, 
TYiml ae n ’ (8) 
7” —= =. (9) 
®im 0; 
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Fig. 4, LAC and LPC of the model consisting of ideal elements. 


Determination of the Error for a Second-Order Equation with Low Damping 





As an example, we now consider the error in the EAC solution of a second-order linear differential equation 
with low damping and constant coefficients, 


y” + 2ywoy’ + wy = 0 (10) 


in the case when the operational amplifiers are three-stage dc amplifiers, The block schematic for the model 
corresponding to Eq. (10) is given in Fig. 3 (RC = 1/w, R'C = 1/2n wo). 


One easily convinces oneself that 








 s 1 2nwo 
{Sap 2n p+ 2nw’ (23) 
» Y,o ws 
ee. PE 12 
° P+ Yia p® + 2nwop + oF : a" 
i 
| Yom | cat (13) 


Figure 4 gives the frequency characteristics corresponding to the block schematic with ideal elements: 
1) | |, 2) arg ®, 3) | ¥,{ , 4) arg Y;, 5)|¥,| and 6) arg Y,. The LAC and LPC of the function Y, (jw) were 
constructed by means of the nomograms for closed systems from the LAC and LPC of the function @ (jw) = 
= wy Yy (jw)/ jw. 


Let the LAC and LPC of the operational amplifiers used in the block schematic of Fig. 3 be given by graphs 
1 and 2 of Fig. 5, These are standard curves for a three-stage dc amplifier with stabilizing circuits, If the 
parasitic parameters are given then, by using the nomograms for closed systems and the formulae of Appendix I, 
we can successively construct the LAC and LPC for the transfer functions of the scale amplifier and the integrators, 
and also for the functions Y; and Y,. Figure 5 shows the frequencies of the scale amplifier: 
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Figure 6 gives the frequency characteristics of the integrator: 
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Figure 7 gives the frequency characteristics 
corresponding to the block schematic with nonideal 
elements: 1) | Y;HyH | » 2) arg Y;HyHs, 3) \¥, | , 
4) arg Y>. 
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All curves constructed are for the following 
values of the basic parameters of the circuit of 
Fig. 3 and of the parasitic elements of the inte- 
grators and the scale amplifier (the parasitic ele- 
ments have been taken as being identical which, 
obviously, does not destroy the generality of the 
discussion), 


1, The angular frequency of the natural 
motions corresponding to Eq. (10) wy = 2° 10*/sec, 
the damping decrement is n = 36-1075, 


am -, 2, The parameters of the elements of the 


60 circuit in Fig. 3; R = 10° ohm, R’ = 1,4+10" ohm, 
C = 5°10" farad, 


3, The leakage impedance on the integrator 
condensers is Rg = 10° ohm, the output impedance 
of the operational amplifiers is r = 10° ohm, the 
37 load impedance for all amplifiers is Ry, = R = 10° 
ohm, the leakage impedance with the summing 
points and ground is R, = 10° ohm, 





Fig. 5, LAC and LPC of the scale amplifier with 
transmission factor k = 1, constructed with an opera~ 


tional amplifier with transfer function A (jw). 4, ‘TES pasuitic capenneee cam berween 


the input and the summing point, C, = 2 nanofarad 

between the output and the summing point of the 
scale amplifier, Cy = 2 nanofarad, between the summing point and ground, C, = 5 nanofarad; the load capacitance 
on the operational element is C; = 20 nanofarad, 


In correspondence with formulae (XXI)-(XXIV), the small deviations from the ideal characteristics [ cf, 
(3) and (4)) will equal 














\H H 
§5(@) = =< COS$s, A 5(@) = | t sin p, + wt (C, — Cp), (14) 
$,(@) = — ft cos ot = bs i RC ‘ (15) 
ye) Se Py % (o) = [Z| sing; + whl, + Sage 

&,(@)=8;(%), %1(@) = z(), (16) 


where g¢ (w) = arg Go/Hg and 9 (w) = arg G,/Hy. 


For a system with very low damping, determination of the maximum of the amplitude characteristic by 
means of the nomograms for closed systems is difficult, since the nomograms are not usually constructed for very 
small values of the phase characteristic of the open-loop system's transfer function, Therefore, to determine the 
magnitude of[¥,,,},we use the obvious expression 





a: Y:h%& 
y¥,--i 714. (17) 
5 Yi Hy H,— 





































Fig. 6, LAC and LPC of the integrator with time constant rt = 1/w, 
constructed with an operational amplifier with transfer function A (jw). 
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Fig. 7. LAC and LPC of the model consisting of nonideal elements, 
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By using the notation of (3) and (4), expanding the individual terms of expressions (17) in series and ignor- 
ing terms of the order of |H/G|? and n? (all the quantities a, y, 53, €4 have order |H/G| ), we obtain 


ot 1—B> 


Ta Se 





(18) 


“ott | 





Since, in the neighborhood of w = wen, the quantity |G/H| >> 1, Le., all the quantities a, y¥ 54 €;* 1, 


it is easy to show that the maximum value of Y; occurs for 





and equals 





7 4 
| ¥am| ds ts (am) + a (@ym) + 71 (®,m) + 27° (20) 


In the case under consideration we obtain, from the corresponding characteristics, 


cos (— 60°) 


- in (— 60° - 
85 (@2m) = Te —4-10°, & < (@2m) ae aL — 6.9-10°5, 
cos (— 20° > 
6 1 (2m) og: Mca ee ) se 5.3-10 s 


in (— 20° 65 > 
Oy (Wom) = Ta + 410% + 2 = 31-107, 


& (@em) = 5.3 ° 105, ‘1 (@2m) = 3.1-10°5, 


By using (19) and (20) we find that Gm » 2°10*/sec and |¥,,,| = 84 db. 


From a comparison of the characteristics of Fig, 7 with the characteristics of Fig. 4, we can draw the 
following conclusions, 


1, Transfer function ¥, has a pair of complex conjugate poles which are translated from their position in 
Y, leading to a variation of the damping decrement and of the frequency of the natural - yomag In correspond- 
ence with (13), the damping decrement, instead of being n = 3.6° 10 “5, equals @ = 3.2 10°5, mh correspondence 
with (19), the frequency equals Wp » 2° 10° sec, 


2, Transfer function ¥, has no new poles, as compared to Y3. 


3, By the appropriate choice of the operational amplifiers’ frequency characteristics one can compensate 
for the EAC's frequency errors with an accuracy determined by the constancy of the operational amplifier*s gain 
at zero frequency, The same result can be attained by using the well-known methods for correcting the phase 
characteristics, or by choosing a sufficiently large operational amplifier gain at the operating frequency. 


An error analysis for the case of very small values of operational amplifier output impedance (< 1 ohm), 
carried out by the method described above, showed that the transfer function ¥, in this case, in addition to a 
shifted pair of complex poles, had new poles, defined by expressions (5). The additional terms in the natural 
motions which correspond to the new poles have the form e~7kWk' cos wt, and are damped by a factor of e 
during time t = 1/n), w, = 2Y Y em/ Wym» Where W,,, is the frequency approximately equal to the operational 
amplifier's cut-off frequency (i,e., the frequency for which | A (jw )| = 1) and Y;,,,, is a quantity defined by the 



















































stability margin of the operational amplifiers, Thus, a larger velocity of damping of the additional terms in the 
natural motions can be provided by a sufficiently high cut-off frequency and by a larger stability margin of the 
operational amplifiers, 








SUMMARY 


A gtaphico~analytic method which is based on the use of logarithmic frequency characteristics allows one 
to estimate, from given graphs of the operational amplifiers’ LAC and LPC, the error, due to the frequency pro- 
perties of the model's elements, in solving linear differential equations with constant coefficients, With this, 
one can determine graphically the fraction of the total error attributable to each source of error in the solution, 
and can also determine the possible measures to take for decreasing and compensating the errors, 


The author wishes to express his gratitude to A. A. Kuz'min and A. A, Vasil’ev for their discussion of a 
number of questions which arose during preparation of this paper, 


APPENDIX I 


Expressions for the Transfer Functions of the Integrator and the Summing Amplifier 





with Account Taken of the Parasitic Elements 





Figure 2 gave the equivalent circuit for the operational element with account taken of the parasitic 
capacities, amplifier output impedance and condenser leakage impedance, By using the method of nodal voltages 
for the amplifier adding points aad the output points, we get 


Rot + PRC, 
> R,t+pRC. * 
U=— , (I) 
{ Fer, + PC) €pn + Com) 


i A(p) + x + prCo| (1 + pRoC,) 
0 








where 
[+1 ‘ 1 ‘ [+1 
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The transfer function. from the summer's kth input to its output can be written in the form 
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The transfer function of the integrator’s k'th input to its output is written in the form 
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Derivation of the Expressions for the Deviations of the Frequency Characteristics | 
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of the Nonideal Elements from those of the Ideal Elements 





We consider the transfer function of some system with a single-channel feedback 





H (p) =H (p)h(p) es 
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We introduce the notation 
G (jo) =|G\e, H (jo)=| Hle*, Hijo) =|Hle, hijo) =jhje*. 


By substituting in (XIV) we get 
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(XIV) 
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If we expand expressions (XVI) and (XVII) in powers of [H/G], limiting ourselves to the first two terms 
of each series, we obtain 
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F=o+8—|% sin () — @) + Br, (XIX) 





where B, and B, are the remaining terms of the series where, for | G/H] = 10, 
H H \2 
By <2.5 Ee Bi<t7| 2. (XX) 


We neglect B, and B, and employ these expressions in the transfer functions of the EAC's operational ele- 
ments, (III) and (IX), Then, by taking (VIII) into account and by ignoring terms of the order of (wr, and 
( wre) we obtain for the summer and for the scale amplifier 
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For the integrator, by taking (XIII) into account and by neglecting terms of the order of (wry) and 
1/(wrs) we get 
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CHOOSING FLEXIBLE FEEDBACK PARAMETERS FOR 
MAGNETIC AMPLIFIERS BY MEANS OF 
INTEGRAL ESTIMATES 


L. V. Safris 


(Rostoy-on-the-Don) 


By means of the method of integral quadratic estimates, we elucidate the 
effect of the parameters of a flexible (velocity) feedback on the duration of the 
transient response in a choke~coupled magnetic amplifier. 


INTRODUCTION 


As is well known, to accelerate the transient responses in magnetic amplifiers, use is made of flexible 
feedback, implemented as a rule by means of an additional winding W, which is connected in parallel with the 
load via condenser C and impedance R, (Fig. 1). 


The recommendations to be found in the literature [1, 2] on the choice of parameters for the accelerating 
contour lead, principally, to the determination, in each individual case, of the value of C for which there would 
occur a critical aperiodic of damping oscillatory transient response, With such a posing of the problem there 
remains hidden the effectiveness of the flexible feedback, in the sense of its effect on the speed with which the 
transient response occurs, Elucidating this effectiveness by constructing curves of the increase in load current 
i, is a very arduous operation. 


We mention that, in some cases, the attainment of a critical aperiodic mode may have little effect as 
far as accelerating the transient response is concerned, 


The method, described in this paper, of choosing the parameters of the feedback path by means of an 
integral quadratic estimate of the form 


J = \(ip— I,)?dt (1) 


ou 98 


(where I; is the steady-state value of load current) possesses the advantage that it allows one simultaneously to 
determine the desirable acceleration of the transient response, As is well known [3], integral estimates are widely 
used for choosing the parameters of automatic control systems, 


We shall consider as optimal those amplifier parameters for which the integral in (1) is a minimum, 


For the basis of our investigation we take the method of linearizing the magnetization characteristic, 
described in [4]. From this same work we also carry over the following basic notation and relationships: W, is 
the number of turns in the control winding; W,, is the number of tums in the ac winding; Wg is the number of 
tums in the current feedback winding; 8 = W,/We; Q is the cross section of the two cores; lis the mean length 
of the magnetic circuit; Rg is the impedance of the control circuit; Ry; is the load impedance; r, is the forward 
impedance of a rectifier bridge arm; R=R;, +1; 4g (Ig) are control currents; 4; (1,) are load currents; 1, is 
the current in the accelerating feedback circuit; ug (U,) are the voltages applied to the input of the controlling 








winding; ey (E.) are the “average” values of emf in the 
ac winding; up (Ug) are the average values of voltage 
on the rectifier bridge input terminals; By is the constant 
(monotonic) component of magnetic induction; Hg is 

the constant component of the magnetic field strength; 
Hie Osffc and yp are coefficients which characterize the 
core properties and which are defined by the magneti- 
zation characteristics; S =o W,/ W,, is the amplifier 
characteristic’s steepness; Rj = yp WiQu/l is the 
amplifier*s internal impedance; kj = R;S/(Rj + R) is the 
current gain; kjg = k;/(1— 6k;) is the current gain with 
inflexible feedback present; kp is the power gain with 
no feedback present; Lg, is the amplifiers input in- 
ductance with an inflexible feedback present; A is a 
symbol which indicates that one has chosen the incre- 
ment of the corresponding quantity within the limits of 
the characteristic’s linear portion. 






























1. Equivalent Circuit and Basic Assumptions 





To compute the integral quadratic estimate, we 
make use of the equivalent amplifier circuit shown in 
Fig. 2, The values of the equivalent amplifier circuit 
parameters were stated above in the list of nomenclature, 
The equivalent circuit is obtained on the basis of the con- 
clusions presented in work [4], in which the following 
basic assumptions are made; 





























a) the values of the magnetic induction and the magnetic field strength do not go beyond the limits of 
the rectilinear region of the characteristic for simultaneous magnetization by constant and variable fields; 


b) the load is active; 
c) the impedance of the control circuit is sufficiently small: 


B* Ro ‘\. 
FR S1): 
d) the rectifiers are replaced by forward constant impedances, and the back impedances are taken to be 
infinitely large; 
e) the transient response occurs “slowly* in comparison with the supply voltage period; 
f) with magnetization of the cores, the increments of voltage at the rectifiers* input terminals and of the 
emf in the ac winding are equal in magnitude and opposite in sign: 
AE. = — AUg, (2) 


i.e., the load curve is approximated by straight line segments, 


By keeping in mind the active load and the low impedance of the ac winding, we can agree that the first 
five assumptions carry over to the case under consideration here, As for the last assumption, it requires special 
verification since, in contradistinction to the cases treated in [4], the load is connected in parallel with a contour 
containing a capacitance, 


Experimental investigations (cf. Appendix 1) allow the following points to be established, 


In the majority of cases, for the steady-state mode, Eq, (2) remains valid with an adequate degree of 
accuracy even when there is capacitance present, As examples, one can take curves 1(C = 0) and 2(C = 4,7 
microfarad) on Fig. 3(cf. Appendix 2), 






































Upy . n With certain values of capacitance, and with low for- 
4 ~ ‘e 8 ward rectifier impedance, the character of the function E.= f, (Ug) 
2 s\ SJ changes somewhat. This case fs illustrated by curve 3 on Fig, 3, 

et vd Here, on individual segments, AU, is significantly larger than the 
0 TRG * corresponding AE,,, The reasons for such a character of the 
8 \ xs At . curves are considered in [5, 6). 

‘ 
‘ | bE \ x ~ In what follows, in those cases when, in a definite range 
» WN\ of variation of Ug, equality (2) is significantly violated, we shall 

4 a Ys replace it by the equality 





























X 
‘; i AE, = — aAUg, (2a) 


Ce Fe where a will be determined experimentally. 
Fig. 3. 








The assumption of Eq, (2) is the basis for considering the 
equivalent circuit of Fig, 2 as being valid in the steady state, 


We shall assume that expression (2) can be extended to the transient response, and we use the circuit of 
Fig. 2 for investigating the steady-state processes while remembering that the final results obtained with such 
assumptions have agreed well with experimental data, 


As was established in [4], the increment of the constant component of the magnetic induction is defined 
by the expression 


y Ae, 
AB, = » 4H, — 7Wwo° (3) 





For the equivalent circuit of Fig. 2, with Eq, (2) being taken into account, Eq, (3) takes the form: 


ao ee) 
ea el (Ai, + Ai.) + W,Q—— : 


ABy = + BaAig+— wa . (4) 








To determine the load current, Ai;, we use the following equations, which follow from the equivalent 
circuit: 








RL AiL+ 2r, (Aiy, + Ai + w,Q (4 Bo) 
S (Aig + BA, + 2, 4%) = Ai, + Aig+ R, , (5) 
WQ2 ne + =| dgdt + Rdig= Ri diy (6) 
W.Q— : ae) + Ry Aig = uo. (7) 


In (4) and (5), B, =(W, +W)/ Wy = W,/ Wo. 


We note in passing that, by discarding the right-most term in Eqs, (4) and (5) and by replacing py. by pj, = 
= AB, /AH,, we obtain equations analogous to those which were taken in [1, 2] as the basis of the assumptions 
regarding the constancy of the input inductance and the validity of equating the variable and constant components 
of the ampere-tums in the transient response, 


If we apply the Laplace transform to Eqs, (4)-(7) with zero initial conditions and a stepwise varying 
voltage up, and then solve these equations jointly, we find that 
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by —by *t) )p+(t— =] 


Ait (p) = kip ARS BE +! 


agp" + ap+i 














b= Ro — E+ itn): 


by = FE fe [m Co — 6) — “| + sae [m(Pa—P) — Rem (®) 





a, = 6 [1 —BRo( a — eye) | + t# (1+ RR) + 
+ C [Ro — Gekip — ma)Ry, (10) 


ay = C {v0 heneline [ee—p(1—54)|+2@,—p)}+ 
+ Ralt + BRo(TeaRR, Rkg)|+ nRy (1+ oe Ry) t 
+t, [Ro Pa B) (Gakio — ms) 6 — Bal (1 — Bk) + Ra(1+ gap) + 
+R, (™+7-5°R)}}- sa 


For brevity in writing Eqs, (9), (10) and (11), we introduced the following notation: 














a ae thd = t = Wee 
a) a ~ ee? "Se: hoe? "8 oe? 
tT u,R Tt 
a io lta as gant 
‘o-i—m, “sak “IKE, 


The quantities rz, and Tyg are the components of the control circuit's time constant (tag + T,,, = 
= Leg / Rg) which, respectively, do not, and do, depend on the load impedance, 


To simplify the analysis, we assume that, with small values of Rg, the gain does not depend on the 
magnitude of capacitance C, Actually, as was noted in [2] the presence of the flexible feedback engenders a 
certain increase in k, due to the even harmonics, 


For ideal cores and diodes, when r, = 0, m =1,n =0 and R; » R, the coefficients of Eq, (8) take the form: 





j BR BR 
by = RC — igh b, = — keh; RCrp, (9a) 
a, = typ (1— a) — Cl Cake — mo) A — Rd, (10a) 
BaRo BR 
ay = Crp | E— Pabo — mp) + Ra(1— *E)). (ila) 


2. The Expression for J and its Analysis 





To compute the integral J, corresponding to the amplifier's stable mode, we use the method presented in 
[7]. 


We shall be interested only in the transitional component of the load current, so that the term b,/ a, in 
(8) is dropped, As a result, we obtain the expression 


. be +b) 
Ai (P) = bo Rp ap + op Ft 


where be = by-byay/ aq and by = 1-by/ a», 





(12) 


































By retaining (with added primes) the notation of [7], we 


ant A’ = a,a,—1, 


A; _ A, Ao - aid, + a2, B, - bo, B, — b?. 


Consequently,] 





a (k z) a® (a; — bi)* + (a2 — bs) (13) 


iB Ry 24) 


We note that if C = 0 then b, = 0 and ag = 0, As a rule, 
BRo/ (Rj +R)ky « 1 and B*Ro/ (Rj +R) « 1, We can therefore 
assume that ay~by & a ¥ Tyg + Tag. Then, 
Uo\* Tp + 
Jc=o = (Ain Ri.) ak . 
Fig. 4, 
For subsequent analysis, it is convenient to consider the 
dimensionless quantity 


ye J __ 22 (a1 — bi)® + (a2 — 63)? 
Jon a\a2 (Typ T Tap) 








; (14) 


The coefficient @ allows one to judge how the duration of the transient response after introduction of the 
accelerating feedback changes from its value when no such feedback is present, 


Expressions (8)-(11) and (14) and all their consequences remain valid when there is no rigid feedback (of 
current), if it is borne in mind that 6 = 0 and, consequently, Tyg + Tage kig and mg must be replaced by, 
respectively, Ty, Ta» kj and m, 


If the numerical values of the amplifier parameters are substituted in (14), one can then determine the 
quantities C, B, and Rg which correspond to the least value of @ and for which, consequently, one can expect 
the best flow of the transient response, 


We begin our analysis with the particular, and simplest, case of an ideal amplifier, Since, in any amplifier, 
the power is BRo /Kig Ry = B kig Re Mig Ry « 1, it is then easily established from (9a) and (11a) that ag-b, = a, 
Thus, for the case under discussion, 


9 — (a— b;)* + as 
7 aT yg : 





(14a) 


By substituting the values of by, a, and a, from (9a), (10a) and (11a), we can write expression (14a) in the 
following way: 





R 
[tp — ¢ (By *ig — ™g) A, |? + Crp) it (Bakig — ™p) + Ral 
0 iis : 


TB (typ — © [(B hig — mg) a R,}} . (15) 





We first explain under what conditions the flexible feedback can provide an acceleration of the transient 
response, i.e,, when the condition @ < 1 will hold, If we take the difference between the numerator and the 
denominator of expression (15), the condition 6 < 1 is transformed to the inequality 


8.R 
CR, (Bgkis — mp) [Batis — mg) CR L— vp (1 ae ie) <0. (16) 














It is easily verified that, with the values 


Bak, —m 
Batis — mp — 1— Bk, <0 (17) 





BaRy/kig RL « 1 and inequality (16) does not hold, 
Since m = 1 and ky » Wy / W,,, condition (17) is 
equivalent to the condition W,/W,,=1, In other 
words, the effect of the flexible feedback will not 
be realized if the current branching off through the 
capacitor after mutliplication by a factor of 8 ,k; 

Fig. 5, induces in the load a component which is less than 

4i,. With such a choice of parameters, the in- 

clusion of the flexible feedback leads to a more of less delicate retardation of the transient response, 





It isnecessary to mention that the conditions for choosing parameters which were given in [2] are inadequate, 


since the appearance of a periodic component in the curve for current iy, does not always lead to an acceleration 
of the transient response, 


As an example, Fig. 4 shows the oscillograms of the load current i, obtained in the mode when 8 ,ky = m 
(cf. Appendix 3), As is obvious, the appearance of a periodic component (Fig, 4b) has no practical effect on 
the duration of the transient response, which remained the same as when there was no flexible feedback (Fig. 4a). 


Numerical verification by means of expression (15) for a number of amplifiers, and also the experimentally- 
founded statements given in [2] allow one to assume that, all other things being equal, an increase in R, leads 
to a retardation of the transient response, Therefore, it is always advantageous to have a minimal value of R, 
(the resistance of winding W;). 


Thus, the problem of choosing the parameters of the flexible feedback reduces to determining the optimal, 
mutually corresponding, values of C and 8,. 


In the general case, the answer can be obtained by solving the system of equations 00/@C = 0 and 0@/ 08, =0, 
In view of the complexity of the expressions for C and 8 we use the following route, 


From the equation 86/8C = 0 we find the optimal value of capacitance as a function of 8: 
et Typ ( ss ry (“ase — mz) Ry _ A 
Cont = (aki, —™,) RL— Ry 1 3 Bakig—™p) Ry + 


R 
‘alaig Cate mo) + Re] Cake ~ R1) 
Ty \(Bakig — ™%,) Ry }* ‘ 




















+ (18) 


Due to the smallness of R, in comparison with (6 ,kig — Mg) Ry (except in the neighborhood of points where 
Bakig = mg)sthe radicand can be simplified somewhat, and we take 


a Typ wh Balto 
Con = Bfg=mii = 7, (1 Vy rn) (19) 


By substituting (19) in (15), and taking certain simplifications into account, we find the relationship for 
C#Cy 
pt 








1) <2 Bho BARo R, 
2) = keh igh, * @akig—™) Fy 





(20) 


Expression (20) allows one to compute the very best value of 8,, corresponding to the minimal value of @. 


As an example of such a construction, Fig. 5 shows, on curve 1(cf, Appendix 4), the representation of 






































ese” 
Fig. 6. 


the function 6 = f; (Bgky/m) for C = Copt The same figure 
shows, by curve 2, 6 = f (8 ,ky/m) for C = Copt for the same 
amplifier with no rigid feedback (6 = 0), 


A comparison of these curves established that the flexible 
feedback is significantly more effective when a rigid feedback 
is also present, 


As is obvious from (20), for R, = 0, and all other con- 
ditions being equal, the least duration of the transient response 
will be attained for 


oy / ta PR 1— Bh 
2 : 


kp 





Here we have taken into account that the least value of 
Bq (for Co,, > 0) equals 1/k;, It follows from this expression 
that an increase in the rigid feedback coefficient 6 entails a 
decrease in 6, i.e., amplifies the accelerating action of the 
flexible feedback, 


On the other hand, curves 1 and 2 of Fig, 5 attest to the 
fact that an extreme increase of 6,,even with the proper 
choice of capacitance, leads to an increase in 6. 


We recall that we have been considering the case of an 
ideal core and anideal diode, 


We now explain how the value of the coefficient 6 varies 
if we take the actual amplifier parameters into account, 


To this end, we determine the expression for 9 by means 


of (14) by substituting the numerical values of the amplifier para- © 


meters in (9), (10) and (11), Then, from the equation 06 /8C =0 
we find Cop as a function of Bg. Finally, by assigning dif- 
ferent values to Bg, we compute 6 = fy (Baky/m) for C =Cop, = 
= & (Bakj/ m). 


Curves 3 and 4 of Fig, 5 correspond to these last two 
functions (cf, Appendix 4), 


As is obvious, curve 3 has a minimum which is significantly 
less sharply expressed than the one for curve 1, it being dis- 
placed to the region where §,kj/m = 2to 3, Within the limits 
of these values of 6,k;/m there must occur a greatest velocity 
of transient response which will be less than for the ideal 
amplifier. 


Curve 3 agrees well with experimental data, Figure 6 


shows oscillograms a, b, c, d, e and f which illustrate the course of the transient response for different values of 
8 ,kj/m when C =Con (cf. Appendix 3), The values of capacitance were chosen experimentally, The para- 
meters of the amplifier tested (cf, Appendix 1) were close to those for which curves 3 and 4 of Fig, 5 were cal- 
culated (cf, Appendix 4), The oscillograms correspond to points a, b, c, d, e and f on this figure 


Thus, to the degree that the parameters of the amplifier under consideration deviate from those of the ideal 
amplifier, the effectiveness of the flexible feedback will be decreased, 





If one wishes to avoid the difficult computations described above, one can take, as a first approximation, 
Baky = 2 to 3 for high-quality cores, and 8,kj = 3 to 4 for the rest, An approximate value of the corresponding 
Copr can be found by means of the following approximate formula which follows from (18) for B gR/ kigR; « 1: 
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Copt = GA mA, (21) 





In conclusion, we turn our attention to the fact that this 
same expression can be obtained if we start from the following 
considerations, We find a value of C such that the first derivative 
of the frequency characteristic, i, (jQ), in the neighborhood of 
the point N = 0 would equal zero, Here, 2 denotes the fre- 
quency of the controlling signal. 


By using expression (8) we find that diy, (jQ)/d(jQ) = 0 
for 2 = 0 if ay~by = 0, By taking (9a) and (10a) into account, 
we find an expression for C which coincides with (21), 


8. Peculiar Features of Certain Modes of 





Operation 





The discussion given above remains qualitatively valid 
also for a < 1 in formula (2a), 


Quantitatively, naturally there are changes since the 
taking into account of expression (2a) leads to the change of a 
number of the coefficients of Eqs, (4) and (5), which depend on 
a(cf. Appendix 5), 





In addition, the computations are significantly com- 
plicated, since the coefficient a depends on the magnitude of 
capacitance C and on the range within whose limits the varia- 
tions in Ae, , occur, 


Fig. 8. 


This circumstance makes it mandatory, in the choice of parameters for the feedback loop in this case, to 
use expression (14), while keeping in mind the possibility of deviation of the actual duration from the computed 
one, these deviations being, as a rule, on the side of accelerating the transient response, 


In particular, if a< 1,8 = 0,6, =0 andR, = 0, then expression (15) takes the form: 


6 =a(14 =). 


When CR, « at,, 9 # a< 1, This means that, by connecting a capacitance across the load, one can, 
in certain cases, obtain an acceleration of the transient response even when no flexible feedback is present, 


As an illustration, we give the oscillograms of Fig. 7(cf. Appendix 3), In taking off the oscillogram of 
Fig, 7b, the process occurred more rapidly, in spite of the fact that the load was shunted by a capacitance which 
was absent in the circuit corresponding to the oscillogram of Fig, 7a. Analogous conclusions, obtained by a 
different route, are contained in [8]. 


In the case of an inductive load (R;, Ly), the analysis is complicated due to the raising of the orders of 
Eqs, (4), (5) and (6), and also due to the necessity of taking into account the switching in the diodes, the character 
of which, in its turn, depends on the capacitance, Since the expressions obtained are quite horrible, they are 
not included in the present work. 


The use, with an inductive load, of the feedback with parameters chosen under the assumption that L; = 0 
gives satisfactory results, However, it is necessary to bear the following in mind with this, Even with no feed- 
back present (6 = 0 and 8, = 0) and with the proper choice of capacitance, autooscillations of the relaxation 
type may be observed, Figure 8 (cf. Appendix 3) shows the oscillogram of such autooscillations. 












1, By means of an integral quadratic estimate, one can choose the optimal parameters of the feedback 
path and avoid the errors which could arise if the parameters were chosen from the condition that a critical 


SUMMARY 


aperiodic transient response be obtained, 


The increase of speed of the transient response due to the introduction of a flexible feedback will be 
greater when there is a rigid feedback present and for values of 8 ak; = 2 to 4, 


2, When the load is shunted by a capacitance and there is no feedback present, there may be observed, 
in certain cases, an acceleration of the transient response and the generation of autooscillations. 


APPENDIX 1. 


Description of the Amplifier Tested 





The amplifier cores were spun toroids of EN80Me (NII TVCh); Q = 2 x 0,325 cm’, 1 = 0,13 meter, py, 7 62 
millihenries/ meter, o = 0,98, p,, = 0.6 millihenries/ meter, y = 1,38, The amplifier windings have the follow- 
ing data: W,, = 190, Wy = 50 and 200, W, = 185, 6 = 3,7, The amplifier’s internal impedance is Rj = 3500 ohm, 
its characteristic steepness was S = 0.26 and 1.05. The amplifier load was Ry = 184 ohm, L; = 25 henries; the 
diodes were type DG-Ts25, r, = 3 ohm. 
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The curves of Fig, 3 were obtained from the specimen described in Appendix 1, The voltages are ex- 




















pressed in average values and were measure by a rectifying instrument 








Curve for L; = 0, 
on C. uf L 
Fig. 3 
| APPENDIX 3. 
2 3.7 | 10.3 | 4.7 
- _ _ 1.25. The following parameters and modes of operation correspond to 
the oscillogra ms; 
Figure | Baki 
number B B, C.yf |, 41,.ma | Initial magnetization ma 
4 a _ _ _— 80 20 1.07 
b — 4.0 42.0 80 20 1.07 
a 3.7 _ 10 Due to feedback 0 
b 3.7 3.7 285.0 10 Ditto 1.0 
c 3.7 4.3 60.3 10 * «8 1.45 
6 d 3.7 5.7 18.3 10 Hae 1.5 
e 3.7 7.7 9.5 10 on 2.05 
f 3.7 10.3 6.0 10 2.75 
i a i iss 20 10 0 
b — _ 1.25 20 10 0 
8 — te 1.25 _ 10 0 























Note: On Fig, 6b, the magnitude of capacity was chosen arbitrarily; on Fig, 8, Ly = 25 


henries, 


APPENDIX 4. 


The following parameters correspond to the curves of Fig, 5: Ro = 2.4 ohm; Ry = 185 ohm, R; + 2r, = 190 
ohm, m = 0,95, n = 0,075, R, = 4 ohm, 1/(1— 6 kj) = 16, kj, = 4.1, Ty, * 3.7 millisecond, Tag = 6,4 milli- 


second, 








For curves 1 and 2, rg = 0, m=1,n=0, Tag = 0, 
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For a # 1, Eqs, (4) and (5) take the form: 





Ww 
Pe (Aig + BAI Lt B Ai) + 


l 
j ; d (ABs) 
av Ry Ai; + 2rq (Aj; + Aig) + WeQ—Z 


to W,0 
SR; (Aig + BAiL+BAi) = aR, Ai» + (Ry + 2ar,) (Ai, + Aid + aW,Q 


ABy= 





d (ABs) 
a 


LITERATURE CITED 
{1] M, A, Rozenblat, Magnetic Amplifiers [In Russian], page 684 Sovetskoe Radio (1956). 


[2] O, L Aven, “Decreasing magnetic amplifier lags by the introduction of flexible feedback,” [In Russian], 
Automation and Remote Control (USSR ) 18, 2 (195'7).* 


{3} A. A. Krasovskii, “Integral estimates and the choice of parameters for automatic control systems,” 
{In Russian], In the book, Automatic Control Fundamentals [In Russian] (edited by V. V. Solodovnikov). Mashgiz 
(1954), 


(4) L. V. Safris, "Magnetic amplifier design by linearized magnetization characteristics," [In Russian], 
Trudy Rostovskogo Instityta Inzhenerov Zh,-D. Transporta 20, Transzheldorizdat (1956), 


[5) V. L Molotkov, *The operation of a magnetic amplifier on a diode bridge with capacitive filters," 
In Russian), Trudy Leningradskogo Politekhnichesogo Instuta Im, Kalinina 194 (1958), 


[6] N. A. Kaluzhnikov, “On the design of a choke~controlled magnetic amplifier connected to a single- 
phase diode bridge,” [In Russian), Automation and Remote Control (USSR) 19, 3 (1958).* 


{7} A. A. Krasovskii, "On the degree of stability of linear systems,” [In Russian}, Trudy VVIA Im, 
Zhukovskogo (1948), 


(8) V. L Molotkov, “Investigation of the transient response in a magnetic amplifier with capacitive 
filters at the output of a rectifier bridge,” [In Russian], Nauchno-tekhnicheskii Informatsionnyi Byulleten’ LPI 
Im, M, L, Kalinina 10, Radiofizika (1957), 


Received November 3, 1959 


* See English translation . 
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Ways are considered for increasing the effectiveness of information trans- 
mission by telemetry on the basis of taking into account such spectral characteristics 
of the measured parameter as the dependency on frequency of oscillation of the 

amplitude of the oscillations and the measurement errors, 






The paper presents a method of transmission with automatic changeover of the 
coding system, Estimates are made of the increase in speed of transmission and of 
the possible decrease in bandwidth, 








Increasing the effectiveness of transmission allows one to decrease the bandwidth or to decrease the trans- 
mission time, and is attained, for example, on the basis of taking the properties of the communications inte 
account in the design of the coding system, In information theory, Shannon has given a general definition which 
includes the statistical aspects of information transmission, The effectiveness of transmission is defined as the | 
ratio of “the actual speed of transmission to the channel capacity" [1]. For judging practical transmission systems 
it is advantageous to use the concept of relative effectiveness (efficiency) n, where we understand by this the 
ratio of the maximum speed of information production by a given source, Rs, to the maximum speed of trans- 
mission of information analogous in its statistical and spectral properties by means of a coding system to be em- 
ployed, Rc, 1 =Rg/R.. Sometimes the increase in efficiency is conveniently estimated in the form of the gain 
($) in speed of transmission (R): 9p %=(1— (Ro/Rcp)) 100, or Sp =Roo/Rc, where Ro» and Rc are the original 
and the newly obtained values of speed of information transmission by the coding system, However, information- 
theoretical estimates (of transmission rate, for example) are ordinarily averages (over individual transmissions, 
over time) so that the use of them for comparing a number of transmission systems, telemetering systems in 
particular, is not always indicative. The specific feature of telemetry is that the lags connected with optimal 
coding are ordinarily inadmissible and, consequently, it is impossible to increase the average transmission rate 

by this method,* Therefore, the increase in efficiency is more frequently estimated directly in terms of the 
physical indicators of the transmission system, for example, as the gain in bandwidth when one switches from 

the initial bandwidth W, to a newly attained bandwidth W, Le., 








W ’ , Ww 
dw % =(1— 9) 100% or w= 7: 


1. Spectral Characteristics of the Parameter to be Telemetered 





It is possible in many cases to increase the relative efficiency of information transformation by taking 
into account the spectral characteristics of this information in designing the coding system.** In the given case, 
we consider as spectral characteristics the dependence on frequency f of the amplitude of the oscillations A (f), 






*To the author’s knowledge, this property of telemetered messages was first formulated by V, A, Il"inoi, 
* *Here, the expression "coding system” is understood in the broad sense, e,, as the transformer (cipherer, 
modulator) of one form of information-bearing signal to another form convenient for transmission along the 
communication channel, 
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the admissible absolute error A (f) in the measure- 
ments and the number of gradations, N (f) = A (f)/ 

4 (f), mecessary for transmission, These char- 
acteristics can be determined from an analysis, not 
only of the message source, but also of the receiver. 
Indeed, the properties of the receiver which render 
unnecessary the transmission of some elements or other 
of the messages must be taken into account in deter 
mining the characteristics of the messages to be 
transmitted, Thus, the spectral characteristics can 

be given to the input of the coding system, 
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Fig. 1, Spectral characteristics of certain parameters Figure 1 shows the most typical spectral char- 
subject to telemetering. acteristics which are conventionally divided into 

four types, We note that the second type of char- 
acteristic corresponds to a decrease in the amplitude of the oscillations, and the third to a increase in admissible 
error, with increasing speed of variation (frequency of oscillation) ofthe telemetered parameter, The fourth case 
is a combination of the second and third cases, 








Ordinarily, telemetry devices are designed in correspondence with characteristics of the first type. How- 
ever, if the actual characteristics of the telemetered parameter are different from characteristics of the first type 
then, in many cases, it is disadvantageous to design a device for the maximum frequency of signal variation and 
the maximum number of gradations, 


One method of matching channel (or coding system) and signal spectral characteristics consists of divid- 
ing the frequency band occupied by the signal into portions with characteristics close to uniform, and then de- 
signing several coding systems in correspondence with these portions [2, 3]. 


We propose another method of matching channel and signal characteristics for telemetry, this method 
consisting of the construction of an automatic switchover between coding systems, automatically varying the 
number of transmitted gradations, at the expense of varying the speed of operation, as a function of the speed 
of parameter variation,* Automatic switchover is essentially a change of a component of signal volume V ** 
although, in many cases, for example in pulse-code modulation, the space itself, V = 2fTN, is changed, 


2. Increasing Efficiency for Characteristics of the Second Type 





We now consider the interval between pulse transmissions, corresponding to a parameter change with fre- 
quency f (Fig, 2), The parameter values x; and x; +4, separated in time by 1/ 2f, can differ only within the 
limits of the zone + A (f), The actual gain in information at each ordinate can reach the magnitude AI = 
= log [A (f)/49] = log N (f). Since the number of gradations which distinguish one transmission from another 
does not exceed N (f) < No (Fig. 1,11), it is then advantageous to transmit, at each transmission, precisely that 
number N (f) of gradations within whose limits lie the possible parameter variations at the given frequency, In 
telemetry, however, the transmission of just one change is, as a rule, inadmissible, since it is necessary to have 
correction indications to exclude the possible observational errors, This leads to a device of such a structure 
that the transmission of parameter variations is joined with the more infrequent transmission of absolute values, 


The block schematic of a device in which matching of source and channel is implemented by means of 
several, two for example, coding systems can be the following. A slow-acting device with pulse-code modulation 
(PCM) transmits, with a long period, the accurate values.of the parameter, By a separate channel are trans- 
mitted the rapid changes, for example, by the method of A-modulation or by one of its variants, distinguished 
by the fact that, in the absence of parameter variations, no signal is transmitted along this channel, Figure 3 
shows the block schematic of a device which implements both these functions, Without going into details, we 
note that the reyersible counter is very convenient in such code devices, and eases the isolation of deviation 


*The character of such transmission has analogies in certain biological processes, for example, in certain human 
reactions, 

**Such a possibility was stated by A, A. Kharkevich [4], The definition of signal volume used here is some- 
what different from that given by A. A. Kharkevich, which is closer to the established definition of speed of 
transmission, 
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Fig. 3. Block schematic of a multichannel 
telemetry code device transmitting, as an 
additional feature, the rapid variations of 
one of the parameters by an individual 
channel, 1) the null-organ (the output is re- 
presented conventionally by contacts); 2) 
the hinary reversible counter; 3) transforms 
code to voltage; 4) the channel switcher; 5) 


the repeater-index and 6 the inquiry switcher, 


The subscripts in parentheses are the channel 
numbers, 


pulses whose values can differ by units, The necessity, in 
principle, of three channels for transmitting one tele- 
metered quantity sharply decreases the desirability of such 
a method, 


A more promising method of telemetering is the 
construction of automatically changed-over coding systems 
which transmit deviations and absolute parameter values 
over one individual channel, Of the possible variants, we 
mention the following code transmission: for the rapid 
variations of the parameter, there are two~pulse code 
transmissions of one of two combinations which assume the 
values +1 or — 1 in the system of A-modulation,* For 
transmitting a constant component, when the parameter 
varies slowly, each such two-pulse transmission consists 
of one of two other combinations which correspond to 1 
or 0 in the binary code, 


3. Increasing Efficiency for Characteristics 














of the Third Type 





Characteristics of the third type correspond to cases 
when the admissible telemetering error increases with in- 
creases in the frequency of oscillations (i,¢,, the speed of 
variation) of the parameter, and there is one and the same 
maximum value of these oscillations, Le., A (fay = Age 
In telemetering devices which transmit such information 
it is also advantageous to use automatic changeover of 
coding systems, 


Thus, for example, during relatively stable operation 
of a dispatcher-controlled energy system, one may record 
accurate data by telemetric means, process these data and 
then carry out the acts called for by these data, In the 
case of an abrupt variation in the controlled parameters, 
due, for example, to an accident, the dispatcher cannot 
judge the then occurring events with the same original 


accuracy, so that the transmission of telemetry data can be rougher in this case, 


At the same time, it is important to the dispatcher to obtain these data with the least possible delay, 


Automatic changeover telemetry devices can be implemented with any type of modulation, The most 
easily realized are telemetry devices with pulse-width, pulse-time and pulse-code modulation (PWM, PTM, 


PCM). 


In the case of PCM, several methods of implementing discrete automatic tuning exist, One of these con- 
sists of the development of a coder which transmits, for slow parameter variations, code transmissions consisting 
of ny places of a binary code and, for rapid variations, shorter transmissions of ng < ny places of the code, A 


signal of the transition to the rapid, but rough, transmission might be, for example, the absence of a synchronizing 
pulse, The receiver must have a decoding device which discriminates the short and long transmissions and decodes 





*The values of these A-transmissions can also differ from units, which corresponds to characteristics of the fourth 
type, Moreover, the elementary transmission may be, not two-, but multi-pulse (cf, section 3), 


Multi-channel transmission with time separation is implemented by such a method, The specific feature 
here is the necessity of transmitting initially the first elements of the codes for all channels successively, then 


the second elements, etc, 





































AE RZ i them, Another possible variant is analogous to that described 
in section 2 with this difference, that for rapid parameter 
- Sei. at variations, the n, pulses transmit, not the parameter deviations, 
a ? but the full scale, with an error of 1/2"— 2), With this, two 
Pie A or combinations in this transmission are used for the transmission 
. of 1 or 0 in the binary code with n, places (ny > ng) used in 
transmissions for slowly varying parameters, 

























With PWM or PTM, proportional automatic tuning is 
carried out, Figure 4a shows an example of the change in 
duration of the measurement cycle (transmission) T in the case 
of PWM for two velocity changes of the parameter (r is the 
measurement time interval), 











Figure 4b gives the block schematic for an automatic 
changeover transmitting device with PWM (or PTM), The 















nog transducer voltage, which is proportional to the measured para- 
channel meter, is compared with the saw-tooth voltage from generator 
2, Null-organ 1, which reacts to the difference of these volt- 
ages, develops pulses whose width is proportional to the trans- 
Fig. 4, Pulse~width telemetry device ducer voltage, With parameter variations, differentiating 
with automatic change-over (tuning). a) block 3 generates a voltage (proportional to the parameter's 
the cycle variation for abruptly changing first derivative) which increases the frequency of the generator, * 
parameters; b) the sending set: 1) the So as not to distort the measurement cycle, the change in 
null organ; 2) the saw~tooth voltage saw-tooth voltage generator frequency is carried out at the 
generator; 3) the differentiating block; moment when this cycle terminates, 
4) the relay element; 5) the output pulse 
sharper; c) the receiver: 1) the de- To exclude the delays in telemetry transmission which 
modulater of ieasurement time interval are possible with this, it can be provided that a sufficiently 
1; 2) the demodulator of cycle time T; large second derivative of the parameter could interrupt a 
3) the divider for r/T. cycle already begun, after which a new cycle would commence 






with a new rate of increase of the saw-tooth voltage, An in- 
dication that an incomplete cycle had been interrupted can be 
transmitted, for example, by the simplest code transmission, 








For the implementation of these functions there is a relay element, 4, in the scheme, whose executive 
organ is conventionally represented by contacts, The relay element is controlled by the differentiating block 
(the second derivative) and by the generator (at the moment of cycle termination), The invariance of the slope 
of the saw-tooth during the measuring cycles is provided by a memory device in the form of a condenser, 








The receiving device (Fig, 4c) contains two demodulators, 1 and 2, for the separate demodulation of the 
time interval, proportional to parameter r, and the current value ofcycle time T. For obtaining indications, 
it is necessary to divide these quantities by means of special divider 3, for example, a ratiometer, 








When the pulse transmission repetition rate is increased, there is a decrease in the telemetry accuracy of 
the code variant of the telemetry device due to a decrease in the number of transmitted places in the binary 
code, and a decrease in the pulse-width variant due to a decrease in the absolute value of the time deviations. 






To estimate the gain in efficiency of such systems, we consider an example, Let the telemetry device 
with PWM transmit a parameter whose characteristics are given in Fig. 1, DL 






_ Asa preliminary, we note that a nonchangeover coding device, designed for upper frequency f,, and 
number of gradations No, must transmit, per second, 





a bits we _degrees 
Fe 2/ plogs Na — or V 27No pram 






*The differentiating block operates on the modulus of the derivative, i,e., its output signal is independent of 
the sign of the parameter variation. 
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In our example, we take the case of a broken-line function 4 (f), corresponding to a PWM device with 
automatic tuning in which the repetition frequency of pulse transmissions is made proportional to the frequency 
of parameter oscillation, starting with some value, f,, of this frequency. With this, the number of bits per second 


will depend on the frequency, Ro (f) = 2f log N (f), and the number of gradations transmitted per second does 
not change, Le., V (f) = Vo 


In this somewhat idealized example, the rate at which the source generates information equals its speed 
of transmission by the coding system, n =1, The gain in transmission rate is Sp = Ro,/ max Ro = log Nog/ log 
N (ff), since the value of R (f) is not constant in the range 0 < f< fys and max Ro (f) = Ro (fy ) = 2ftylog N (fy). 


We now determine the gain in bandwidth, If, with a nonchangeover device, the bandwidth in the channel 
is proportional to fj; then, in our case, this bandwidth is proportional to frequency f, with easily implemented 
conditions on the invariance on the edge steepness in the process of automatic changeover, The gain in band- 
width can be significant: Sw = Wy/W = fy/ fh. 


SUMMARY 





The advantages of special measures directed towards increasing telemetry efficiency are defined, certainly, 
by a number of factors and, in the first instance, by the simplicity of their technological implementations, To- 
gether with such well-known methods for fuller use of communications channels as improving the quality of 
elements and components (stability of filters, generators, etc), increasing telemetering efficiency on the basis 

of taking the spectral characteristics of the measured parameter into account allows an additional manifold de- 
crease in bandwidth or increase in tele-transmission speed, 


To match the information source with the channel in telemetry, it is advantageous to use an automatically 
changed-over coding system which varies the number of transmitted gradations and the rate of telemetering as 
functions of the speed of parameter variation, 


It is necessary to mention that although the basic types (second and third) of spectral characteristics con- 
sidered contain decreasing functions of numbers of gradations on frequency, i,e., N (f), nonetheless the difference 
in characteristic A (f), and also 4(f), determines the specific features of the transmissions systems for parameters 
of each type, 


Whether some spectral characteristics of parameters or others are typical cannot be definitively decided 
today. However, it seems most probable that, in industrial telemetry, the most widespread will be the para~ 
meters with characteristics of the fourth type where, with increases in frequency, there are permitted, not only 
increases in the measurement errors, but also decreases in the amplitude of parameter oscillations, 
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SWITCHING CIRCUIT OPERATION DURING TRANSITION 
PERIODS. 


V. N. Roginskii 


(Moscow) 


A method is given for disclosing the possible disruptions in switching circuit 
operation during the periods of operation or release of the individual relays in the 
circuit, and recommendations are made for avoiding these disruptions, 


In the course of switching circuit operation, the changes in individual relay states do not occur instant- 
aneously and, during the “transition periods,” it is possible that the contact circuits will be disturbed due to the 
fact that individual contacts of one and the same relay do not react (make or break) simultaneously, Thus, for 
example, if, at operation of a relay, a break contact is initially broken after which a make contact is closed, 
there will then exist a period during which both contacts will be open, If the contacts react in the reverse order, 
there will then be a period of time during which both contacts will be closed, Either one or the other condition 
might lead to a disruption of proper circuit operation since, in the first case, there might be a momentary open 
circuit and, in the second case, a false circuit might appear, This, in its turn, might lead to the incorrect 
operation or release of a relay or other device on which the given circuit operates. 


The nonsimultaneous making and breaking of the individual contacts of one relay might be caused, on the 
one hand, by the design of the contact group, and can be systematic and, on the other hand, by the mechanical 
properties of the individual springs, and may be random, Analogous phenomena are also observed in cases when 
two or more relays in the circuit simultaneously change their state and, due to temporal separations of their 
operations or releasings, their contacts do not react simultaneously, 


Switching theory [1, 2), at the beginning of its development, studied circuits only in their “static™ states, 
and did not consider the transition periods, Recently, a number of works, [3, 4, 5) have appeared which consider 
individual aspects of switching circuit operation in the transition periods, In the present work, we present a 
general analytic method for accounting for switching circuit disruptions during the transition periods, and we 
consider methods for compensating these disruptions, 


1. The Appearance of Switching Circuit Disruptions in the Transition Periods 





The transition processes in switching circuits will be considered as the severing of mechanical connections 
between individual contacts of one relay, i.e,, these contacts will be treated as contacts of different " ficticious" 
relays, analogouslyto what was done byGr. Moisil and Gh, Ioanin (3). 


First of all, we consider the case of disruptions which can appear in a switching circuit when the state of 
some relay X changes, 


It can be shown that homologous contacts (all make, or all break) cannot disrupt circuit operation during 
the transition periods and, therefore, analysis is required only for circuits in which make and break contacts of 
relay X enter simultaneously, 


For the analysis, we use the method of elimination of variables [6] by expressing the circuit f (x, y,..., Z) 
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Fig. 1. 
in terms of the contacts of X [1, 2): 


f(z, 2. y,..-,2) =2f/(A, 0, y,...,2)+2/(0, 1, y,...,2). 


(1) 
Since we are only interested in contact x, we shorten expression (1) to the following form: 
f(x, 2) = xf (1. 0) + z/(0. 4), (2) 
where f (1, 0) and f (0, 1) are functions which include all remaining variables but x, 
It can be shown that 
7(1. 0)7(0. 1) = /min (3) 


is an exact lower bound of the function f (x, x) since it contains only those subconstituents [2] (i,¢., constituents 


of the circuit which consists of contacts of relays y,..., Z) which enter into the circuit simultaneously in con- 
junction with x and with x. On the other hand, 


f(A. 0) + /(0. 1) = fmax (4) 
is an exact upper bound, since it contains all the subconstituents entering into the circuit. 


In the transition period, as was already stated, there is a severing of the mechanical connection between 
make and break contacts of one relay, X. 


It is completely obvious that no disruptions in operation will be observed if, for all possible combinations 
of states of the individual contacts of relay X in the circuit, no new circuits appear and no existing circuits are 


disrupted, and this corresponds to the assertion that so new subconstituents must appear and no existing sub- 
constituents must be expunged, 


We now analyze the case when the circuit contains two nonhomologous contacts of relay X: make contact 
x and break contact x, which we shall consider as contacts of the individual relays X; and Xg. Thus, we have 
implemented the transition from the circuit f(x, x) to the circuit f(x}, xX). If the scheme is described by 


structural formulae or by the matrix of immediate conductions, then this replacement entails the transformation 
of these formulae or matrices, 


As was already stated, circuit disruption is possible in the case when both contacts tum out to be simultan- 
eously open or closed, 


In the first case, for explaining the possibilities of momentary open circuits, we substitute the values 
X, = X_ = 0 in the structural formula, If the expression thus obtained, f (0, 0), equals f,,jn, this indicates that 


none of the circuits will be disrupted, If one obtains the inequality® 


*The symbol f, = f;, or f = f;, (f includes f,) indicates that § contains all the constituents which occur in f, 
and also some quantity of additional constituents [6, 7]. 








' then this indicates that there can be disruptions in circuit 
operation, 


By analyzing the expressions for f(0, 0) and f, we can 
determine in which circuit states open circuit can occur, 
For this, we expand f (0, 0) and fin in terms of the sub- 
constituents (for all variables except x), and those subcon- 
stituents of f,,in which do not occur in f (0, 0) define the 
states in which the circuit is open during the periods of 
operation or release of relay X, Such “openable circuits 
we denote by g. By definition, it follows that 





& = Ienis O-0). (6) 


For finding the constituents of the “openable circuits" 
we can recommend the use of the coordinate method of 
writing Boolean functions [8], where one notes, by points for 
example, those quadrants which correspond to f,,,.4p, and 
then strikes out those quadrants which correspond to F (0, 0), 
The nonstricken quadrants with points correspond to those 
circuit states in which open circuits will appear, 
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In the second case, when both contacts are simultan- 


. b eously closed, i.e,, when x, = x, = 1, false circuits can 
appear if the following inequality is obtained: 
Fig. 3, 
f(4.4) > Tf prax = 1 (1-0) 4-7 (0. 4). (1) 


Those states in which false circuits can appear may be found by a method analogous to the one given pre- 
viously if one notes by points those quadrants corresponding to f(1, 1), and then strikes out those of them which 
correspond to frnax. 


The sum of the subconstituents of f(1, 1) which do not occur in fax we call the "false circuits," and 
denote by g, where 


= aay! (1: 4). (8) 


As an example, we consider the circuits shown in Fig. 1, for which we can write the structural formulae; 


fi (2, 2) = zb+ za, fy(x, 2) = (z+ a) (z+) +6) 
fa(x, 4) =2b+2a+ab, fa (x, 2) = (x + a)(z+b)(a+ d). 
For all the circuits we get 
Snas™ € (10) 
hnay = 4+ b. (11) 
By replacing x by xq and x by x2, we get 
fi(%, #2) = 2b + aa, Is (%1» Ty) = (41 + @) (Xe +b) 
Jo(X1, Le) = 2yb+ T2a-+- ab, Aa (4%, Ze) = (41 + a) (x2 + b) (a + p. (12) 


1376 








2) 








} 
+ 


Al 


— 


ae 
ee | 
dra 


a 
' 





L 





+ 


Fig. 4. 


*. 


iL. 


r? 


t. 


From whence, 


/,(0, 0) =9, A(i, 1)=a+p4 
/,(0, 0) = ab. f.(4, 1) =a+ b, 
/3(0, 0) =ab /s(1, 1) = 4, 

14(0, 0) = ab; fa(A, 1) =a + 6. 


(13) 


By comparing these expressions with (10) and (11), we can see 
that an open circuit is possible in the circuit of Fig, 1a and a false 
circuit can appear in the circuit of Fig. 1c. There will be no dis- 
ruption of operation during the transition periods of the circuits of 
Fig. 1b and d, 


By comparing f,,,4,, and f, (0, 0) we obtain g = ab and, from 
a comparison of f, (1, 1) and f,,4,5 We get g = ab, Le,, in the 
circuit of Fig, 1a an error is possible when relay A operates and 
relay B does not operate and, in the circuit of Fig, 1c, an error is 
possible when A does not operate and B does, 


We note that, generally speaking, it is not possible in practice 
to replace x and x in the structural formulae first by x, and by x, 
and thereafter by zeroes or ones, It is possible to substitute ones or 
zeroes instead of x and X directly in the structural formulae, 


For class N circuits such a substitution must be carried out 
in the structural matrix of immediate conductions (or in the 
corresponding determinant), after which one can expand the deter- 
minant, 


We consider, for example, the circuit shown in Fig. 2, for which 
the structural matrix will be 


z 


The conductance between poles 1 and 2 is expressed by the determinant 


From which we obtain 


10 a 
M= 01 bed (14) 
a br 1 ¢ 
ade i 
0 az } 
/u= ba i cl]. (15) 
d ci 
} (2, z) = acd + abx + dz, 
{(1, 0) =ab+acd, (0, 1) =d, 
min = 2d + acd, (16) 
f max = 2 + @. (17) 
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For the transition period, we get from (15) that 


0aQ0 Oai 
{(0, ON) =|O1c}], fi, 1p)—=lbt ec]. (18) 
dec 1 de i 
From whence 
1(0, 0) = acd < fin (19) 
f(i, 1)=d+ab+ be>f Ay (20) 


Consequently, in the transient periods in the given circuit it is possible that, in certain states, both false 
circuits and open circuits can occur. 


Figure 3a shows the finding of the open circuit for the given case 

g = abcd, (21) 
Figure 3b shows the finding of the false circuit 

g = abcd, (22) 


2, Measures to Prevent Circuit Disruptions in the Transition Periods 





If it is observed that the circuit disruptions in the transition periods can influence circuit operation as a 
whole, i,e., if they occur in states of the circuit which are used [2, 9], then it is necessary to take corresponding 
measures, Among these measures are: 


a) the use of sequence fixing of the actions of the individual contacts of a relay (transfer contacts, or 
successively acting contacts); 


b) changing the structure of the circuit so as to exclude nonhomologous contacts from the circuit; 


c) applying lags, covering the duration of the circuit disruptions due to transition processes, to the relays 
on which the given circuit acts; 


d) the introduction of covering circuits which would remain closed or open during operation or release 
time of the relay whose switched contacts are in the main circuit, 


The first three measures require no special explanation, so we shall consider only the conditions which 
must be satisfied by a covering circuit, 


For the case of momentary open circuits in the transition period, the covering circuit is connected in 
parallel to the circuit under investigation, and must include, as a minimum, the subconstituents of all states in 
which an open circuit is possible, and can also include all the remaining subconstituents entering into fppjin, Le. 
the shunting covering circuit (> must satisfy the inequality 


<6 </mir (23) 


The solution of this inequality can be written in the form [6] 





/ min 0, 0 


where w is an arbitrary circuit. 
In the simplest cases, the covering circuit can be found from the well-known equivalence 


ra + 2b = za + zb-+ ab, (25) 


to which corresponds, for example, the transition from the circuit of Fig, 1a to the circuit of Fig, 1b, 
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Fig. 6, 


In the second case, when false circuits appear, the covering circuit must annihilate them, For this, the 
covering circuit # must be connected in series with the circuit being investigated, and must satisfy the inequality 


: P>0>I max’ (26) 
1.¢., 
_ ° 14 
=f maxt 0p = fmax+ & = fmax t+ SS . (27) 


In the simplest case, one can use the equivalence 
(a + a) (x + b) = (x + a)(z + b) (a + d), (28) 


to which corresponds the transition for the circuit of Fig, 1c to the circuit of Fig. 1d. 
For the circuit of Fig, 2, we obtain from (16), (17), (21) and (22) 





8 gp 
§ = abed + 490 yd + ab+ 
a+b+c+d 
+ 0 


In this case, the simplest covering circuits are 
p = abd, p = d 4 a. 


The corresponding circuit is shown in Fig. 4a and, in Fig. 4b, the same circuit, transformed by the general 
methods of contact circuits, 


38. Dynamic Disruption 





If there are make and break contacts of some relay X in a circuit, where the total number of contacts of 
this relay in the circuit is greater than two (i,e,, there are several homologous contacts) then, together with the 
| disruptions of which we spoke above, additional pulses can appear in the circuit during the transition periods [4], 

Thus, for example, we consider the circuit of Fig. 5, where aa x and x* are the numbered contacts of relay X, 
In the case when relays A, B and C operate (D does not operate), the circuit is closed when relay X operates, 























sabi hin if the contacts so react that initially contact x! closes, then contact x" breaks and only then does con- 
tact x* make, the circuit will initially close, then open and then close again, i.e., an additional pulse appears, 
One may easily convince oneself that such a pulse does not appear for any other sequence of contact reactions 
when relay X operates, 


The disclosure of such disruptions can be done by the method described in this paper by considering each 
of the contacts of relay X as the contact of an independent ficticious relay and by investigating the possible se- 
quences in the action of these relays, It can be shown that there can be disruptions only in the cases when between 
the reactions of two homologous contacts (make contacts, for example) there react one or several contacts of the 
other type (break contacts), Thus, for the circuit of Fig. 5, “dynamic” disruptions are pean only for the se- 
quences x¢--;° or xe- x - x4, In other words, if the circuit is written in the form f (x', x’, x’) then, in the 
first case, the sequence of states f(0, 0, 0)— f(1, 0, 0)— £(1, 1, 0)— f(1, 1, 1) is chosen and, in the second case, 
the sequence f (0, 0, 0)— f(0, 0, 1)— £(0, 1, 1)— f(1, 1, 1). 


To disclose those states in which, for a given sequence of contact reactions, the appearance of "dynamic 
disruptions” is possible, one can use, for example, a table, constructed for each sequence, of circuit states 
corres mespona to - re subconstituents, Thus, Fig, 6 shows two such tables for the sequences of reactions 

%* and x*- x - x, Here, the blank spaces show that the circuit corresponding to the given constituents 

oo be open, and ri the plusses will be closed, at each cycle engendered by the making or breaking of contacts 
of relay X. The sequence ~ +— +(or +— +7) shows that, for the states of the circuit elements corresponding 
to the given constituent, an additional pulse is generated upon the closing (opening) of the circuit in the process 
of operating or releasing relay X. For the circuit of Fig, 5, abcd and abed are such states, In a number of other 
cases momentary closures and breaks are possible, but these are brought to light by an investigation carried out 
in accordance with section 1 of this paper. 


If, for some reason, it is impossible to exclude some of the contacts of relay X from the circuit when there 
are “dynamic” disruptions, one can use covering circuits, analogous to those previously discussed, 


4. Simultaneous Change of State of Several Relays 





If there are cycles in the circuit during which several relays are operated or released simultaneously then 
the contacts of these relays may react, not simultaneously, but with some dispersion in time, determined not 
only by the properties of the contacts themselves, but also by the times of operation or release of the individual 
relays, 


Analysis of such circuits and the disclosure of the possible appearances of disruptions can be carried out 
by the same methods as for contacts of one relay. 


In the cases where one observes the possibility of disruptions of operation due to the nonsimultaneous change 
of state of different relays, one may, in addition to the methods mentioned above, introduce a delay in the re- 
lays so that the sequence of their operations is fixed, or one may duplicate the contacts of one relay by the con- 
tacts of another relay. 
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ON ONE METHOD OF ANALYZING COMPLEX 
SWITCHING CIRCUITS* 


V. F. D* yachenko 


(Moscow) 


A method is presented for the analysis of switching circuits which is based on 
the finding of the intermediate relay windings and executive elements carrying current, 
and their interconnections, for each state of the circuit's contacts, 


Two-terminal, class P and N switching circuits with multiwinding relays and 
parametric dependencies are considered, 


1. Introduction 


The problem of analyzing switching circuits amounts to the determination, for nonsequential circuits, of 
the state of the executive circuit as a function of the states of the receiving elements and, for sequential circuits, 
the determination as well of the sequences of action of the intermediate and executive elements, For this 
analysis one must have, in addition to the structural scheme, data on the relays and other elements (resistances, 
number of turns, ampere~turns for operation and release, etc) which enter into the circuit, For multi-winding 
relays it is necessary to have data on the interactions of the individual windings, 


The existing analytical methods [1-5] are based on the finding of the contact circuits which act on each 
intermediate and executive element of the switching circuit, Thus, for example, for a class P switching circuit 
[(1], the structural formula 


FF az $(@;, Gs, ..-, Ge, Za, Zqr-+ +1 En, Ay, Xq,..., Xn), (1) 


where Xj is a symbol for the reacting organ of an intermediate or executive element (j = 1, 2,..., 2), a4, Agreees 4k 
are contacts of the receiving relays, x4, Xgree0» X,, are contacts of the intermediate relays of fj = fj (Ag, Agrees Ay, 
Xq» Xgreces x) is the contact circuit acting on element Xj» must be transformed to the form 


E=f,X,+foXe+ +++ + f/nXn, (2) 
i,¢,, must go over to the group of circuits of the form £jX;.- 


For complex circuits, however, it is necessary in many cases to have recourse to additional graphical trans- 
formation of the circuits prior to separating out, by the existing methods, the contact circuits acting on each 
intermediate or executive element, For nonplanar circuits, and sometimes for circuits with parametric depend- 
encies, it turns out to be necessary to use the existing methods to isolate the circuits acting on the individual 
elements of the circuit, The analysis is still further complicated for circuits with multi-winding relays, For 
example, the general form of the structural formula for a class P circuit with multi- winding relays is written as: 


F = f(a, Go,.-+, 2k, Z, Zq,+++sTns Xi, Bac 0k Oe ps AS x. <a (1a) 


*Presented at a joint seminar of the Laboratory for Wire Communications of the AN SSSR and the Telephony 
Dept. of the MEIC on January 15, 1958, 
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In individual cases one can go from this formula to the formula 


F = fay Xi-+ fag Xi +--+ fea Xi + fog KSA + + fy Xb t- faa XS + vos (2a) 


where Xj (j = 1, 2,...,n; i121, 2, 3,...) is the symbol for the ith winding of element i and fi = fy (Ag, Agrooes ays 
Xq» Xgeeee» Xp) is the contact circuit acting on element x}. 


In this case, the action of multi-winding relay X; depends not only on whether or not the circuit (fjj) for 
some winding (Xj) is open or closed, but also on the interaction of the magnetic fluxes induced by the individual 
windings. In such cases, one should have recourse to the electrical computations of the currents flowing through 
the windings and of the magnetic fluxes induced by them. 


2. Amalysis of Nonsequential Class P Circuits 





In the present work we consider two-terminal switching circuits, In order to determine whether or not an 
intermediate relay or executive element of the circuit operates, we must find, for each state of the circuit's 
contacts, the windings of the intermediate relays and executive elements through which current flows, and then 
determine the interconnection of the windings in the circuit. 


For the analysis of class P switching circuits, we write the circuit in the form of structural formulae (1), 
(la). With this, we adopt the convention that the symbol for a winding in which the current flows from end to 
beginning is written with a tilde (for example, Xg) while the symbol for a winding in which the current flows 
from: beginning to end has no tilde (for example, x}, Thereafter, instead of the contact symbols, we substitute 
the values of their structural conductances (1 for a make contact and 0 for a break contact) for a given state of 
the elements, and we transform the expression obtained, bearing in mind that 


0.W =0, 
04W=W, 

1.W=W. @) 
14W=1, 


where W is the symbol for a relay winding or an impedance, After this, the structural formulae can assume one 
of the following forms, 


1, F=0, This says that the entire circuit between terminals of the battery is open, so that no current 
flows through any circuit element, 


2, F=1, This means that terminals are short-circuited, and that no current flows through the relay wind- 
ings. 


3. F=f cxi, a xi, | “s. Xe node R). This asserts that the relay windings and impedances which 


correspond to the symbols actually in the formula carry current, where the sumbols +, and brackets (or parentheses) 
joining the symbols of these elements indicate that these elements are interconnected, By knowing the relation- 
ships of the parameters of these elements, we can determine which relays will operate for a given state of the 
circuit’s contacts, For this, it is necessary to take into account, for multi-winding relays, the total action of all 
the current-bearing windings, A multi-winding relay will operate (given the condition that each winding induces 
the same number of ampere turns) if unequal numbers of fundamental windings (where the current flows from 
beginning to end) and counteracting windings (where the current flows from end to beginning ) are carrying 
current, and will not operate if equal numbers of basic and counteracting windings are carrying current, even 
though the relay as a whole is carrying current. 


If the interconnections between the circuit elements are such that the switching in of one element com- 
pletely upsets the operation of another, i.e,, the conditions of order of conductance are satisfied,* we may then 


*The concept of order of conductance characterizes the action of some circuit element (an impedance, a relay 
winding), with conductance G, on the operation of a relay A. It is considered that, if conductance G is connected 
in series with relay A and so decreases the winding current that the relay neither operates nor holds but, with the 
two in parallel, relay operation is not disturbed, then the order of conductance of G is less than the order of con- 
ductance of A (G< ,) The concept of a higher order of conductance (G , ) is introduced in an analogous fashion, 
(Footnote continued on following page) 
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use the additional transformations for analysis of the switching circuit: 


































Acb -|-Aza = B >A» (4) 
Acb Bua = Acp- 


If the interconnections between circuit elements are otherwise, and do not satisfy the conditions of order 
of conductance, we can then introduce other rules, Thus, for example, if impedance R, when connected in series 
with relay A, results in the latter not operating, but holding, then impedance R can be replaced by a contact, a, 
of relay A, If a parallel connection of relay A and impedance R causes A, not to operate, but to hold, then im- 
pedance R can be replaced by a contact a, 


3. Analysis of Class N Circuits 





In the analysis of a class N switching circuit [1], its structure is written in the form of a structural matrix 
(8, 9) 


M =||o,;| @ (=1, 2, 3,..., 9) (5) 


where oj,j is the directed path of current flow from node i to node j, with the convention that a; ;=1, The 
relay winding symbols which enter into the expressions Oj must include a tilde if the beginning of the winding 
is connected to node j, and must be without a tilde if the "connection is to node i, For multi-winding relays, 
each winding must necessarily bear its ordinal number, which is written in the form of a superscript, 


If there are no rectifier elements in the circuit, the matrix will be symmetric with respect to the principal 
diagonal, except that the winding symbols on apposite sides of the diagonal will be peeior WRN by their bear- 
ing a tilde or not, ie,, if, for example, a; ; = ax, then oj i= = ax!, and if oj, ; =aY + x?, then a i aY + X?, 


For the analysis of switching circuits with diode elements, we can consider that the impedance of a diode 
element in the conducting direction equals 0, and in the cut-off direction, equals @, Therefore, in setting up 
the structural matrix, a diode element is replaced by 1 if the direction of the current flow path coincides with 
the diode’s conducting direction, and by a 0 for the reverse direction. In this case, the matrix will be non- 
symmetrical, 


For example, for the circuit of Fig. 1, the matrix is written in the form: 





1 0 aX? b-- dR, 
M aa 0 1 Z R, + dy ¥ (6) 
aX* Z 1 x* 
b+-dR, R,-+ dY X? 1 | 
For the circuit of Fig. 2, the matrix assumes the form: 
As abZ aB} i 
abZ 1 bR A} (7) 
aB! bR 4 (1 + B) A? 


i Al A*(B*+ 0) 4 


To find, from the matrix, all possible current-conduction paths in a circuit, it is necessary to strike out 





*(Footnote continued], by being applied to the case when a parallel connection of G and A causes relay A to 
be shunted and its operation disturbed, In the case when neither a parallel nor a series connection of conduct- 
ance G and relay A affect the latter's operation, one considers that they have identical orders of conductance 

(G2,a) For more details on this, cf., [6, 7]. 


In theory, the determination of the orders of conductance of circuit elements (impedances, relay windings) 
can be carried out prior to circuit structural analysis, and independently of it. 
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the column whose ordinal number coincides with the number of the node con- 
nected to the battery’s plus terminal, and the row whose ordinal number coin- 
cides with the number of the node joined to the battery’s minus termimal, (Such 
nodes are sometimes called the circuit's terminals), We expand the determinant 
obtained by the rules for expanding Boolean determinants [8, 9], Le., we take 
all terms with a + sign, and take into account that, for contacts, a-a = a and 
a-a=0, The expansion of Boolean determinants gives, not only the symbolic 
expressions for the nonself- intersecting current-conducting paths, but also the 
different redundant terms, which correspond to sequences of nodes which do not 
enter into the current-conducting paths, These terms are easily removable, In 
; addition, in expanding a determinant which contains relay windings and ele- 
Vy. ments of finite conductance, one must take the following transformations into 


‘2 se ; E account: 














4 
i i 
. Xj-X; =0, (8) 
. F,+ F,-F, = F,, (9) 
at where F,; and F, are arbitrary two-terminal circuits (in which, in addition to ele~- 
Z ments with finite conductances, relay contacts may also enter), 








By thus expanding the determinant, we obtain an expression whose terms 
represent the symbolic expressions of the elementary nonself-intersecting current- 
conducting paths in the circuit between the + and— terminals of the battery, 
Further transformations are carried out in the same way as for the structural 
formulae of class P circuits, After substitution of the values of the structural conductances of the contacts for a 
given state of the circuit in the formula, one must transform the resulting expression so that, for each winding 
and each resistance, there would remain just one symbol in the formula, after which one can analyze relay 
operation as was done above, If it is impossible to lead the formula to a form where there is just one symbol 
for each winding, this indicates that the elements are connected in a bridge circuit, In this case, additional 
quantitative investigation is necessary, 


ok 2. 


We consider the switching circuit shown in Fig, 1, where it is known of the elements that all relay windings 
have the same order of conductance and, moreover Ry ¢ yxix2z- The matrix for the immediate conductances 
of this circuit was given above, (6), By striking from this matrix the first column, corresponding to the battery's 
plus terminal, and the second row, corresponding to the battery’s minus terminal, we obtain a determinant whose 
expansion gives the elementary current-conducting paths from battery + to battery— : 


0 aX! b+dR, 
F (a, b, d) = Z 1 X? |= 
R,+dyY X? 1 


= aX! X? R, + adX! X*Y + bZX* + dR, ZX* + 
+ bR, + bdY + dR,R, + aX%Z. (10) 


We now determine the state of the circuit's relays if, for example, receiving relays A and D operate, and 
relay B does not operate: 


F(i, 0, 4) = X'X°R, + X1N*Y 4+ XZ = 
= X1{X9(R,, + Y) + Z) = X*[X*Y + Z}. 


From this expression we can draw the conclusion that relays X, Y and Z will operate, If only receiving 
relay A operates, we get 


F(1, 0, 0) = X1X*R, + R,ZX* +- RR, + XZ = 
= X1(X*R, + Z) + R(ZX* 4- R). 



















It is impossible toso transform this expression that there would be only one 
symbol for winding X* and, consequently, this winding is connected in a bridge 
circuit. 





4, Analysis of Sequential Circuits 








For the analysis of sequential circuits, it becomes necessary to establish the 
temporal sequence of circuit element actions, i,e., to obtain a switching table, 
For this, one first of all writes the circuit*s structure in the form of structural 
formulae (1), (1a) (for class P circuits) or in the form of a structural matrix as in 
Fig. 3. (5) (for class N circuits), It is also necessary to know which circuit elements are 
receiving, and which intermediate and executive. 


The analysis of sequential circuits is carried out analogously to the cases 
already considered, with the sole difference that the element states are now taken 
in the same order as they occur in the process of circuit operation, With this, if it 
be assumed that intermediate relay X; operated (x; = 1) in some state and, upon 
substitution of the contacts® structural conductances for their symbols in the structural 
formula, it turned out that the symbol for this relay’s winding did not occur in the 
right member, this would indicate that the given state (given cycle of the se- 
quence) was unstable, and that it was necessary to go to the following cycle, in 
which relay X, must not operate, Analogously, if it were assumed that Xj did not 
operate (x; = 0), but the conditions for its operation were created, then it must 
operate at the following cycle. In the stable cycles there must be a correspond- 
ence between the states of the relays and the current flow through their windings, 
In this case, a transition to the following cycle must occur as the result of a change 
in an external stimulus, On the basis of such an investigation one constructs a switching table, If the sequence 
of € xternal stimuli is unknown (the operation of the receiving elements), it is recommended that analysis of the 
sequential circuit operation begin with the state when all the receiving elements are nonoperated, This will be 
the zero cycle, which may be either stable or unstable, 











Fig. 4. 


In conclusion, we consider three examples of the analysis of the double pulse divider circuits shown on 
Figs. 2, 3 and 4, where I is the receiving (pulsed) relay, A and B are intermediate relays, and Z is the executive 
relay, 


In the circuit of Fig. 2, the relay parameters and the resistance are so chosen that parallel or series con- 
néction of the windings of relays A, B, Z and resistance R does not disturb the operation of these relays, Relay 
A has two windings and, when current flows through both windings in different directions, relay A releases, 





By striking out the first column and second row in the matrix, (7), for this circuit, we get 


abZ ab’ i 
F=1|bR 4 A? (11) 
A! A*B? 4 


By expanding this determinant, we obtain 
F (i, a, b) = abZ 4-iA'-+ abB'R 4- aA! A* B'-|- ibs? BPR = 
= abZ +- i(A! + bA® B2R) 4- a (bB'R +. A* A* B}), (12) 
We assume that, in the initial state, all relays were nonoperated, i.e, 


F,(0, 0, 0) = 0. 


This will be a stable state of the circuit, To set the circuit in action, it is necessary that there be an 
external stimulus, i,e,, it is necessary that receiving relay I operate, After this occurs, we shall have, in the 
first cycle 


F(A, 0, 0) - A}, 
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This state is unstable, since there is no correspondence between the relay state (a = 0) and the flow of current 
through its winding, Relay A operates, and the circuit transfers to the following stable state: 


F,(1, 1, 0) = A! + A! A? B! = A}, 
To transfer to the following state, receiving relay 1 must release, Then we get 
F,(0, 1, 0) = A’ A? Bt. 
In this unstable state, relays A and B operate, with which 
F,(0, 1, 1)=Z24+ BR+AAB=Z-+ B(R-4- AA). 


At this cycle we have a source for executive relay Z but the circuit state will be unstable, With the 
application of an external stimuli, receiving relay I again operates, with which 


F,(1, 1, 1) = Z4- A’ + A® B*R + BR + Al AB! = Z + Al+ (A*B? + BY R. 


It is clear from thts expression that the currents in the windings of relay A will flow in opposite direction, 
so relay A must release, i,e,, the circuit transfers to the following stable state: 
F,(1, 0, 1) = Al 4- A*BeR. 
After release of relay I we get 


F,(0, 0, 1) =0. 


Relay B releases, and the circuit returns to its initial state: 
F,(0, 0, 0) = 0. 


If an external stimulus is again applied to the receiving relay, the sequence of circuit operation is re~ 
peated, On the basis of the analysis just given, we can construct the switching table (Table 1), 


In the circuit of Fig, 3, the elements’ parameters are so chosen that the order of conductance of windings 
A? and B* is less, and the order of conductance of winding B* is greater, than the order of conductance of winding 
Al, i.e., (A7B*) < ai and Bt > Ale Since this is a class P circuit, we immediately write its structural formula: 


r (i, a. b) [i +a (A? }3*) -A'} (DB' Ate} A’). 








. TABLE 2 
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F, (0, 1, 0) =(A2B%). 41: A= A2B?, 


F(A, i, i)= Boas + Al= B', 
F, (1, 0, 1) = B'uas+ Al= B', 
F, (0, 0, 1) = 0, 
F, (0, 0, 0) = 0. 
The circuit has returned to the initial state, On the basis of these results, we construct the switching 











] 
A ~|/+1+/4+]-4 +{/+]+ 
B >~i+i+i + —>|i+ 
Al P © 4 © . 
At Om MO) © 
Bi e © © © 
B* - 10 
TABLE 3 
oO; 1 Ti eT e) eT 10 | 11 
T —EE —S 
A ~itit+i+ic +iricr 
B >i +ti+tit+ ~~ it 
Z ~|+/ + 
By reasoning as before, we obtain 

F, (0, 0, 0) = 0, 

F,(1, 0, 0) = A?, 

F,(1, 1,0) = A’, 


F,(0. 1, 1) = (A®B%) <4: (BY as + A) = (A2B Ga Al = APB, 


In the circuit of Fig. 4, the relay parameters and the impedances are so chosen that parallel or series 
connections of windings A, B, Z and resistors Ry, Re and Rg does not disturb the operation of these relays, 


The structural matrix for this circuit has the form: 


*Note, For clarity, we have written, in the lower portions of Tables 1 and 2, the relay windings, on the left, 
and, in the columns corresponding to the cycles of circuit operation, we have noted by dots those windings 
carrying current at the given cycle and, by circled dots, the windings through which current sufficient for relay 




































We strike out the first column and second row of matrix (13), to obtain the determinant: 








abZ ia0d0 
O1ab5 
PeatO «et as (14) 
R, b A1 0 
R,b BO 1 
By expanding this determinant, we get 
F (i, a, b) = abZ + aAR,+ aBR,+ iaAR,+ iaBR,+ 
+ ibR,+ ibR,AB + ibR,AB + ibR, (15) 


Reasoning as before, we obtain 
F,(0, 0, 0) =0 
F, (1, 0, 0) = AR,+ Ry 
F, (1, 1, 0) = AR, +R, 
F;(0, 1, 0) = AR, + BR,, 
F,(0, 4, 1) = Z+ AR,+ BR,, 
F, (4, 1, 1) =Z + Ri+ BR,, 
F,(1, 0, 1) = R,+ BR,, 
F, (0, 0, 1) = 0, 
F, (0, 0, 0) = 0. 


The circuit has returned to the initial state, Table 3 was constructed on the basis of the analysis given, 


The method considered for analyzing mixed switching circuits broadens our capabilities of analyzing the 
operation of a more general class of circuits which are widely used in contemporary devices of automation, re- 
mote control and Communications, 
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EXPERIMENTAL DETERMINATION OF POWER POSITIONING 
SERVO SYSTEM CHARACTERISTICS 


M. B, Tumarkin 
(Khar’ kov) 


A method is considered for the investigation and testing of positioning servo 
systems which permits the direct recording of the systems’ velocity characteristics, 
The distinguishing feature of this method is the application to the system's input of 
a forcing stimulus with constant acceleration, 


Power “positioning” servo systems have, in recent years, been widely used for the automation of many 
types of machines, Thus, for example, several thousand machine stands are being used today in which servo 
systems are employed for the mechanization and automation of the processing of irregularly shaped surfaces, as 
well as for other purposes, 


In a large proportion of the cases, such systems consist of an executive mechanism, a controlling organ 
and a feedback path shunting the two, In these systems the executive mechanism motion must have an aperiodic 
character with tractive forces of the order of hundreds of kilograms and more, with velocities of as much as 
several meters per minute, Under these conditions, the system must provide a high degree of accuracy in follow~ 
ing position or trahslation, where the size of the discrepancies, as a function of the system’s purpose, must lie 
approximately within the limits of 0,03 to 0.3 mm, 


A very cogent question today is that of methods of testing and investigating these systems, One seeks 
methods which would allow the basic characterizing parameters of the systems, necessary for their rational 
technological use, to be determined, 


Among the basic system characteristics should be placed its velocity characteristic, this being the depend- 
ence of the velocity v of the executive mechanism on the magnitude of the discrepancy h, The velocity char- 
acteristic of a system mirrors the actual features of the phenomena which occur in the system during its operation, 
and is capable of being used to estimate the following (tracking) error, the system's zone of insensitivity (dead 
zone) and the system's tractive properties (behavior of the system under loads), These data are sufficient to allow 
the designer and technician to judge the possibilities of using the system in concrete conditions of its employ- 
ment, 


From the velocity characteristic one can determine the system’s velocity gain, and also the form of the 
system's nonlinearities, these being necessary data for the design analysis of the system, 


The development of methods for taking off the velocity characteristics of power positioning servo systems 
is of considerable interest, Essentially, this question reduces to the choice of the character of the forcing stimulus 
(input function) to the system, y (t), which would allow one to obtain such characteristics by the simplest and 
most reasonable methods, 


From this point of view, let us make an estimate of several standard forms of forcing stimuli, 


The unit pulse stimulus, Generation of an instantaneous unit stimulus to the controlling element is not a 
difficult task, For this purpose one ordinarily uses a plunger which, on priming, compresses a spring after which, 











upon release of the index pin, it provides a very rapid trans- 
lation of the controlling organ, It is assumed that this action 
is implemented instantaneously (Fig. la). The actual stimulus 
(shown by the dotted line) differs somewhat from the the- 
oretical, 














The unit pulse method is one of the standard ones for 
the investigation of automatic control systems [1)}. It is used 
for taking off the characteristics of the transient response, 

i and is capable, to a definite degree, of indicating the system's 
behavior from the point of view of its stability. In parti- 
cular, the so-called "time constants" are determined, these 
entering in the form of coefficients in the differential 





Fig. 1. 


equations describing the system, 


For the investigation of the systems under consideration here, the unit pulse method is of limited appli- 
cability, and is ordinarily used for generating equal stimuli for the investigation of system stability (5). Taking 
off system characteristics by this method is an extremely laborious task which, in many cases, does not achieve 
its aim. However, investigation of system stability by the unit pulse method is not always successful, since the 
system might have an oscillatory mode, not at the moment when a new position is reached, but at some inter- 
mediate position for motion of the system with an acceleration which, which such experiments, might not be 
expressed, 


To illustrate this statement, we adduce the example of the investigation of a hydraulic servo system with 
a controlling valve of positive overlap, Investigation of such a system by the unit pulse method showed that it 
was stable, but when the system worked in the continuous tracking mode, unstable portions were observed, 


Therefore, for investigating the systems under considerarion here, this method is of only supplementary 
value, 


Rectangular-pulse forcing functions. A rectangular pulse is applied periodically to the system's controlling 
element (Fig, 1b). 





Such a method is employed in the investigation of linear (or near-linear) automatic control systems with 
slowly flowing processes [2]. 


For the systems considered here, this method is not applied; it is essentially a modification of the unit 
pulse method, and {fs of little avail in obtaining velocity characteristics of systems, 


Constant-velocity forcing functions, The stimulus y* = v = const, carrying out its translations with uniform 
velocity, is applied to the controlling element. When the forcing motion is implemented by a cam, it may 
easily be given a periodic character (Fig. 1c), With this, the graph of the stimulus consists of segments with 
derivatives of constant modulus and alternating sign, 





The greater part of the investigations of the servo systems under consideration here are carried out with 
such stimuli [3, 4, 6, 7]. 


In this case, one obtains the system's velocity characteristic by successively changing the forcing velocity. 
Each velocity, plus the corresponding discrepancy, is recorded, giving one point on the velocity characteristic, 


With such a method of investigation, the accurate taking off of the system's characteristic is a very tedious 
and arduous experiment, When the transition is made from one velocity to another, it must be verified that all 
the remaining system parameters (oil pressure, temperature, load, etc) are unchanged, which is not always easily 
accomplished, Moreover, such an investigation gives no picture of system behavior when it moves with variable 
speed which, in many cases, is particularly important, The velocity characteristic cannot be directly given, but 
must be obtained only after the experimental data is processed; this introduces a definite subjectivity into the 
results obtained, 


Thus, the method of investigating servo systems by forcing functions with constant velocity, although 
sufficiently simple, cannot be acknowledged as ideal, 
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Fig. 3. 








Sinusoidal forcing functions, The particularly widely used method of harmonic sinusoidal forcing functions 
allows one to obtain the amplitude and phase-amplitude frequency characteristics of the system, from which one 
can determine the system's transfer function plus a number of other characteristics, and from which one can also 
estimate the system's stability. 





In such a method of investigation, the harmonic stimulus y = a (1— cos wt) (Fig. 1d) is applied to the con- 
trolling element, 


The frequency method of investigating automatic control systems is widely known and widely employed 
[1]. The question arises as to the effectiveness of this method for investigating servo systems of the type con- 
sidered here, ' 


It is necessary to start from the purposes for which the investigations of these systems are intended, As a 
result of an investigation, one must be able to obtain the system characteristic which will enable one to establish 
the technical possibilities of the system, The investigation must explain the influence of various factors on this 
characteristic, its correspondence with theoretical deductions, etc, It has already been stated that the velocity 
characteristic is such a characteristic, giving, as it does, the speed of tracking motion as a function of the 
magnitude of the discrepancy (for a given load) . 


With a sinusoidal input, such a characteristic can only be obtained after the proper processing of the results 
of an investigation, This processing requires the expenditure of a significant amount of time, and is very tedious, 


In investigations of the systems considered here, the capabilities provided by this method (determination 
of cut-off frequency, regions of stability, etc) are not used, since the working velocity of executive motion does 














not exceed 50 to 70% of the maximal, and the motion itself has an 
aperiodic character, 








Constant-acceleration forcing functions, From what has al- 
ready been said, it follows that the enumerated types of forcing 
functions do not satisfy the aim of obtaining directly the velocity 
characteristics of the servo systems under consideration here, 








One solution of the problem is the use of a method suggested 
by the author, This method is based on the introduction in the 
system of a forcing function which varies with constant acceleration, 
y* = w = const, 


The controlling element has applied to it a stimulus of the 
type y = wt'/2, generated by a cam of the proper profile. It is 
advantageous to give the cam a symmetric profile, so that half of 
its rotation would correspond to a forcing function with positive 
acceleration (y” = + w) and the other half, to a negative acceleration 


(y" =— w) (Fig. le). 


With such a forcing function, one can immediately record 
the system’s velocity characteristic in the velocity range from 
v =0 tov = yy, for different directions of motion, 





The cam profile (Fig. 2) is described by the curve p = py + 
+ a9”, where @ is the angle of rotation, The velocity of the forcing 
Fig. 4. function will be v = p*w = 4man@, and the forcing function accelera- 
tion will be w =p” =8 nan’, where n is the number of cam rotations 
per unit time, 





If we take the length of travel of the forcing motion, corresponding to a half-rotation of the cam, as equal 
to L, then 
L a 


a=, v = — Lnd, w = 8Ln?. 
v 


By varying the number of rotations, n, of the forcing function, one varies the range of velocities of the 
forcing function, within which the velocity characteristic is taken off, There is, with this, a corresponding in- 
crease in the acceleration of the forcing function, which agrees with the actual conditions of system operation, 
since an increase in speed of the executive motion in a servo system ordinarily is accompanied by an increased 
acceleration, This is related to the requirement that accuracy of tracking (following) be maintained, 


If it is necessary to increase the velocity of the forcing function without a proportional increase in the 
acceleration, one changes the cam profile by varying the angle corresponding to the necessary, or chosen, travel, 
L, of the forcing motion, 


With a given profile and a known number of cam rotations, it is not necessary to record the forcing 
function, It suffices to note, on the oscillogram, the beginning of motion, when v = 0, and the termination, 
when v = Vmax. In this case, the time marking on the oscillogram gives, in the proper scale, the velocity of 
motion, The system's velocity characteristic will be fixed on the oscillogram by recording the magnitude of 
the error (discrepancy), The oscillogram will also give the system's dead zone, and an indication of the char- 
acter of its motion, whether aperiodic or oscillatory, 


Figure 3 gives a functional schematic of the equipment for investigation of systems by the method 
described here, 


The plane cam 3, with the profile considered above, acts on controlling element 2, which controls servo- 
motor 1, This latter translates carriage 5, on which is positioned mechanism 4 which transmits the rotation of 
cam 3, The magnitude of the discrepancy h is recorded by transducer 6-7, Cam rotation is noted by marker 8, 
Mechanism 4 regulates the number of cam rotations for obtaining a given speed range. 
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The method given here for investigating positional servo systems by obtaining the speed characteristics 
directly is an objective method, allowing one to estimate the systems’ behavior and operations in conditions 
which approximate to the actual ones, 


The oscillogram of a system's speed characteristic is shown on Fig. 4, 


The oscillogram shows the presence of a dead zone (of the order of 0,05 mm) in the system's char- 
acteristic, this dead zone being related to the operation of a separating valve in the system, The oscillogram 
also shows that system operation occurs with “hunting” as a result of insufficient system damping. 


The method presented here can be used, not only for laboratory investigations, but also under industrial 
conditions for testing and finishing of servo systems. 
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THE USE OF MAGNETIC AMPLIFIERS IN INTEGRATING 
CIRCUITS 


S. B. Negnevitskil 


(Moscow) 


Two basic circuits employing magnetic amplifiers as integrators are con- 
sidered, 


Approximate expressions are derived for the transfer functions of these ampli- 
fiers, both without the amplifier time constant being taken into account, and with 
it. 


The disadvantages and advantages of integrating magnetic amplifiers are given. 


Cases frequently arise in technology when it is necessary to integrate dc voltages, As is well known, 
passive RC and RL circuits can be used for this purpose, However, the capabilities of passive circuits are very 
limited, 


In practice, therefore, active integrating devices are ordinarily employed, namely, the so-called integrating 
amplifiers, The most widely-used of these is the integrating amplifier with a parallel feedback path, the block 
schematic of which is shown in Fig, 1. 


The expression for the transfer function of such an integrating amplifier is usually given with the assumptions 
that the input impedance Ry, and the amplifier’s natural time constant r, are equal, respectively, to infinity 
and zero, When vacuum-tube amplifiers are used, such assumptions are very close to reality. If magnetic 
amplifiers (MAs) are used, their comparatively small input impedance and their time constants should be taken 
into account, 


By initially taking only the input impedance into account, we obtain the following expression for the 
transfer function of the integrating amplifier, constructed in accordance with the circuit of Fig. 1: 


Vout ky | 
Ro 








Here, the amplifiers voltage gain is ky = Upyt /Up, the system's time constant is 


and the loop time constant is r = RC, 


By substituting Rg» oo in (1) and (2), we obtain the well-known expression for a yacuum-tube integrating 
amplifier, It is clear from (1) and (2) that taking the input impedance into account does not changes the form 
of the transfer function, but leads only to a decrease in the system's time constant and, by this, to an increase in 
the error of integration (for the same time of integration), 
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The voltage gain is not always the characteristic feature of magnetic amplifiers, It is advantageous to 
introduce for them the concepts of current gain 


(3) 


and magnetizing force gain 


The coefficients ky» kj and kp are interrelated by the expressions 


kyW in 
ee 











k=—, k&p= ° = 





With this, the time constant of the system, (2), can be written in the form 











Rot+k, 
T= R™ (6) 


Ordinarily, in the practical applications of MAs, kj > Ro. With this, 
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W(?)=RER1+T?' (7) 
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T= Re ° 
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(8) 


For Tp > 1, the system approximates to an ideal integrator with transfer function 
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Moreover, it is frequently the case that R » Ro», and then 






T =k,C. (10) 








It is interesting to note that the numerical value of the coefficient ky, with the dimensions of ohms, 
corresponds to the input impedance of the amplifier in a circuit with 100% negative, fixed, series-connected 
feedback, 











The special features of an MA allow one to implement an integrating MA with no galvanic connection 
between input and output, With this, one applies to one winding the input voltage or current, proportional to the 
quantity to be integrated and, by means of the second winding, one implements the negative feedback (Fig. 2). 
An integrating amplifier constructed in accordance with the circuit of Fig, 2 will be called an integrating 
amplifier with magnetic feedback; its transfer function is written in the form 


U 1 
W (7) = “ptt map ER a1) | 


t= RC, T = (Ro + Wokg) C. 


































































































ay’ With pr « 1, which ordinarily is the case, the expression for the 
A system's transfer function takes the same form as for the circuit with 
Lew 7 parallel feedback (1): 
u > in Y out 
ss hat - 1 
o sod W (p) = kp Tp 5T (13) 
For pT » 1, when 
Fig. 1, 
kp 
W (p)= sre (14) 
a . 
the system approximates to an ideal integrator, 
I i EN —— You In practice, the following relationship ordinarily holds 
Woky =k, > Ro, (15) 
vie, 2 and then expression (12), for the time constant T of an integrating 
& 


amplifier with magnetic feedback, coincides with (10); with this, one 
has in expression (15) the current gain by the feedback winding W». 


It should be mentioned that, in vacuum~tube integrating amplifiers, in order to provide negative feedback 
it is necessary that the polarity of the output voltage be the reverse of the polarity of the voltage at the ampli- 
fier’s input, This explains why, in practice, for changing the sign of the output voltage of a vacuum-tube 
amplifier, an additional special inverting amplifier is connected, In an integrating MA with magnetic feed- 
back, the output voltage may have any polarity necessary, since a negative feedback can be provided by the 
proper connection of the terminals of the feedback winding. 


This property is an advantage of integrating MAs with magnetic feedback, particularly for use in complex 
interconnected circuits, 


We now consider the approximate effect of the amplifier’s natural time constant r,, We assume that the 
amplifier consists of two links, an amplifying one and an inertial one, and has the transfer function 


u 
%, =TFpr, 
The transfer function of the integrating amplifier constructed in accordance with the circuit of Fig. 1, 
with r, taken into account, will have a comparatively complex expression; however, in practice, the relation- 


ships r> fr, and ky/(1+R/R,) > 1 usually hold, and then, approximately 


i 





k 
W (p) = ——- 


: , (16) 
‘+R d+7p) (1+ 72?) 

As should be expected, the presence of r, leads to the appearance of a second-order term in the deno- 
minator of the transfer function, 


One easily obtains from (16) the zone of possible integrating frequencies, For r, = 0, this zone lies 
between the limits of 1/T and o, and for r, # 0, the upper frequency will be bounded by the quantity ky/ ra. 
It is important to emphasize that, as is clear from (16), the effect of the time constant r, is to decrease pro- 
portionately the gain, This circumstance widens the possibilities of MA employment in integrating circuits, 


It should not be assumed that an integrating MA with magnetic feedback is best in all cases, Depending 
on the concrete conditions; internal impedance of the transducer, the possibility of connecting the transducer 
to a matched load, the limits of variation of the input voltage, the amplifier parameters, etc., it is advantageous 
to use one integrating circuit or another, 











The detailed make-up of both integrating MA circuits is not the problem of the present paper, As an 
example, we consider a case where the problem {s only solved by the use of a MA with magnetic feedback, 


Let there be a transducer with maximal emf of E;,, = 0,02 volt and intemal impedance of R, = 10 ohms, 
We assume that on the magnetic amplifier, with account taken of the least possible diameter of the winding 
conductor, one input winding can be mounted, with Wi, = 6350, Rg = 950 ohms, ky = 952500 ohms, or two input 
windings can be mounted, with Wi 4 = 300, Ry = 10 ohms, ky = 45000 ohms and Win 2™ 5000, Re = 750 ohms 
and kyg = 750000 ohms, 

The amplifier magnetizing force gain is kp = 150, 

The input-output characteristic of the MA is linear up to the limit of Upy; = 70 volts, 


It is required to obtain a system with time constant T = 200 seconds and maximal output voltage of Upy: = 
= 45 volts during an integration time of 1 = 200 seconds, 


We consider the possibilities of using an amplifier with parallel, and with magnetic, feedback, 


Integrating Amplifier with Parallel Feedback 





As the input winding, we use the winding with W,,, = 6350 and R, = 950 ohms, and we assume that the loop 
impedance R is also equal to 950 ohms, Since, in the given case, the condition that ky >> Rg holds we can, by 
using (8), determine the necessary magnitude of capacitance: C = 420 microfarads, 


With a step- function connection of the maximum input signal, the output voltage after an integration time 
of 200 seconds will, in correspondence with (9), equal 


t 


If a lower value of R is taken, a larger value of capacitance C will be required, while if a larger magnitude 
for R is chosen, the already small value of output voltage is decreased, 


Thus, the given MA, used as an integrating amplifier in accordance with the circuit of Fig. 1, can provide 
the required magnitude of T with a comparatively large capacitance C, but cannot, in the given conditions, pro- 
vide the necessary U oy. 


Integrating Amplifier with Magnetic Feedback 





In this case, the amplifier must have two input windings, As the input winding for the transducer we take 
the winding with W, = 300 and, as the feedback winding, W, = 5000, Since relationship (15) holds, the necessary 
value of capacitance can be found from (10): C = 266 microfarads, 


We now determine the input magnetizing force: 
E 
in 
Fn = R, + R, W, = 0.3aw. 


With a step- function connection of the input voltage, the output voltage after the integration time will, 
on the basis of (14), equal 


t 
UVour = An Wwe= 45y. 


It is clear from this example that, with the given concrete conditions and for one and the same input 
quantities, the use of the integrating amplifier with magnetic feedback allows one to use a capacitance smaller 
by a factor of 1,6, and to obtain an output voltage greater by a factor of 4,5 than in the case of parallel feed- 
back, 


Experimental work with integrating MAs shows that the assumptions made are valid, and gives a com- 
paratively good agreement between experimental data and results computed from the formulae in this paper, 









SUMMARY 





Magnetic amplifiers have a host of specific special features, thanks to which integrating MAs have their 
disadvantages and advantages in comparison with vacuum-tube amplifiers, One of the fundamental disadvantages 
of integrating MAs is the decrease in the band of possible integrating frequencies, since an MA has a com~- 
paratively large time constant. Moreover, in an integrating amplifier with parallel feedback, the system time 
constant is decreased by a factor of (1 +R/R,); this decrease will be significant in cases where an MA is used, 
since its input impedance R, is comparatively small, 





The advantages of an integrating MA are the greater reliability of operation, the lower expenditure of 
energy, the lower influence of the condenser bleeder resistance in the feedback circuit and of the impedance 
of the isolation circuit, since the comparatively low input impedance of the MA causes a correspondingly low 
value of impedance in the remaining circuits, In an integrating MA with magnetic feedback, the polarity of 
the output voltage is not rigidly related to the polarity of the input voltage. In such an amplifier, the input 
quantity may be applied from a current generator, while an integrating amplifier with parallel feedback is use- 
less in principle in such a case, The input impedance, in comparison with all other schemes, has a com- 
paratively low value and can be easily matched with the transducer’s internal impedance, By means of an inte- 
grating amplifier with parallel feedback, one can integrate the sum of several signals from different sources with- 
out interconnecting them galvanically, The absence of any galvanic connection between amplifier input and 
output simplifies its use in complicated interconnected circuits, 
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CHRONICLE 


SCIENTIFIC SEMINAR ON PNEUMOHYDRAULIC AUTOMATA 


A. I. Semikova 


The Third All-Union Conference~Seminar on Pneumohydraulic Automata was held May 18-19, 1959 at 
the Institute for Automation and Remote Control of the Academy of Science of the SSSR, 


The conference was attended by 308 people from 18 cities; a total of 98 organizations working in the do- 
main of pneumo-automata were represented, 


In opening the session, seminar director M A. Aizerman noted the important role played by the annual 
conferences in matters of organization of work and in becomingcognizant of new developments in the given 
domains of science, Very important also are the personal acquaintanceships made during the course of these 
seminars, 


The conference was opened by a plenary session at which several papers were presented, after which the 
work of the seminar was carried in two sections: pneumo-automata and hydro-automata; a concluding plenary 
session was open to general discussion, 


The first paper presented was that of Z, Ya, Beirakh on an electronic-pneumatic autoregulation system in 
which electrical measuring and summing devices were joined with a pneumatic executive mechanism which fs 
more reliable than an electromechanical one, As transducers and amplifiers, the electronic portion of the 
system control, by means of an electropneumatic relay, a piston pneumatic servomotor, 


The system permits the automation of production where there is a risk of fire or explosion, since all 
electrical elements and circuits can be removed to a significant distance from the object itself and from the 
controlling organs, 


V. D. Mironov (Moscow), in the paper, “Sampled~=data electronic-hydraulic regulators,” spoke of the 
electronic-hydraulic variants of mass-produced sampled-data electronic-mechanical control systems, The 
speaker provided descriptions of two variants of sampled-data electrohydraulic executive mechanisms, The 
first variant contains a unique power element, namely, a geared hydroturbine controlled by a flexible valve, 
The second variant (simplified) is a membrane executive mechanism, controlled by two electrovalves, 


The paper of A, A, Tale and T, K, Berends (Moscow), "New elements of digital pneumatic devices,” con- 
tained descriptions of developed system elements which permit the construction of complex pneumatic switching 
circuits and any finite automata, The pneumovalves and pneumorelays permit the implementation of the ele- 
mentary operations of instantaneous logic (negation, repetition, conjunction, etc) and any primitive (non- 
sequential ) circuit, Moreover, the one~cycle delay elements constructed from this same apparatus show the 
possibility, in conjunction with the instantaneous elements, of being used in the construction of nonprimitive 
(sequential) circuits, The continuous memory cell which has been developed has an independent value in its 
own right. 


Based on the elements developed, devices were constructed which have important industrial value: iso- 
drome regulators with automatic changeover, controlling color-detectors (chromographs), etc, 


E, S, Arzumanov (Kirovokan) spoke of multiple-projection pneumo-electrical and pneumatic recording 
devices, 











The Section on Pneumo-Automata 





The first speaker in the section of pneumo-automata was L, A. Zalmanzon (Moscow) who, with the paper, 
"On the characteristics of small-size intervalve chambers," discussed an investigation of the characteristics of 
chambers in which, due to the constriction of the air flow by the walls, the pressure along the chamber length fs 
significantly changed, The paper contained a generalization of the characteristics, and cited a method for de- 
signing certain types of chambers, There was mention of the possibilities of using small intervalve chambers as 
correctors or as transformers of one pressure as a function of another, 


1, V. Vaiser (Moscow), in the paper, “Analysis of the possibilities of conversion of pneumo-automata instru- 
ments to lower supply pressures," considered the possibility of converting pneumo~automata instruments to supply 
pressures of 100 mm of water, which would provide a number of advantages, both technological and economical, 


The paper provided the results of experimental testing and comparison of the operation of basic pneumo- 
system elements with supply pressures of 1 atmosphere and of 100 mm of water. Considered were * nozzle- flapper*® 
type relays, membrane blocks, dead-end and flow-through chambers and long lines, 


The results of the experimental investigation allowed the conclusions to be drawn that conversion to supply 
pressures of 100 mm of water would not impair the operation of pneumosystem elements and, in addition, would 
essentially simplify the instruments’ air supply. 


P, M, Shanturin (Moscow) provided information on the design development of pneumo~automata instru- 


ments at the design branch of the TsLA, particularly of the design development of pneumatic instrument elements 
made of plastic compositions, 


Plastic composition pieces, prepared by pressing, have a heightened sturdiness and are readily amenable to 
mechanical processing. Mass-production assembly of parts is greatly simplified, since testing and fitting are 
reduced, Instruments prepared of plastic compositions can be used in chemically active (corrosive) media, are 
frost-resistant and heat-resistant and, moreover, their preparation is economically advantageous. 


V. G. Ponomarenko (Khar*kov), in the paper, *The use of pneumatic control apparatus for automating the 
benzene~scrubbing branch of a coal-tar chemical factory," described a scheme for automating the benzene 
section and the operation of a connected-control device which, with quality transducers present, allows a quality- 
control system to be constructed and which permits the unifying of control of all parameters in one aggregate, 


In the communication on their joint work by N, L, Kharas, V. R. Anders and T. K, Berends (Moscow), there 
was a description of a pneumatic control block of a controlling chromograph designed for the continuous in- 
spection of the gas flux make-up directly at the technological stands, An industrial chromograph, supplied with 
a pneumatic control block, opens new possibilities for complex automation, The control block was constructed 
of pneumatic digital apparatus, It contains a timer (a command instrument ) and an executive device (two dis- 
junction circuits and a continuous memory cell), 


A. A, Tagaevskaya (Moscow), in her communication, spoke of a developed attachment to an extremal 
regulator for its use on objects in which the quantity to be optimized depends on the position of n controlling 
organs, To realize the attachment, one used elements from systems of pneumatic discrete logical devices, 


A. L Semikova (Moscow), in the communication, “Experimental investigation of the characteristics of jet 
elements of pneumo-automata,” revealed the results of an experimental investigation of the effect of changing 
elements’ geometric dimensions on their output characteristics, 


R. A. Auzan (Moscow) gave a paper on the subject, "Dyamic characteristics of pneumatic regulators and 
their adjustment on industrial objects." The paper presented the results of an experimental investigation of the 
dynamic characteristics of iass-produced pneumatic control apparatus, 


The data of this work allowed one to characterize and judge the dynamic quality of standard blocks — the 
unified system aggregate (USA), the small unified system aggregate (SUSA), the system of the Chermet Central 
Laboratory for Automation (CLA)— and to develop instructions for the adjustment of USA regulators for a large 
class of industrial objects, The results were illustrated on one of the standard USA blocks, Estimates were given 
of the speed of response of several old type transducers and compensated transducers, 


N. D. Lanin (Moscow), in his paper, spoke of pneumatic analog computers, The operational elements used 
in the machine do not have any limitations on integrating time, for all practical purposes, 


1402 





Circuits were given for a differentiating link and for a delay block. 


V. L Chernyshev spoke of work, carried out jointly with V. N. Dmitriev, on the investigation of a piston 
pneumatic executive mechanism with a constant load on the piston rod, The communication contained the 
presentation of a graphicoanalytic method for computing the speed of a piston pneumatic executive mechanism 
with a constant load on the piston rod, The speaker provided a comparison of experimental and computed data, 


V. M. Gorokhov (Khar*kov) gave a paper on the subject,"A piston pneumatic servo drive for general in- 
dustrial use,” The new design of the pneumatic drive differs from the well-known ones by its simplicity, com- 
pactness and high sensitivity, 


In a joint communication, V. V. Volkov, E. Yu. Gutnikov and M, P, Kogon (Sverdlovsk) spoke of the auto~ 
mation of long-travel pneumatic drives, 


In the communication of G, T. Berezovets on his joint work with post-graduate student Chzhou Tszin-Len 
(Moscow), the results were given of an experimental study of the characteristics of membrane pneumatic re- 
gulators, The experimentally-determined characteristics were given of the translation of rigid center and of the 
change in pressure drop as functions of changes in effective membrane area, as well as the results of a study of 
membrane operation for membranes prepared of different materials, An estimate was given of the effect of 
membrane characteristics on the quality of pneumatic device operation . 


A. L Kazeev (Kirovokan), in his communication, gave recommendations on methods of designing and 
choosing controlling valves, 


E, A. Andreeva (Moscow) reported on the results of an investigation of power and discharge characteristics 
of “nozzle- flapper” type elements for the flow of compressible viscous liquids, A comparison of experimental 
and calculated curves was given in the communication, 


In a joint communication, A. G, Arkad’ev and M, L, Podgoetskii (Moscow) spoke of pneumatic transformers, 
transforming small forces of 0,2 to 10 grams to standard compressed-air pressures of 0,2 to 1 atmosphere, and of 
the development of a pneumatic transducer with an isodrome amplifier in the feedback line, The transducer has 
a dynamic tuning element, 


The Section of Hydro-Automata 





V. N. Veller (Moscow), in the paper, “Theory and design of high-sensitivity piston mechanisms,” presented 
a solution to the problem of increasing the sensitivity of piston mechanisms by the use of self-centering pistons, 
By means of special flow rings and drainage flanges on the piston’s working surface, the piston is maintained 
strictly coaxial with the cylinder during operation. 


L Yu, Klugman (Odessa), in the paper, “Certain considerations in the choice of hydro-copying machine 
systems,” considered hydro-copying machines and their basic elements, namely, valves, executive mechanisms 
and energy sources (pumps). 


The quality of a system is defined by the ratio of the cylinder diameters of the given system to the dia~ 


meter of the cylinder of the system of minimal dimensions, Recommendations on the choice of hydro-copying 
systems were given, 


Yu, E, Zakharov, in the communication, “On hydrodynamic force and the discharge coefficient in valves,* 
spoke of computing the hydrodynamic force and the discharge coefficient by the method of “jet flow” theory. 


A, F. Arkhangel’skii spoke of a developed small-size pump with high productivity of high pressures for the 
hydrosystem of a 250 hp tractor. 


In his address, M, B, Tumarkin presented a uniform approach to the study of hydraulic systems by three 
steps; 


1) system consideration is carried out from the points is view both of its design and of the control of the 
executive mechanisms; 























2) the system is judged with respect to six basic characteristics (speed, traction, speed gain, error, coef- 
ficient of rigidity and coefficient of insensitivity); 


3) system testing is carried out by taking down the characteristic in the form of the complete curve when 
a signal with constant acceleration is applied, 


In his address, N. S. Gamynin gave the characteristics of hydraulic aggregates used in servo systems, The 
speaker provided the velocity, frequency and load characteristics, as well as the transfer functions, of hydro- 
system elements (valves, pump drives) and of the systems as wholes, 


The paper of V, M, Dvoretskii (Moscow) was devoted to a device for obtaining small stable discharges of 
fluids, the principle of operation being control by the pressure drop on a constant valve, 


V. D. Konstantinov, in the communication, *A hydro-aggregate as the object of synchronous generator 
rotational speed control,” considered a hydro-aggregate consisting of a drive and a synchronous generator, gave 
various designs, deduced the differential equations and estimated the lag time. 


In the concluding plenary session, G. I, Papernii (Leningrad) presented a paper, “A pneumatic system for 
the automatic control of combustion and supply processes in boilers," in which the speaker dealt with a pneumatic 
system, developed in the TsKTI, for the automatic control of the combustion and supply processes in boilers, The 


development of the pneumosystem was based on the principle of force compensation, as being the most pro- 
mising, 


Industrial testing under actual conditions of use of the pneumatic automatic control system demonstrated 
the efficiency of all its elements, as well as their high sensitivity. 


The paper of V. L. Korobochkin (Moscow) contained the results of an investigation into the dynamics of 
hydraulic drives of copying machines, 


In conclusion, there was an address by V. S. Prusenko (Moscow) with some examples of the construction of 
scheme of cascaded and interconnected control of thermal processes based on the instruments of a pneumatic 
unified system aggregate designed for general industrial use, 


The comrades who took part in the discussion were unanimous in noting the great value of holding annual 
conferences, a number of remarks and desires were expressed for better organization of the exchange of ex- 
perience among the organizations taking part in the work of the conference, 


Several comrades, namely, P, M, Shanturin (Moscow), A. F, Arkhangel’skii (Chelyabinsk) and L. A. 
Sopochkin (Lisichansk), expressed their desires that the organization of the annual conferences be improved, that 
the publication of the Proceedings (Trudy) of the conference be accelerated, that there be a preliminary distribution 
of the subjects of the papers, that exhibits be organized, that there be a special bulletin of our seminar. 


Participants I, F. Kozlov (Moscow), V. G. Ponomarenko (Khar*kov) and others also expressed the desire for 
a broader enlistment in the work of the annual conferences of workers in factories and related branches of in- 
dustry, on whom depends the realization of the designs developed, 


The chief designer of the “Tizpribor* factory, M, I, Zhutovskii, noted the parallelism in the work of dif- 
ferent organizations, and expressed his opinion on the necessity for the coordination of work, The work of 
coordination should be headed up, in his opinion, by the IAT AN SSSR and TsNIIKA, 


A. P. Kopelovich stated the opinion that it is necessary to coordinate, with the consumer, the technical 
conditions on the instruments developed up to the time they are mass-produced, 


Z. Ya, Beirakhom also expressed the desire to bring KB“ closer to the factories, or to strengthen the KB at 
the factories by changing the system of payment. 


B. F, Stupak mentioned the necessity of organizing the issuance of catalogs by factories-producers. 


It was proposed, moreover, that the last two pages of the journal "Priborostroenie” be reserved for a dif- 
ferent type of advertisement and informational material. 


L, L. Feigel’son, G, L, Papernyi, Yu, E, Zakharov and others then spoke of their work and noted the in- 
adequacies of the currently existing information on current work and results obtained, In connection with this, 
the great value of these annual conferences was once again asserted, 


In conclusion, the conference tumed to letters to the TsK KPSS on questions related to the present state 
and prospective development of pneumohydraulic automata, 
*Publ, Note: KB, Design Office. 
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